2. 化学问题的量子线路资源抽象
在量子计算化学中,一个分子问题最终能否在真实硬件上运行,并不取决于化学式写得多漂亮,而取决于它被翻译成量子线路之后消耗了多少资源。资源估算工程的第一步,就是建立一条从“分子哈密顿量”到“量子线路”的抽象链路,并在这条链路的每个环节上明确可计量的资源指标:量子比特数、门数、线路深度、测量次数。只有把这四个量算清楚,才能判断某个分子体系是否落在当前或近期硬件的可行区间内。
一、四个核心资源指标的含义
量子比特数是线路所需的逻辑量子比特总量,它决定了问题能否被装进硬件。对化学问题而言,比特数主要来自两部分:一是表示电子占据状态的寄存器,二是为编码哈密顿量、执行算法(如量子相位估计或变分方法)所需的辅助比特。以 Jordan-Wigner 变换为例,N 个自旋轨道需要 N 个量子比特;若采用 Bravyi-Kitaev 变换,比特数相同但局域性更好,会显著影响后续门数。
门数是线路中所有量子门的总和,通常分为单比特门和双比特门。双比特门(如 CNOT、CZ)在超导硬件上误差高、耗时长,因此工程上更关注双比特门计数。门数直接决定线路的执行时间与累积错误。
线路深度是线路中可并行执行的门被压缩后,从输入到输出的最长依赖路径长度。它反映的是“时间维度”上的资源:深度越大,退相干影响越严重。门数多但并行度高的线路,深度可能反而小于门数少但串行执行的线路。
测量次数(shots)是获取期望值所需的重复制备与测量次数。由于量子测量具有随机性,每个可观测量的期望值都需要多次采样估计。测量次数与哈密顿量分解出的项数、每项的方差以及目标精度直接相关。
二、从哈密顿量到线路的抽象链路
分子哈密顿量经过二次量子化后,通常写成费米子算符的线性组合。要把它变成量子线路,需要经历几个关键步骤:先选择费米子到量子比特的映射(Jordan-Wigner、Bravyi-Kitaev 或 parity 编码),把费米子算符替换为泡利算符的张量积;再对哈密顿量做分组,把可对易的泡利项合并,减少需要单独测量的组数;然后为每一项设计对应的量子线路,通常由基变换门和纠缠门构成;最后把算法外层(如 UCCSD 拟设、相位估计的受控演化)与哈密顿量项线路组合起来,形成完整线路。
这条链路上,每一步都会改变资源指标。映射方式影响门数和局域性;分组策略影响测量次数;拟设的截断程度影响线路深度;而辅助比特的引入方式则影响总比特数。资源估算不是对最终线路做一次静态统计,而是要在设计选择之间做权衡。
三、一个可操作的估算示例
以氢分子 H₂ 在 STO-3G 基组下的最小基组为例。H₂ 有 2 个电子、4 个自旋轨道(2 个成键、2 个反键)。采用 Jordan-Wigner 映射时,需要 4 个量子比特。经过对称性约简(冻结核心轨道、利用粒子数守恒和自旋对称性),实际可压缩到 2 个量子比特。哈密顿量在 JW 映射下展开为约 15 个泡利项,分组后通常可合并为 3~5 组可对易项。
若采用 UCCSD 拟设,H₂ 的激发算符只有单激发和双激发各一项,对应线路中约 2~4 个 CNOT 门,深度约 10~20。若采用更简单的硬件高效拟设,深度可压到 5~8,但精度需要额外优化。测量次数方面,若目标化学精度为 1.6 mHa(约 0.0016 Hartree),每组期望值的采样次数通常在 10⁴~10⁵ 量级,总测量次数为组数乘以单组采样数。
这个例子说明:同一个分子,在不同映射、不同拟设、不同分组策略下,资源指标可以相差数倍。资源估算的价值就在于把这些选择量化,而不是凭感觉判断“能不能跑”。
四、估算时容易忽略的细节
第一,辅助比特的隐性开销。相位估计需要额外的 ancilla 寄存器来存储相位,其数量与目标精度成对数关系;变分方法虽然不需要 ancilla,但参数优化本身需要大量线路执行,等效测量次数会显著增加。
第二,双比特门的实际代价。在超导硬件上,CNOT 门保真度通常比单比特门低一个数量级,且持续时间长。因此门数统计中应单独列出双比特门数量,而不是只看总数。
第三,深度与并行度的关系。哈密顿量分组后,不同组的线路可以并行执行,但同一组内的项往往需要串行。估算深度时,应按“组内串行、组间并行”的模型计算,而不是简单地把所有门数相加。
第四,测量次数的方差依赖。每个泡利项的测量方差与其系数平方成正比。系数大的项需要更多采样,系数小的项可以少采样。若采用加权采样策略,总测量次数可以低于均匀分配。
五、把资源指标写成可复用的估算表
为了让估算结果可比较、可复现,建议为每个分子体系建立一张资源估算表,至少包含以下字段:分子名称与基组、自旋轨道数、映射方式、约简后量子比特数、哈密顿量泡利项数、分组数、拟设类型、双比特门数、线路深度、单组采样次数、总测量次数、目标精度。表中每个数值都应标注来源(理论计算、文献参考或合理假设),避免把估算值误当作实测值。
例如,对 H₂/STO-3G 可以填写:自旋轨道 4,JW 映射,约简后 2 比特,泡利项约 15,分组 4,UCCSD 拟设,双比特门约 4,深度约 15,单组采样 20000,总测量 80000,目标精度 1.6 mHa。这些数值属于教学示例中的合理估算,实际运行时会因编译优化、硬件拓扑和采样策略而有所浮动。
当这张表建立起来之后,不同分子之间就可以横向比较:LiH 需要多少比特、多少双比特门;BeH₂ 的深度是否超出当前硬件的相干时间;某个候选分子是否值得投入更多采样。资源抽象的意义,正是把化学问题转化为工程上可度量、可比较、可决策的指标集合。