基于准绝热量子几何优化工作流的高通量虚拟筛选

《Journal of Chemical Theory and Computation》：Toward High-Throughput Virtual Screening via Quasi-Adiabatic Quantum Geometry Optimization Workflow

　　

这项工作致力于利用量子计算方法学应对post-Hartree–Fock（后HF）理论水平下的分子几何优化挑战。研究人员提出了一种基于变分量子本征求解器（VQE）算法的具有潜在准绝热性的新方法，用于高效探索势能面（PES）。该方法将量子线路模拟器与经典优化程序相结合，并利用Hellmann–Feynman定理（含Pulay修正）以及有限差分格式来计算力。研究人员的实现采用并行化框架内的GPU加速量子线路运算，通过对波函数进行热启动（warm-starting）并在后续每一步几何优化中对其重新优化，从而获得计算效率。针对小分子体系的基准测试结果证明了这种混合量子—经典方法的可行性，并与已建立的量子化学方法（包括CASCI和DFT优化）进行了详细对比。研究人员讨论了该方法在药物发现应用的量子化学背景下的当前局限性、计算优势以及潜在的未来发展方向。

**论文解读：基于准绝热量子几何优化工作流的分子结构优化研究**

**一、研究背景与意义**

在计算机辅助药物发现（CADD）领域，高通量虚拟筛选（HTVS）是过滤海量化合物的核心技术手段，其目的在于从巨大化学空间中识别能够与特定生物大分子靶标结合的先导化合物。虚拟筛选流程的第一步通常是基于轻量级SMILES字符串生成候选分子的三维结构，这一步骤往往依赖经典分子力场对势能面（PES）进行局部最小化搜索。然而，基于经典力场恢复的分子结构忽略了一些关键的化学细节，特别是量子核效应（如零点能），这在富含氢原子的分子中尤为重要。更关键的是，经典力场所描述的势能面与真实量子势能面存在偏差，导致后续用于结合亲和力预测的打分函数（SF）产生大量假阳性和假阴性结果，从而浪费宝贵的实验资源。

从量子化学的角度看，分子几何优化在post-HF理论水平下是一项计算挑战。传统的post-HF方法（如耦合簇CCSD(T)）尽管精度高，但计算成本随体系规模呈指数级增长；而密度泛函理论（DFT）虽然计算效率较高，但在处理强关联电子体系（如过渡金属配合物、金属蛋白等）时存在固有局限。当前最先进的方法通常依赖DFT进行常规应用，或在需要更高精度时采用昂贵的CC方法。量子计算的兴起为解决这些挑战提供了新途径，变分量子本征求解器（VQE）已成为近期量子化学应用中最常用的量子算法之一。基于此，研究人员提出了一种基于VQE的准绝热量子几何优化（QGO）方法，旨在为HTVS提供更精确的输入结构，最终提高虚拟筛选的可靠性。

**二、研究内容**

本研究开发了一套集成的混合量子—经典框架（QGO），用于post-HF水平的分子几何优化。该框架将量子线路模拟（GPU加速）用于电子结构计算，将经典优化算法用于化学空间几何参数更新。具体包括三个核心组件：基于有限差分近似的核力计算、绝热波函数追踪策略以加速收敛，以及非绝热探索策略以防止变分优化陷入局部极小值。研究人员还特别处理了分子轨道相位一致性问题，确保在有限位移几何下构建的量子哈密顿量之间的一致性。以Dompe Keys（DKs）分子碎片为测试体系，系统评估了该方法在亚甲亚胺、乙醛和乙酰胺三种代表性有机分子上的性能，并与CASCI、CASSCF、CCSD和DFT/B3LYP等标准量子化学方法进行了对比。

**三、主要技术方法**

研究人员从CCCBDB数据库获取初始几何构型（均采用Restricted-HF/STO-3g优化），在cc-pVDZ基组下构建CAS(6,6)活性空间，经Jordan–Wigner变换映射到12个量子比特。核心算法为VQE-D框架，采用CNOT梯式硬件高效拟设（HEA），深度为4。力计算采用中心有限差分法（位移步长为δ），每个梯度需6N_atoms?12次（内坐标）VQE求值；通过GPU加速（CudaQ或Qiskit statevector模拟器）和joblib多进程并行化，实现了位移几何间波函数的热启动与轨道相位一致性校正。优化器采用L-BFGS-B，50次VQE初始化搜索与每3次宏迭代的非绝热分支探索。

