2025-11-12T10:37:09.722784

Classical simulation of circuits with realistic odd-dimensional Gottesman-Kitaev-Preskill states

Calcluth, Hahn, Bermejo-Vega et al.
Classically simulating circuits with bosonic codes is challenging due to the prohibitive cost of simulating quantum systems with many, possibly infinite, energy levels. We propose an algorithm to simulate circuits with encoded Gottesman-Kitaev-Preskill (GKP) states, specifically for odd-dimensional encoded qudits. Our approach is tailored to be especially effective in the most challenging but practically relevant regime, where the codeword states exhibit high (but finite) squeezing. Our algorithm leverages the Zak-Gross Wigner function introduced by J. Davis et al. [arXiv:2407.18394], which represents infinitely squeezed encoded stabilizer states positively. The runtime of the algorithm scales with the negativity of the Wigner function, allowing for efficient simulation of certain large-scale circuits - namely, input stabilizer GKP states undergoing generalized GKP-encoded Clifford operations followed by modular measurements - with a high degree of squeezing. For stabilizer GKP states exhibiting 12 dB of squeezing, our algorithm can simulate circuits with up to 1,000 modes with less than double the number of samples required for a single input mode, in stark contrast to existing simulators. Therefore, this approach holds significant potential for benchmarking early implementations of quantum computing architectures utilizing bosonic codes.
academic

Classical simulation of circuits with realistic odd-dimensional Gottesman-Kitaev-Preskill states

基本信息

  • 论文ID: 2412.13136
  • 标题: Classical simulation of circuits with realistic odd-dimensional Gottesman-Kitaev-Preskill states
  • 作者: Cameron Calcluth, Oliver Hahn, Juani Bermejo-Vega, Alessandro Ferraro, Giulia Ferrini
  • 分类: quant-ph (量子物理)
  • 发表时间: 2024年12月17日 (arXiv预印本)
  • 论文链接: https://arxiv.org/abs/2412.13136

摘要

本文提出了一种经典模拟算法，用于模拟基于Gottesman-Kitaev-Preskill (GKP)编码的奇维量子电路。该算法特别针对具有高压缩度但有限压缩的码字状态这一最具挑战性但实际相关的场景进行了优化。算法利用Zak-Gross Wigner函数，该函数能够正表示无限压缩的编码稳定子状态。算法的运行时间与Wigner函数的负性呈比例关系，使得某些大规模电路的高效模拟成为可能。对于具有12dB压缩的稳定子GKP状态，该算法能够模拟多达1000个模式的电路，所需样本数不到单输入模式所需样本数的两倍，这与现有模拟器形成鲜明对比。

研究背景与动机

问题定义

经典模拟玻色编码量子电路面临巨大挑战，主要原因是需要模拟具有许多（可能无限）能级的量子系统，计算成本极其昂贵。现有的模拟算法在处理实际相关的有限压缩GKP状态时存在显著局限性。

重要性分析

  1. 实用价值: GKP编码在微波腔耦合超导电路、囚禁离子和光子平台等实验系统中得到实现，是容错量子计算的重要候选方案
  2. 理论意义: 为验证早期玻色量子处理器的实现提供基准测试工具
  3. 技术挑战: 传统方法无法处理高压缩GKP状态的大规模电路模拟

现有方法局限性

  1. Fock基展开方法: 由于GKP码字涉及无限多能级，暴力模拟方法对少数几个GKP量子比特就变得不可行
  2. 准概率分布方法: 编码GKP状态的Wigner函数负性很大，导致运行时间爆炸
  3. 高斯态分解方法: 不适用于实际相关的高压缩GKP状态，因为这些状态具有大的非高斯性

核心贡献

  1. 提出了多模Zak-Gross Wigner (ZGW)函数: 将单模ZGW函数推广到n模情况，并证明其满足修正的Stratonovich-Weyl公理
  2. 建立了ZGW函数与Gross Wigner函数的关系: 证明了CV GKP状态的ZGW函数与编码逻辑状态的Gross Wigner函数之间的对应关系
  3. 开发了高效模拟算法:
    • 对于理想稳定子GKP状态提供弱模拟算法
    • 对于有限压缩GKP状态提供概率估计算法
  4. 实现了显著的性能提升: 对于12dB压缩的稳定子GKP状态，模拟1000个模式的电路仅需不到单模式两倍的样本数

方法详解

