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