Digital Discovery | 量子计算驱动的高精度非绝热分子动力学模拟
近日，北京大学袁骁课题组与合作团队提出了一套面向量子计算的实用化非绝热分子动力学框架，首次实现了Sub-microhartree精度级别的势能面计算与高效、稳定的非绝热动力学模拟。该成果以”Efficient Quantum Simulation of Non-Adiabatic Molecular Dynamics with Precise Electronic Structure”为题，发表于综合化学领域的著名期刊《Digital Discovery》，并入选2026年2月第2期封面论文，第一作者为北京大学博士生李天翊。
研究背景
光合作用、光催化、电荷转移——这些重要的化学过程背后，隐藏着电子态之间的快速跃迁。这种”非绝热效应”是理解和调控光化学反应的关键钥匙，然而精确模拟这些过程却是一项艰巨挑战：求解多个电子态的势能面，其计算复杂度随系统规模呈指数增长，经典计算机在这一领域早已触及能力边界。量子计算凭借其独特的量子叠加与纠缠特性，被寄予突破这一困境的厚望。但现实中，如何设计一套足够通用、高效且精确的量子算法，使其能够适配多样化的化学体系，并与成熟的非绝热分子动力学理论无缝衔接，仍是一个核心议题。
方法概述
北京大学袁骁课题组联合多家研究机构，构建了一套”混合子空间量子电子结构求解器+Landau-Zener势能面跃迁动力学”的融合计算框架。在电子结构计算层面，团队采用变分量子本征求解器（VQE）及其子空间变体（SSVQE）在选定的活性空间内高效搜索基态或多个正交参考态，随后通过量子子空间展开（QSE）方法将哈密顿量投影到由参考态及其激发构成的子空间中，经典对角化得到激发态能量。针对不同化学体系的特性，研究者还开发了灵活的算符扩展策略（QSE*），显著提升了复杂体系的计算精度。
在动力学层面，团队将量子电子结构求解器与Landau-Zener势能面跃迁（LZSH）框架深度整合。LZSH方法的精妙之处在于，它仅依赖能隙及其导数信息即可处理态间跃迁，无需计算昂贵的非绝热耦合向量，天然适合与量子算法协同工作。针对势能面在解离区域可能出现的数值不连续性，研究者创新性地引入”曲率驱动的跳跃校正”机制：通过监测能隙曲率的相对变化，自动识别并过滤由数值噪声引起的伪交叉，同时保留物理意义明确的真实跃迁事件，为量子计算非绝热分子动力学的稳定性提供了保障。
研究成果
该框架在多个具有代表性的分子体系中展现了性能。在H₃⁺离子的解离路径上，扩展的QSE方法对S₁和S₂激发态均能控制在”Sub-microhartree精度”($10^{-6}$量级)。即使在锥形交叉区域采用0.0005 Å的极细位移步长，QSE*仍保持稳定的高精度输出。
对于更具挑战性的乙烯分子（C₂H₄）”金字塔化”几何扭转轨迹，SSVQE-QSE混合方法精确再现了S₀与S₁态之间的锥形交叉行为。SSVQE展现出稳健的态追踪能力，能够在交叉前后持续锁定目标电子态，充分证明了该框架处理有机分子光化学过程的潜力。
在非绝热动力学的集成测试中，H₃⁺的50条轨迹量子模拟器仿真显示，应用曲率校正的VQE-QSE*方法产生了平滑的布居演化曲线，与精确参考结果趋势高度一致：S₂态稳定退布居，S₁态先升后降，S₀态逐步累积——这正是非绝热弛豫过程应有的物理图像。
在计算效率方面，研究团队实现了电子结构与轨迹两个层面的并行优化。QSE矩阵元的估计可完全并行化，利用哈密顿量的厄米特性将矩阵估计开销减半，自旋适配版本进一步实现4倍算符缩减；表面跳跃集成中各轨迹相互独立，可分布式并行执行。量子模拟器上的基准测试显示，并行QSE相比串行版本实现约12倍加速，Hellmann-Feynman方法计算核力比有限差分法快约41倍，整体框架可为典型光化学应用提供1-2个数量级的计算时间缩减。
研究意义
这项研究不仅是算法层面的技术创新，更为量子计算在计算化学领域的实际应用探索了道路。研究团队发展了从量子电子结构计算到分子动力学模拟的完整技术路线，并行框架可无缝迁移至真实量子硬件。更重要的是，该方法已在离子体系（H₃⁺）、有机分子（C₂H₄）、含杂原子体系（CH₂O）上验证有效性，展现出广泛的通用性与强大的扩展潜力。未来，这一框架有望在光催化机理研究、光敏材料设计、药物光化学行为预测等领域发挥作用，推动量子计算从理论验证走向实际应用。
论文信息
Efficient Quantum Simulation of Non-Adiabatic Molecular Dynamics with Precise Electronic Structure, Digital Discovery, 2026, Issue 2 (Cover Article)