任务定义

模拟如下量子电路的测量结果概率分布：

  • 输入: 任意GKP编码的量子比特，具有反压缩参数Δ
  • 演化: 由整数辛矩阵描述的酉操作和任意位移（包括所有编码Clifford操作）
  • 测量: 位置基的模测量 MZ(s)=1dℓ∑ne−iℓs⋅nT(0,n)M_Z(s) = \frac{1}{d\ell}\sum_n e^{-i\ell s \cdot n}T_{(0,n)}

ZGW函数定义

定义1: 多模奇维ZGW函数定义为： Wρ^(η)=Tr(ρ^A^η)W_{\hat{\rho}}(\eta) = \text{Tr}(\hat{\rho}\hat{A}_\eta)

其中相点算子为： A^η=1(2π)n∑a∈Z2neiℓ[a,η]+iπaXTaZT^a\hat{A}_\eta = \frac{1}{(2\pi)^n}\sum_{a \in \mathbb{Z}^{2n}} e^{i\ell[a,\eta]+i\pi a_X^T a_Z}\hat{T}_a

这里 η∈[0,dℓ)2n\eta \in [0,d\ell)^{2n}，dd是正奇数，T^a=eiπaXTaZ/dT^aXT^aZ\hat{T}_a = e^{i\pi a_X^T a_Z/d}\hat{T}_{a_X}\hat{T}_{a_Z}。

关键理论结果

引理1: ZGW函数与Gross Wigner函数的关系 Wρ^(η)=Wˉρˉ(u)(t)W_{\hat{\rho}}(\eta) = \bar{W}_{\bar{\rho}(u)}(t) 其中 η=ℓ(u+t)\eta = \ell(u+t)，u=1ℓη mod 1u = \frac{1}{\ell}\eta \bmod 1，t=1ℓη−ut = \frac{1}{\ell}\eta - u。

定理1: 辛演化下的ZGW函数变换 对于整数辛矩阵SS描述的高斯酉演化，演化后状态的ZGW函数为： WU^Sρ^U^S†(η)=Wρ^(Sη−t)W_{\hat{U}_S\hat{\rho}\hat{U}_S^\dagger}(\eta) = W_{\hat{\rho}}(S\eta - t)

定理2: 测量概率的计算 逻辑算子 Z^L⊗n\hat{Z}_L^{\otimes n} 的测量概率分布为： Tr(ρ^M^Z(s))=∫dηXWρ^((ηXs))\text{Tr}(\hat{\rho}\hat{M}_Z(s)) = \int d\eta_X W_{\hat{\rho}}\left(\begin{pmatrix}\eta_X\\s\end{pmatrix}\right)

模拟算法

理想稳定子GKP状态的弱模拟

  1. 从有效概率分布 Wρ^0(η)W_{\hat{\rho}_0}(\eta) 中采样向量 η\eta
  2. 在辛矩阵SS和位移向量tt作用下变换：Sη−tS\eta - t
  3. 应用线性位移算子：Sη−t+cS\eta - t + c
  4. 测量结果由向量的后半部分给出

有限压缩GKP状态的概率估计

使用准概率分布采样方法： Pr(η)=1MU^S,cρ^0U^S,c†∣WU^S,cρ^0U^S,c†(η)∣\text{Pr}(\eta) = \frac{1}{M_{\hat{U}_{S,c}\hat{\rho}_0\hat{U}_{S,c}^\dagger}}|W_{\hat{U}_{S,c}\hat{\rho}_0\hat{U}_{S,c}^\dagger}(\eta)|

所需样本数为： N=2ϵ2Mρ^02log⁡(2/δ)N = \frac{2}{\epsilon^2}M_{\hat{\rho}_0}^2\log(2/\delta)

其中Mρ^M_{\hat{\rho}}是ZGW函数的负性（负体积）。

实验设置

理论分析框架

本文主要进行理论分析和数值计算，重点关注：

  1. ZGW函数负性的计算: 对于不同压缩水平的GKP状态
  2. 样本复杂度分析: 基于负性的缩放行为
  3. 与现有方法的比较: 理论复杂度对比

评价指标

  • ZGW对数负性: log⁡Mρ^\log M_{\hat{\rho}}，衡量模拟复杂度
  • 样本数缩放: 相对于单模式情况的样本数增长
  • 模式数容量: 可有效模拟的最大模式数

实验结果