**四、研究结果**

**4.1 亚甲亚胺（Methylene Imine）**

QGO优化结果复现CASCI参考平衡结构，键长偏差在0.003–0.017 ?范围内（远低于0.02 ?化学精度阈值），键角差异均在1°以内。相比HF/STO-3g初始几何，相关方法（QGO与CASCI）给出略长的键长（最多约0.03 ?）和2–3°的小角度调整，反映了电子相关和基组效应的纳入。这一结果表明绝热波函数追踪结合有限差分力计算能够有效捕获分子的相关平衡构型，且精度与CASCI相当。

**4.2 乙醛（Acetaldehyde）**

QGO复现CASCI几何，键长差异在0.002–0.009 ?，键角偏差在1°以内，两者几乎完全一致。与初始HF/STO-3g结构相比，两种方法均预测稍短的C–O键与更长的C–C和C–H键，这和电子关联效应及cc-pVDZ基组灵活性提升相一致。CASSCF结果与CASCI接近但键长略增，而CCSD和B3LYP呈现预期的系统键长延长趋势，这源自其相关处理方式的不同。该结果表明，QGO能同时捕获局域和全局的分子结构特征。

**4.3 乙酰胺（Acetamide）**

QGO与CASCI的键长偏差为0.002–0.015 ?，键角差异约1–2°，定量精度优良。两种方法均预测C=O和N–H键略为缩短、C–C键适度延长。值得注意的是，QGO复现了CASCI中酰胺键长（如C2–N1键1.392 ?）的变化趋势，而这种酰胺共振导致的浅势能面和微小梯度正是量子优化最需要稳健性的场景。CASSCF结果与CASCI一致但呈系统性小幅偏移，CCSD与B3LYP同样显示键长伸长规律。

**4.4 计算效率**

GPU加速显著缩短了VQE模拟的墙钟时间，并行执行相比串行（时间过长而无法显示）实现了数量级加速。在GPU V100和H100上的VQE循环运行时间均远优于AMD EPYC 9554 CPU，部分测试在CPU上超出20000秒时间限制无法完成，而在GPU上可正常运行，表明GPU加速是实现实用化量子几何优化的重要支撑。

**4.5 已识别的挑战**

当前实现存在多方面局限：有限差分力计算带来大量额外VQE优化需求，尽管并行化缓解了墙钟时间，但渐近复杂度未变；在存在多个近简并电子态或接近键解离极限的区域出现收敛困难；硬件高效拟设（HEA）虽提供参数灵活性但化学表达能力有限；小活性空间可能不足以描述完整电子关联。活性空间选择和基组效应也限制了达到化学精度的能力。此外，波函数之间SCF解的相干性对力的计算精度有重要影响，相位不一致会引起力的较大失衡。

**五、讨论与结论**

这项工作证明了基于VQE的post-HF分子几何优化的可行性。准绝热策略与GPU加速实现的结合，能够通过跨几何优化步骤回收中心波函数，显著加速计算，且该方法适用于与药物相关有机分子碎片，突破了先前文献（如H₂、H₃⁺、H₄等）仅限简单元素体系的局限。该方法原生适配高性能计算（HPC）架构，具有两个层次的并行性：多起点搜索与梯度分量的独立计算。QGO研究促进量子计算方法在药物发现中应用的潜力，以经典可行的计算开销获得post-HF精度，为未来在强关联体系、激发态优化和反应路径表征等方向的研究奠定基础。未来工作将重点发展解析梯度法替代有限差分力计算、探索ADAPT-VQE等自适应拟设构造策略、扩展到更大分子体系，以及纳入量子纠错技术。尽管当前实现仍限于经典模拟，但该框架为未来量子硬件可用时的部署提供了基础。