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

현실적인 홀수 차원 Gottesman-Kitaev-Preskill 상태를 이용한 회로의 고전 시뮬레이션

기본 정보

논문 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 함수의 음성(negativity)에 비례하므로, 특정 대규모 회로의 효율적인 시뮬레이션이 가능하다. 12dB 압축을 가진 안정화자 GKP 상태의 경우, 이 알고리즘은 최대 1000개 모드의 회로를 시뮬레이션할 수 있으며, 필요한 샘플 수는 단일 입력 모드에 필요한 샘플 수의 2배 미만이다. 이는 기존 시뮬레이터와 극명한 대조를 이룬다.

연구 배경 및 동기

문제 정의

보손 부호화 양자 회로의 고전 시뮬레이션은 많은 (잠재적으로 무한한) 에너지 준위를 가진 양자 시스템을 시뮬레이션해야 하기 때문에 거대한 도전에 직면해 있으며, 계산 비용이 극히 높다. 기존 시뮬레이션 알고리즘은 실제로 관련된 유한 압축 GKP 상태를 처리할 때 상당한 한계를 가진다.

중요성 분석

실용적 가치: GKP 부호화는 마이크로파 공동 결합 초전도 회로, 포획된 이온 및 광자 플랫폼 등의 실험 시스템에서 구현되며, 내결함성 양자 계산의 중요한 후보이다
이론적 의의: 초기 보손 양자 프로세서 구현의 검증을 위한 벤치마크 도구 제공
기술적 도전: 전통적 방법은 높은 압축 GKP 상태의 대규모 회로 시뮬레이션을 처리할 수 없다

기존 방법의 한계

Fock 기저 전개 방법: GKP 부호어가 무한히 많은 에너지 준위를 포함하기 때문에, 무차별적 시뮬레이션 방법은 소수의 GKP 큐비트에서도 실행 불가능해진다
준확률 분포 방법: 부호화된 GKP 상태의 Wigner 함수 음성이 크므로, 실행 시간이 폭발적으로 증가한다
가우스 상태 분해 방법: 실제로 관련된 높은 압축 GKP 상태에는 적용되지 않는다. 왜냐하면 이러한 상태들은 큰 비가우스성을 가지기 때문이다

핵심 기여

다중 모드 Zak-Gross Wigner (ZGW) 함수 제시: 단일 모드 ZGW 함수를 n-모드 경우로 일반화하고, 수정된 Stratonovich-Weyl 공리를 만족함을 증명
ZGW 함수와 Gross Wigner 함수 간의 관계 확립: CV GKP 상태의 ZGW 함수와 부호화된 논리 상태의 Gross Wigner 함수 간의 대응 관계 증명
효율적인 시뮬레이션 알고리즘 개발:
- 이상적인 안정화자 GKP 상태에 대한 약한 시뮬레이션 알고리즘
- 유한 압축 GKP 상태에 대한 확률 추정 알고리즘
현저한 성능 향상 달성: 12dB 압축을 가진 안정화자 GKP 상태의 경우, 1000개 모드의 회로 시뮬레이션은 단일 모드의 2배 미만의 샘플만 필요

방법론 상세 설명

작업 정의

다음 양자 회로의 측정 결과 확률 분포를 시뮬레이션:

입력: 역압축 매개변수 Δ를 가진 임의의 GKP 부호화 큐비트
진화: 정수 심플렉틱 행렬로 설명되는 유니터리 연산 및 임의의 변위 (모든 부호화 Clifford 연산 포함)
측정: 위치 기저의 모듈로 측정 $M_Z(s) = \frac{1}{d\ell}\sum_n e^{-i\ell s \cdot n}T_{(0,n)}$

ZGW 함수 정의

정의 1: 다중 모드 홀수 차원 ZGW 함수는 다음과 같이 정의된다: $W_{\hat{\rho}}(\eta) = \text{Tr}(\hat{\rho}\hat{A}_\eta)$

여기서 상점 연산자는: $\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$

이때 $\eta \in [0,d\ell)^{2n}$ , $d$ 는 양의 홀수, $\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_{\hat{\rho}}(\eta) = \bar{W}_{\bar{\rho}(u)}(t)$ 여기서 $\eta = \ell(u+t)$ , $u = \frac{1}{\ell}\eta \bmod 1$ , $t = \frac{1}{\ell}\eta - u$ 이다.

정리 1: 심플렉틱 진화 하에서의 ZGW 함수 변환 정수 심플렉틱 행렬 $S$ 로 설명되는 가우스 유니터리 진화에 대해, 진화된 상태의 ZGW 함수는: $W_{\hat{U}_S\hat{\rho}\hat{U}_S^\dagger}(\eta) = W_{\hat{\rho}}(S\eta - t)$

정리 2: 측정 확률의 계산 논리 연산자 $\hat{Z}_L^{\otimes n}$ 의 측정 확률 분포는: $\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 상태의 약한 시뮬레이션

유효 확률 분포 $W_{\hat{\rho}_0}(\eta)$ 에서 벡터 $\eta$ 샘플링
심플렉틱 행렬 $S$ 와 변위 벡터 $t$ 작용 하에서 변환: $S\eta - t$
선형 변위 연산자 적용: $S\eta - t + c$
측정 결과는 벡터의 후반부로 주어짐

유한 압축 GKP 상태의 확률 추정

준확률 분포 샘플링 방법 사용: $\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)|$