2025-11-13T18:37:11.036211

Uniformity Bias in Ground-State Sampling Induced by Replica Alignment in Quantum Monte Carlo for Quantum Annealing

Maruyama, Ohzeki, Tanaka

Quantum annealing (QA) with a transverse field often fails to sample degenerate ground states fairly, limiting applicability to problems requiring diverse optimal solutions. Although Quantum Monte Carlo (QMC) is widely used to simulate QA, its ability to reproduce such unfair ground-state sampling remains unclear because stochastic and coherent quantum dynamics differ fundamentally. We quantitatively evaluate how accurately QMC reproduces the sampling bias in QA by comparing the final ground-state distributions from the QMC master equation and the SchrÃ¶dinger equation. We find QMC tends to produce uniform ground-state probabilities, unlike QA's biased distribution, and that this uniformity bias strengthens as annealing proceeds. Our analysis reveals that this bias originates from replica alignment -- the dominance of configurations in which all Trotter replicas coincide -- caused by the energetic suppression and entropic reduction of kink configurations (replica mismatches). These findings clarify a fundamental limitation of discrete-time QMC in faithfully simulating QA dynamics, highlighting the importance of replica correlations and transition rules in achieving realistic ground-state sampling.

academic

Uniformity Bias in Ground-State Sampling Induced by Replica Alignment in Quantum Monte Carlo for Quantum Annealing

基本信息

论文ID: 2510.10566
标题: Uniformity Bias in Ground-State Sampling Induced by Replica Alignment in Quantum Monte Carlo for Quantum Annealing
作者: Naoki Maruyama, Masayuki Ohzeki, Kazuyuki Tanaka
分类: quant-ph cond-mat.stat-mech
发表期刊: Journal of the Physical Society of Japan
论文链接: https://arxiv.org/abs/2510.10566

摘要

量子退火(QA)在横向磁场作用下经常无法公平地采样简并基态，这限制了其在需要多样化最优解问题中的适用性。虽然量子蒙特卡罗(QMC)被广泛用于模拟QA，但由于随机和相干量子动力学的根本差异，其重现这种不公平基态采样的能力仍不明确。本研究通过比较QMC主方程和薛定谔方程的最终基态分布，定量评估了QMC重现QA采样偏差的准确性。研究发现QMC倾向于产生均匀的基态概率，与QA的偏置分布不同，且这种均匀性偏差随着退火过程的进行而加强。分析揭示这种偏差源于副本对齐——所有Trotter副本重合的配置占主导地位——这是由扭结配置(副本不匹配)的能量抑制和熵减少引起的。

研究背景与动机

问题定义

核心问题: 量子退火在处理具有简并基态的问题时存在不公平采样现象，即不同基态被采样的概率不相等
模拟挑战: QMC作为模拟QA的主要工具，其随机马尔可夫演化与QA的相干量子演化存在根本差异
实际需求: 许多应用场景(如SAT滤波器、机器学习模型)需要均匀获取多个最优解

研究重要性

理论意义: 理解QMC模拟QA的准确性和局限性对于量子计算理论发展至关重要
实用价值: 明确QMC的偏差机制有助于改进模拟方法，提高量子退火算法的实用性
商业应用: D-Wave等商用量子退火设备的性能评估和优化依赖于准确的模拟工具

现有方法局限性

薛定谔方程: 提供最忠实的QA动力学描述，但数值求解仅限于小系统
传统QMC: 虽然能处理大系统，但与QA在动力学本质上的差异可能导致行为偏差
缺乏定量分析: 之前研究缺乏对QMC重现QA不公平采样能力的系统性定量评估

核心贡献

发现了QMC的均匀性偏差: 首次系统性地证明了QMC倾向于产生比QA更均匀的基态采样
揭示了偏差的物理机制: 识别出副本对齐是导致均匀性偏差的根本原因
提供了定量评估框架: 建立了比较QMC和QA基态分布的量化指标
分析了Trotter数的影响: 阐明了离散化参数对采样公平性的影响机制
比较了不同转移规则: 研究了Metropolis和heat-bath方法对模拟准确性的影响

方法详解

任务定义

输入: 具有简并基态的Ising模型输出: 各基态的采样概率分布目标: 比较QMC和QA在基态采样上的差异

量子退火模型

使用横向磁场量子退火，时间依赖哈密顿量为： $\hat{H}(t) = \frac{t}{\tau}\hat{H}_0(\{\hat{\sigma}^z_i\}) - \left(1-\frac{t}{\tau}\right)\sum_{i=1}^N \hat{\sigma}^x_i$

其中 $\hat{H}_0 = -\sum_{i=1}^{N-1}\sigma_i\sigma_{i+1} - \sigma_1 + \sigma_N$ 是目标问题哈密顿量。

QMC实现

Suzuki-Trotter分解

将量子系统映射为经典系统，有效哈密顿量为： $H_{eff}(\sigma) = \frac{s}{M}\sum_{k=1}^M H_0(\sigma^k) - J^* \sum_{k=1}^M \sum_{i=1}^N \sigma_{i,k}\sigma_{i,k+1}$

其中 $J^* = \frac{1}{2\beta}\ln\coth\left(\frac{\beta}{M}(1-s)\right)$ ， $M$ 是Trotter数。

主方程动力学

系统按照离散时间主方程演化： $P(\sigma, t+\Delta t) = \left[1-\sum_{\sigma'(\neq\sigma)}w_{\sigma\to\sigma'}\Delta t\right]P(\sigma,t) + \sum_{\sigma'(\neq\sigma)}w_{\sigma'\to\sigma}\Delta t P(\sigma',t)$