课程信息

本讲是”第一性原理计算方法”(Lectures on First-Principles Calculations)的引言,授课教师为清华大学物理系王冲(Chong Wang)。

  • 主教材:R. M. Martin, Electronic Structure: Basic Theory and Practical Methods(Cambridge University Press)
  • 先修要求:量子力学、固体物理
  • 考核方式:编程项目(Coding Project)50% + 文献综述与展示(Literature Survey and Presentation)50%
  • 编程项目:给 toy model 写 Hartree-Fock / DFT 代码,要求提交源代码 + 实现报告 + 结果;即使结果不对,提交部分结果并附上原因猜测也可得分;允许使用 AI 工具
  • 文献综述展示:每组两三人,15 分钟展示,主题给出建议但不强制

今天的材料

开讲用几个材料引出”为什么需要第一性原理计算”:

  • 负泊松比(Negative Poisson’s Ratio):存在内禀的负泊松比材料(拉长时横截面反而变宽),这类”反常”力学响应需要从电子结构层面理解;
  • 光伏效应与倍频(Photovoltaic Effect & Second Harmonic Generation):太阳能收集的效率、频率转换系数,都取决于材料的电子结构;
  • 负热膨胀(Negative Thermal Expansion):如 ZrW$_2$O$_8$、ScF$_3$,加热反而收缩。

这些例子想说明:材料的各种宏观性质(力学、光学、热学)归根结底由电子结构决定。

理论中的一切(Theory of Everything)

体系的自由度是电子和原子核,物理为非相对论量子力学,基本方程是薛定谔方程 $H\psi = E\psi$(薛定谔因建立波动力学获 1933 年诺贝尔物理学奖)。

对电子和原子核的非相对论哈密顿量为(此处用原子单位 $\hbar = m_e = e = 1$):

原则上,解这个方程就”解出了一切”——但 $N$ 体波函数 $\psi(\mathbf r_1,\ldots,\mathbf r_N)$ 无法直接处理,这是整个学科的出发点。

指数墙

描述 $\psi(\mathbf r_1,\mathbf r_2,\ldots,\mathbf r_N)$ 所需的参数个数为:

其中 $3N$ 是变量数(忽略自旋),$p$ 是描述单电子波函数到所需精度的参数个数。取 $p = 3$、$N = 100$:$3^{300} \approx 10^{150}$!即使是 $N=100$ 个电子也完全无法存下波函数。

W. Kohn(1999):”一般说来,当 $N>10^3$ 时,$N$-电子波函数不是一个合法的科学概念。”存下 $\psi(\mathbf r_1,\ldots,\mathbf r_N)$ 需要的字节数为 $B = q^{3N}$,其中 $q$ 是描述单电子波函数到所需精度所需的字节数。若用网格存波函数,$q = 100$,4 个电子就要 $10^{13}$ 字节 = 10 TB。

目前基于波函数的方法(QMC、DMRG 等)都在尽力绕过这堵指数墙,但对大多数体系仍然难以胜任。

密度泛函理论的基本思想

关键观察:$N$-电子波函数包含太多信息。我们真正关心的是物理可观测量(电子密度、能隙等)。

密度泛函理论(DFT)的核心主张

  • 基态电子密度 $n(\mathbf r)$ 是唯一需要的信息;
  • 所有物理可观测量都是 $n(\mathbf r)$ 的泛函。

这大大降低了问题的自由度:从 $3N$ 个坐标降到 3 个坐标的函数。W. Kohn(”for his development of the density-functional theory”)与 J. A. Pople(”for his development of computational methods in quantum chemistry”)共享 1998 年诺贝尔化学奖。

涌现物理与 DFT 的位置

涌现(emergent)物理:科学具有层级结构,每一层建立在下层之上(P. W. Anderson,1977 年诺贝尔物理学奖)。”理解涌现态需要新的变量与原理。”

  • DFT 是唯一从量子力学出发预测材料性质的例子——它打破了量子力学与凝聚态物理之间的界限;
  • DFT 的局限:强关联体系(如分数量子霍尔效应)中单粒子图像失效,DFT 无法直接处理;
  • DFT 的正向与反向应用:正向——从原子结构到材料性质的映射;反向——材料设计(material design),根据想要的性能反推原子结构。

课程大纲

主线(对应 Martin 教材):

  • 概览与背景主题:Chs. 1-5
  • 密度泛函理论:Chs. 6-9
  • 原子相关的重要预备:Ch. 11
  • 电子结构的确定:Chs. 12-15
  • 从电子结构到物质性质:Chs. 19-21, 23

附加主题:

  • Hartree-Fock 计算
  • 多体微扰论与激发态

穿插约 2 次实践课(编程项目)和 2 次展示课。

思考题

  1. 为什么说”N-电子波函数在 $N>10^3$ 时不是合法的科学概念”?用 $p^{3N}$ 定量说明。
  2. DFT 为什么选择电子密度而非其他量作为基本变量?密度有哪些”好性质”?
  3. 强关联体系(如分数量子霍尔)为什么会让 DFT 失效?