有限压缩GKP状态的ZGW函数

通过显式计算，得到现实0-逻辑GKP状态的ZGW函数： WGKP((uv))∝ϑ(Γ;z)W_{\text{GKP}}\left(\begin{pmatrix}u\\v\end{pmatrix}\right) \propto \vartheta(\Gamma; z)

其中z=(v/(dℓ),−u/(dℓ),0,0)Tz = (v/(d\ell), -u/(d\ell), 0, 0)^T，Γ\Gamma是由压缩参数Δ\Delta确定的复矩阵。

主要数值结果

对于Δ=0.25\Delta = 0.25（对应12dB压缩）：

  • 单模负性: Mρ^0Δ≈e3×10−4M_{\hat{\rho}_0^\Delta} \approx e^{3 \times 10^{-4}}
  • 千模式模拟: 由于负性的乘性性质，模拟1000个输入模式仅需不到单模式两倍的样本数
  • 性能提升: 相比现有方法提升数个数量级

魔态分析

对于π状态 ∣ψπ⟩=13(∣0L⟩+∣1L⟩−∣2L⟩)|\psi_\pi\rangle = \frac{1}{\sqrt{3}}(|0_L\rangle + |1_L\rangle - |2_L\rangle)：

  • 在所有Δ\Delta值下都表现出显著负性
  • 模拟大规模电路时开销显著，符合预期

相关工作

经典模拟算法发展

  1. Gottesman-Knill定理: 稳定子电路的多项式时间模拟
  2. 准概率方法: 基于Gross Wigner函数的采样算法
  3. 玻色系统模拟: 利用准概率分布正性、张量网络、恒星表示等方法

GKP编码模拟现状

现有方法主要限于：

  • 理想（无限压缩）GKP状态
  • 高斯操作和零差测量
  • 无法处理实际相关的有限压缩情况

本文创新点

相比现有工作，本文首次实现了有限压缩GKP状态的高效模拟，特别是在高压缩regime下的显著性能提升。

结论与讨论

主要结论

  1. 理论贡献: 建立了ZGW函数的完整理论框架，包括多模推广和演化规律
  2. 算法创新: 提出了适用于现实GKP状态的高效模拟算法
  3. 实用价值: 为验证玻色量子处理器提供了强有力的基准测试工具

局限性

  1. 奇维限制: 当前方法仅适用于奇维GKP量子比特，扩展到偶维（如量子比特）仍是开放问题
  2. 操作限制: 限于整数辛矩阵描述的操作，无法处理任意单模压缩
  3. 噪声模型: 未考虑损耗等实际噪声对ZGW函数负性的影响

未来方向

  1. 扩展到GKP量子比特: 开发适用于二维情况的方法
  2. 更广泛的玻色编码: 将框架推广到GKP之外的其他玻色编码
  3. 准概率框架: 发展专门针对玻色编码的全面准概率理论

深度评价

优点

  1. 理论严谨性: 数学推导完整，证明详细，理论基础扎实
  2. 实用性强: 解决了实际相关的有限压缩GKP状态模拟问题
  3. 性能突出: 在高压缩regime下实现了数量级的性能提升
  4. 方法创新: ZGW函数的多模推广和应用是重要的理论贡献

不足

  1. 适用范围有限: 仅适用于奇维系统，限制了普适性
  2. 实验验证缺乏: 主要是理论工作，缺少与实际量子系统的对比验证
  3. 复杂度分析: 对于大规模电路的实际计算复杂度分析不够详细

影响力

  1. 学术价值: 为量子模拟理论提供了新的工具和视角
  2. 实用前景: 对早期容错量子计算实验具有重要指导意义
  3. 技术推动: 可能启发更多玻色编码模拟方法的发展

适用场景

  • 玻色量子处理器的基准测试和验证
  • GKP编码容错量子计算方案的理论分析
  • 连续变量量子计算的经典模拟研究

参考文献

本文引用了68篇相关文献，主要包括：

  • Gottesman-Kitaev-Preskill编码的原始工作
  • 量子模拟算法的经典文献
  • 玻色量子计算的最新进展
  • Wigner函数和准概率分布的理论基础

总体评价: 这是一篇高质量的理论物理论文，在GKP编码量子电路的经典模拟方面取得了重要突破。虽然存在一些局限性，但其理论贡献和实用价值都很显著，对量子计算领域具有重要意义。