第1章:量子力学基础——薛定谔方程、波函数与能带、DFT简介
各位同学好,我是你们这门课的主讲。今天咱们聊聊量子力学在材料模拟里的应用。说实话,我当年刚接触这个领域时,也觉得量子力学太抽象了——什么波函数、能带,听着就像天书。但做了十几年项目后,我发现这些东西其实没那么可怕。你想想看,我们做材料模拟,本质上就是在解薛定谔方程。
1.1 薛定谔方程——材料模拟的起点
薛定谔方程,说白了就是描述微观粒子运动的基本方程。我习惯把它理解成「量子版的牛顿第二定律」。在经典力学里,我们算一个球的轨迹;在量子世界,我们算的是粒子出现的概率。
定态薛定谔方程长这样:
Ĥψ = Eψ
其中 Ĥ 是哈密顿算符,ψ 是波函数,E 是能量本征值。嗯,这里要注意:这个方程不是推导出来的,它是薛定谔在1926年「猜」出来的——但后来被无数实验验证了。
核心要点:薛定谔方程的解给出了体系的能量和波函数。有了波函数,原则上我们可以算出材料的任何性质——电导率、光学性质、力学性能等等。
我在项目中遇到过一个问题:用经典分子动力学算不准半导体材料的带隙。后来改用DFT(密度泛函理论),结果就对了。为什么?因为经典方法忽略了量子效应,而DFT直接解薛定谔方程。
1.2 波函数——它到底是个什么东西?
波函数 ψ(r) 本身没有物理意义,但它的模平方 |ψ(r)|² 代表在位置 r 找到电子的概率密度。我刚开始学的时候总纠结「波函数到底是什么」,后来导师跟我说:你就把它当成一个数学工具,用来算可观测量的。
波函数有几个重要性质:
- 归一化:全空间找到电子的概率总和为1
- 连续性:波函数及其一阶导数连续(除了某些特殊点)
- 单值性:每个位置只有一个确定的波函数值
我的经验:做DFT计算时,波函数的展开基组选择很关键。平面波基组适合周期性体系,而局域轨道基组适合分子体系。我曾经因为选错基组,算出来的能带结构完全不对——浪费了两周时间。
1.3 能带——固体中的电子能级
在单个原子中,电子能级是分立的。但在晶体中,原子周期性排列,电子能级展宽成能带。为什么会这样?你想想看,两个原子靠近时,它们的电子轨道会重叠,形成成键和反键轨道。当N个原子排成晶体,就形成了N个紧密排列的能级——这就是能带。
能带结构决定了材料的导电性:
| 材料类型 | 能带特征 | 典型例子 |
|---|---|---|
| 导体 | 价带和导带重叠 | 铜、铝 |
| 半导体 | 带隙 0.1-3 eV | 硅、砷化镓 |
| 绝缘体 | 带隙 > 3 eV | 金刚石、氧化铝 |
避坑指南:我曾经用LDA泛函算氧化锌的带隙,结果只有0.7 eV,而实验值是3.4 eV。后来改用杂化泛函HSE06,才得到合理结果。记住:标准DFT泛函会严重低估带隙,这是系统误差。
1.4 密度泛函理论(DFT)——从波函数到电子密度
DFT的核心思想很巧妙:用电子密度 ρ(r) 代替波函数 ψ(r₁, r₂, ..., rN) 作为基本变量。为什么这么做?因为波函数有3N个变量(N是电子数),而电子密度只有3个变量——计算量大大降低。
DFT的理论基础是Hohenberg-Kohn定理:
- 体系的基态能量是电子密度的唯一泛函
- 使能量泛函取最小值的密度就是真实的基态电子密度
实际计算中,我们用的是Kohn-Sham方程:
[-ħ²/2m ∇² + V_ext(r) + V_H(r) + V_xc(r)] ψ_i(r) = ε_i ψ_i(r)
其中 V_xc 是交换关联泛函——这是DFT中唯一需要近似的部分。我常用的泛函有:
- LDA:简单快速,适合金属和简单半导体
- GGA:精度更好,适合大多数材料
- 杂化泛函:精度高但计算量大,适合带隙计算
实用建议:做DFT计算时,先跑一个简单的LDA计算看看趋势,再用GGA或杂化泛函精修。我一般用这个策略:LDA扫参数 → GGA算性质 → 杂化泛函验证关键结果。
1.5 本章知识体系
下面这张图展示了本章的核心逻辑:从薛定谔方程出发,通过波函数描述电子状态,进而得到能带结构,最后用DFT实现高效计算。
这张图你看懂了吗?从上到下,我们一步步从量子力学基础走到了实际材料模拟。DFT就是连接理论和应用的桥梁。
1.6 本章小结
好了,咱们把这一章的核心内容捋一捋:
- 薛定谔方程是材料模拟的理论基础,它描述了微观粒子的量子行为
- 波函数本身不可观测,但它的模平方给出了概率密度
- 能带结构决定了材料的导电性,DFT计算中要注意带隙低估问题
- DFT用电子密度代替波函数,大幅降低了计算量,是现代材料模拟的基石
我个人觉得,学DFT最好的方法就是动手算。下一章我们会讲怎么搭建计算模型、选择参数。到时候我会分享一些实战中的小技巧——比如怎么处理收敛问题、怎么判断结果是否合理。
今天就到这里。有问题随时问我,咱们下次见。
公众号:蓝海资料掘金营,微信deep3321