原子是物质的构建模块，但它们并非学校图表中经常描绘的那种坚实、静态的球体。每个原子的核心都有一个原子核，周围环绕着电子云。在大多数化学计算中，科学家将沉重的原子核视为固定的点，而较轻的电子则在周围飞速穿梭，这种简化方法对许多物质都非常有效。然而，当涉及到氢——这种最轻的元素时，这种简化失效了。因为氢原子核如此之轻，它们的表现并不像微小的、静止的弹珠；相反，它们会抖动并穿越能量势垒，这种行为只能用量子物理学来描述。这种被称为“核量子效应”的行为，对于理解质子如何在光合作用等生物过程中移动，以及水分子如何聚集在一起至关重要。然而，在传统超级计算机上模拟这种具有抖动性的原子核量子特性通常是不可能的，因为其数学复杂度的增长速度过快，难以处理。
来自德国航空航天中心和赫姆霍兹乌尔姆研究所的一个研究小组通过将模拟转移到量子计算机上，在解决这一问题方面迈出了重要一步。他们开发了一种新的混合方法，将量子处理器的力量与经典计算机的可靠性结合起来。在他们的方法中，量子计算机充当一个高度专业化的计算器，用于确定原子在任何给定时刻所受的能量和作用力，而经典计算机则负责处理原子的运动。至关重要的是，这种方法并不将某些原子核视为固定的点，而是将其视为可以同时存在于多个位置的量子对象，从而能够更真实地描绘分子的振动和反应过程。
为了测试他们的算法，研究人员模拟了三个特定系统：一个简单的氢分子、一个水分子和一个容德尔离子（Zundel ion，即由两个共享单个质子的水分子组成的簇）。容德尔离子特别重要，因为它作为质子转移的一个极简模型，而质子转移是化学和生物学中的一个基本过程。研究团队发现，如果忽略原子核的量子性质，他们的模拟产生的结果往往是不准确的。然而，一旦纳入这些核量子效应，分子的行为发生了剧变。在容德尔离子的案例中，纳入这些效应消除了通常阻碍质子在两个水分子之间自由移动的能量势垒。质子的运动不再被困在“山谷”中，而是变得实际上是无势垒的，使其能够轻松地来回穿梭。这种转变与实验数据高度吻合，显著减少了对分子振动频率预测的误差。
研究人员还探索了如何使这些模拟适用于当今现有的、仍处于早期且存在噪声阶段的量子计算机。他们测试了不同的编码问题的方法，将一种高精度但复杂的方法与旨在适应当前硬件的更简单、更紧凑的方法进行了对比。他们发现，虽然复杂的方法提供了最精确的结果，但更简单的紧凑方法却出人意料地有效。对于水分子和氢分子，这些精简后的方法产生的结果几乎与复杂方法一样好，且所需的计算资源要少得多。这表明该算法可以在近期的量子设备上运行，即使是那些功率有限的设备。
然而，通往完美准确性的道路并非没有挑战。团队发现，模拟的成功很大程度上取决于他们如何选择表示原子核的量子态。对于水分子，简单的表示方式就足够了；但对于更复杂的容德尔离子，他们需要包含更详细的核轨道描述，以捕捉正确的相互作用。当他们使用更详细的描述时，模拟成功重现了离子振动光谱中特定的双峰模式，这与现实世界实验中观察到的现象相匹配。如果没有这些额外的细节，模拟就会错过这些关键特征。研究结论指出，虽然目前的量子计算机尚不足以模拟大型化学系统的每一个细节，但这种混合方法提供了一条可行的路径。通过仔细选择系统中哪些部分在量子计算机上模拟，哪些部分由经典计算机处理，科学家们可以开始以一种此前无法达到的精度来模拟原子核的量子行为。
技术摘要：量子计算机上的核效应分子动力学
问题陈述
核量子效应（NQEs），如质子隧穿和零点能，对于准确描述质子转移和氢键等过程至关重要。然而，在经典硬件上将这些效应纳入量子化学计算通常在计算上是难以承受的，因为这会导致全构型相互作用（FCI）空间呈指数级增长。虽然量子计算机由于其对系统规模的量子比特需求呈线性缩放，提供了一个潜在的解决方案，但在涉及包含 NQEs 的动力学模拟中展示实际的量子优势仍然是一个挑战。现有的算法往往缺乏对核量子效应的包含、动力学模拟，或者依赖于预先确定的轨迹。
方法论
作者提出了一种新型的混合量子-经典算法，用于通过核-电子轨道（NEO）框架集成 NQEs 的从头算分子动力学（AIMD）。
- NEO 框架： 被选中的轻核（具体为本研究中的质子）与电子一起使用类似于电子原子轨道的核轨道进行量子力学处理。这避免了针对这些特定核的标准 Born-Oppenheimer 分离。
- 混合动力学算法： 该算法遵循基于经典 Born-Oppenheimer 动力学的第二量子化方法。
- 量子组件： 在量子计算机上采用变分量子特征值求解器（VQE）来确定瞬时基态能量（E0）以及作用于核的力。力的计算通过利用海尔曼-费曼定理（Hellmann-Feynman theorem）评估力算符 ∇IH(R) 在基态 ∣ψ0(R)⟩ 上的期望值来实现。
- 经典组件： 分子几何结构使用基于量子处理器计算出的力，通过速度韦尔算法（velocity verlet algorithm）进行更新。
- 活性空间技术： 为了应对噪声中等规模量子（NISQ）设备有限的量子比特数，作者实现了哈特里-福克（Hartree-Fock）嵌入。量子计算机选择一个活性空间，而未参与的轨道则在哈特里-福克水平上进行处理。选择标准包括单轨道熵、沿反应路径的变化性以及对称性保持。
- Ansatz 选择： 研究使用了具有化学启发意义的幺正耦合簇（UCCSD/NEOUCCSD）ansatz 以及硬件高效型 ansatz（实振幅 ansatz 和平铺幺正乘积 ansatz - tUPS），以评估在近程设备上的可行性。
- 基组： 模拟采用了 6-31G 电子基组以及特定的核基组，后者衍生自 DZSPDN（仅针对 s 型的 DZSNB；针对 s 型和 p 型的 DZSPNB）。
核心贡献
- 算法实现： 本工作展示了首次在量子计算机上使用 NEO 方法显式包含核量子效应的混合量子-经典 AIMD 算法实现。
- 活性空间策略： 作者开发了一种专门用于 NEO 动力学中选择活性空间的协议，该协议考虑了轨道变化性和对称性，确保了动力学过程中的物理一致性。
- 硬件效率： 研究表明，紧凑的硬件高效型 ansatz 可以达到与更精确的 UCCSD ansatz 相当的结果，这表明其在近程硬件上执行的可行性。
- 基组敏感性： 论文强调了模拟精度对核基组选择（特别是 p 型轨道对于复杂模式耦合的必要性）和活性空间的依赖性。
结果
该算法通过对氢分子（H2）、水（H2O）和 Zundel 离子（H5O2+）三个系统的模拟进行了验证。
- 氢（H2）： 通过 NEO 引入 NQEs 修正了纯电子方法中高估振动频率的问题，尽管这导致了轻微的低估。硬件高效型 ansatz（tUPS 和 Real Amplitudes）在 NEO 情况下的误差低于 1%，这归功于 ansatz 高估与 NEO 低估之间的偶然误差抵消。
- 水（H2O）： 对于所有简正模（弯曲、对称伸缩、反对称伸缩），NEO 模拟产生的相对误差均低于纯电子方法。涉及质子-质子近距离相互作用的弯曲模仍然具有挑战性，但显示出了改进。
- Zundel 离子（H5O2+）：
- 势能面： 在纯电子情况下，质子转移势表现出带有势垒的双阱结构。在引入 NQEs 后，质子穿梭运动变得实际上是无势垒的，使势能面转变为单阱。这表明共享质子的零点运动产生了显著影响，使振动基态高于纯电子势垒。
- 光谱： 纯电子模拟未能捕捉到 1000 cm−1 附近的特征光谱双峰以及质子转移与水分子摆动模式之间的耦合。
- NEO 成功案例： 使用 p 型核轨道的 NEO 模拟成功重现了光谱双峰以及质子转移模式与水分子摆动/伸缩模式之间的耦合。为了解决较大活性空间的优化问题，必须使用缩减后的 NEOUCCSD ansatz（从 63 个激发中选择 17 个）。
- 准确度： 与活性空间内的精确 NEOCASCI 解相比，NEO 结果显示每个峰的最大误差小于 40 cm−1，并且与实验预解离光谱具有定性一致性。
意义与主张
该论文声称其混合算法成功证明了在量子计算机上模拟包含核量子效应的分子动力学的可行性。作者强调：
- NQEs 对于正确建模 Zundel 离子的质子转移至关重要，特别是通过消除人工能量势垒并实现正确的模式耦合。
- 虽然 UCCSD 提供了最高的准确度，但硬件高效型 ansatz 为近程设备提供了一条可行的路径，前提是必须仔细选择活性空间和核基组。
- 核基组（特别是 Zundel 离子中 p-轨道的包含）的选择比所使用的特定 ansatz 对准确性的影响更为关键。
- 这项工作为未来更大规模系统的研究奠定了基础，并指出通过从正则轨道基组向自然轨道基组过渡，可以进一步提高资源利用率或准确度。
作者对当前的技术状态保持谦逊，承认其结果依赖于无噪声的状态矢量模拟器，且尚未在物理硬件上展示出显著的量子优势。这项工作作为将量子化学方法扩展到动态包含核自由度的概念验证。
每周获取最佳 quantum physics 论文。
受到斯坦福、剑桥和法国科学院研究人员的信赖。
请查收邮箱确认订阅。
出了点问题，再试一次？
无垃圾邮件，随时退订。