1. 课程导入与全景概览

量子化学的核心任务,是在给定原子核构型下求解电子薛定谔方程,从而获得分子或材料的能量、电子密度、轨道能级与各类可观测性质。而这一方程的全部物理信息,都浓缩在一个算符之中——哈密顿量。可以说,量子化学计算的一切后续步骤,本质上都是在“如何表示哈密顿量”和“如何求解哈密顿量”这两个问题上做文章。本课程《量子化学哈密顿量构建与映射》正是围绕这两个问题展开:从二次量子化形式出发,构建分子电子哈密顿量,再通过不同映射方案将其转化为量子计算机可执行的量子比特算符,最终服务于变分量子本征求解(VQE)、量子相位估计(QPE)等算法。

一、哈密顿量是什么:从物理定义到算符形式

在非相对论、玻恩-奥本海默近似下,一个含 N 个电子、M 个原子核的分子体系,其电子哈密顿量可写为三项之和:电子动能、电子与核的吸引势能、电子间库仑排斥能。若采用原子单位(ħ = 1,me = 1,e = 1),其第一量子化形式为:

H = -Σ_i (1/2)∇_i² - Σ_i Σ_A (Z_A / r_iA) + Σ_{i<j} (1 / r_ij)

其中 i、j 遍历电子,A 遍历原子核,ZA 为核电荷数。这一形式直观但难以直接用于量子计算,因为电子是全同费米子,波函数必须满足反对称性。于是我们转向二次量子化:引入产生算符 ap† 和湮灭算符 aq,哈密顿量变为

H = Σ_{pq} h_{pq} a_p† a_q + (1/2) Σ_{pqrs} h_{pqrs} a_p† a_q† a_r a_s

这里 hpq 是单电子积分(动能加核吸引),hpqrs 是双电子积分(电子排斥)。这一步的意义在于:哈密顿量被离散化为一组积分系数与算符乘积的线性组合,而积分系数可由经典量子化学程序(如 PySCF、Psi4)在给定基组下算出。至此,问题从“解微分方程”转化为“处理算符代数”。

二、哈密顿量在电子结构计算中的核心地位

为什么说哈密顿量是电子结构计算的“心脏”?因为一旦哈密顿量确定,体系的所有信息都隐含其中:基态能量是它的最小本征值,激发态对应更高本征值,电子密度可由基态波函数求得,偶极矩、极化率等响应性质则通过微扰展开与哈密顿量的导数关联。经典方法如 Hartree-Fock、CI、CCSD、DFT,本质上都是对同一哈密顿量采用不同近似求解策略;量子算法如 VQE、QPE,则是换用量子硬件去逼近它的本征值。因此,哈密顿量构建的精度与映射的效率,直接决定后续计算的可行性与准确度。

学习本课程时,建议始终把“哈密顿量”当作一个可操作的对象:它由积分系数和算符结构共同定义,而不是一个抽象符号。后续每一章都会围绕它的某个具体环节展开。

三、构建与映射的整体流程

从分子结构到量子线路,完整流程可概括为五个阶段:

  1. 几何与基组输入:给定原子坐标、电荷、自旋多重度,选择基组(如 STO-3G、6-31G)。
  2. 积分计算:用经典量子化学程序计算 hpq 与 hpqrs,并做对称性约化。
  3. 二次量子化哈密顿量:将积分组装成费米子算符形式。
  4. 费米子-量子比特映射:采用 Jordan-Wigner、Bravyi-Kitaev 或 Parity 等映射,把 a、a† 转成 Pauli 算符(X、Y、Z 及其张量积)。
  5. 量子算法求解:在量子线路中制备试探态,测量期望值,优化参数,得到基态能量。

下图展示了这一流程中各阶段的输入、输出与关键对象,帮助建立整体认知框架。

量子化学哈密顿量:构建与映射全景流程 ① 几何与基组 原子坐标 / 电荷 STO-3G / 6-31G ② 积分计算 h_pq / h_pqrs PySCF / Psi4 ③ 二次量子化 a† a 算符组合 费米子哈密顿量 ④ 量子比特映射 JW / BK / Parity Pauli 算符串 ⑤ 量子算法求解 VQE / QPE → 基态能量 E₀ 测量期望值 ⟨H⟩ 并优化参数 各阶段关键对象与输出 分子结构 坐标 / 基函数 电子数 / 自旋 积分张量 h_pq (N²) h_pqrs (N⁴) 费米子算符 a† a 项 反对易关系 Pauli 算符 X / Y / Z 张量积 量子比特数 n 能量本征值 E₀ / E₁ ... 化学精度 1.6 mHa

四、课程学习地图与工具链简介

本课程共 30 章,整体分为四个模块。第一模块(第 1–8 章)建立基础:二次量子化、积分计算、费米子算符代数、Jordan-Wigner 与 Bravyi-Kitaev 映射原理。第二模块(第 9–16 章)深入映射技术:Parity 映射、对称性约化、Qubit 数压缩、泡利算符分组与测量策略。第三模块(第 17–24 章)面向算法:VQE 线路设计、UCCSD 拟设、QPE 与量子线路实现、噪声影响评估。第四模块(第 25–30 章)综合实战:小分子(H₂、LiH、H₂O)全流程案例、结果对比与误差分析、工具链集成与性能调优。

工具链方面,经典积分计算推荐 PySCF(轻量、Python 原生)或 Psi4(功能全面);费米子-量子比特映射与线路构建可使用 OpenFermion、Qiskit Nature 或 PennyLane;量子线路模拟与硬件执行可借助 Qiskit Aer、Cirq 或 Amazon Braket。建议读者在本地安装 PySCF + OpenFermion + Qiskit 组合,后续章节的示例将基于这一环境展开。

本课程不要求读者预先掌握量子计算硬件细节,但需要具备基本的线性代数、量子力学符号(狄拉克符号)与 Python 编程能力。每章会给出可运行的代码片段与参数示例,帮助把抽象算符落到具体数值上。

最后强调一个贯穿全课程的原则:哈密顿量的构建与映射不是两个孤立步骤,而是一个精度与资源相互权衡的连续过程。基组越大,积分越准,但量子比特数越多;映射越紧凑,线路越浅,但可能引入额外对称性约束。理解这一权衡,是后续所有章节的主线。