1. 分子能量面与反应路径量子模拟导论

分子能量面（Potential Energy Surface, PES）是理论化学中最核心的概念之一：它把原子核坐标映射为体系电子能量，从而把“分子如何运动、如何反应”转化为一个可计算、可分析的几何对象。对 N 个原子的非线性分子，若忽略整体平动与转动，其内部自由度有 3N−6 个，因此势能面本质上是高维空间中的标量场。反应路径则是这个高维面上连接反应物与产物的最低能量通道，通常以过渡态为鞍点、以反应坐标 s 为参数来描述。经典计算化学通过量子化学方法（如 Hartree–Fock、DFT、CCSD(T)）在特定核构型下求解电子能量，再通过几何优化、频率分析和内禀反应坐标（IRC）追踪来刻画路径。然而，当体系包含强电子相关、多参考特征或非绝热跃迁时，经典方法在精度与成本之间往往难以兼顾。

量子模拟的介入点正在于此。量子计算机天然适合表示和演化量子态，而电子结构问题本身就是量子多体问题。将分子哈密顿量映射为量子比特形式（如 Jordan–Wigner 或 Bravyi–Kitaev 变换），再通过变分量子本征求解器（VQE）、量子相位估计（QPE）或量子虚时演化（QITE）等算法求取基态与激发态能量，可以在原则上以多项式资源处理经典方法指数级困难的问题。对于势能面研究，这意味着可以在多个核构型上重复求解电子能量，从而“逐点”构建能量面；对于反应路径，则可以在候选路径上采样并比较能量，定位过渡态与中间体。量子模拟并非替代经典方法，而是与经典几何优化、采样和统计方法形成混合工作流。

本课程的整体学习路线围绕“能量面—反应路径—量子算法—工具链”四条主线展开。首先建立势能面与反应路径的数学描述，包括驻点、Hessian 矩阵、鞍点阶数、内禀反应坐标和最小能量路径（MEP）的定义。随后进入量子模拟基础：二次量子化、费米子—量子比特映射、变分原理与 ansatz 设计。在此基础上，重点讲解如何用量子算法计算单点能量、如何扫描核坐标构建势能曲线、如何通过有限差分或量子力学的解析梯度获取力，以及如何结合经典优化器寻找极小值与鞍点。最后，通过工具链实践把上述环节串起来：使用 OpenFermion/PySCF 生成分子积分，用 Qiskit Nature 或 PennyLane 构建量子线路，在模拟器或真实硬件上运行，并用经典后处理完成路径分析。

为了把上述流程落到可执行层面，下面给出一个最小但完整的示例任务：以 H₂ 分子在 STO-3G 基组下的势能曲线为例，说明如何从几何输入走到能量输出。首先用 PySCF 在键长 R = 0.5–2.5 Å 范围内取若干点，每个点计算一次 Hartree–Fock 与 FCI 能量，得到参考势能曲线；随后用 OpenFermion 生成二次量子化哈密顿量，经 Jordan–Wigner 变换得到量子比特形式；再用 PennyLane 构建含 2 个量子比特的 VQE 线路，ansatz 采用硬件高效形式，优化器选 COBYLA，最大迭代 200 次。运行后记录每个键长下的基态能量，并与 FCI 对比。判断标准是：能量误差小于 1.6 mHartree（化学精度）时认为该点可用；若误差偏大，则增加 ansatz 层数或改用 UCCSD 型线路。这个例子虽然简单，但它完整覆盖了“几何输入—积分—映射—线路—测量—误差判断”的闭环，后续章节会在此基础上扩展到多原子分子、激发态与反应路径。

关于工具链的选择，OpenFermion 适合做费米子算符与映射的底层操作，Qiskit Nature 提供了从分子到线路的较完整封装，PennyLane 则在可微量子编程与经典优化器联动方面更灵活。实际项目中可以混合使用：用 PySCF 做经典预处理，用 OpenFermion 做映射，用 PennyLane 做变分优化，最后用经典脚本扫描核坐标并拼接势能曲线。需要特别注意的是，量子模拟得到的能量是电子能量，若要与实验活化能比较，还需加上零点能校正与溶剂化效应等经典修正。