Realistic simulation of universal multiport interferometers with losses within a Gaussian formalism

제목: 가우시안 형식론을 이용한 손실이 포함된 범용 멀티포트 간섭계의 현실적인 시뮬레이션

문제 정의
광학 기술은 양자 우위를 입증하기 위한 주요 후보 기술이지만, 광자 손실은 결함 허용(fault-tolerant) 광 양자 컴퓨터를 실현하는 데 있어 가장 큰 장애물입니다. 이러한 손실을 정확하게 특성화하는 것은 노이즈가 있는 백엔드에서 양자 알고리즘을 클래식하게 시뮬레이션하는 데 매우 중요합니다. 현재의 시뮬레이션 접근 방식은 종종 광학적 손실이 모든 채널에 걸쳐 균일하다는 사전 가설에 의존합니다. 저자들은 이러한 가정이 특정 하드웨어에서의 알고리즘 결과에 대해 잘못된 예측을 초래할 수 있다고 주장합니다. 또한, 큐비트 기반 장치에는 오류 정정 코드가 존재하지만, 연속 변수(continuous-variable, CV) 양자 컴퓨팅에서 광자 손실을 보정하기 위한 확립된 기술은 아직 존재하지 않습니다. 노이즈 환경에서의 CV 시뮬레이션을 위한 기존 소프트웨어 라이브러리는 물질 기반 큐비트(예: Qiskit, Cirq)에 비해 부족한 실정입니다. 또한, 모든 상호작용마다 외부 환경 모드를 도입하여 손실을 모델링하는 방식은 모드 수의 다항식적 증가를 초래하여, 포크 공간(Fock space) 시뮬레이션을 계산적으로 불가능하게 만듭니다.

방법론
저자들은 공간 빈(spatial-bin) 및 시간 빈(time-bin) 인코딩 아키텍처 모두에 적용 가능한, 멀티포트 간섭계의 광학적 손실을 특성화하고 시뮬레이션하기 위한 가우시안 형식론을 제안합니다. 방법론은 다음과 같이 진행됩니다:

1. 이론적 프레임워크: 시스템은 포크 기저(Fock basis) 대신 하이젠베르크 기저(quadratures $\hat{q}$ 및 $\hat{p}$)를 사용하여 기술됩니다. 가우시안 상태의 역학은 평균 벡터와 공분산 행렬에 의해 지배됩니다. 광학적 손실은 물리적 모드와 외부 환경 모드 사이의 상호작용을 빔 스플리터 연산으로 표현하여 모델링됩니다.
2. 불완전한 간섭계의 특성화: 저자들은 Qi 등[52]에 의해 기술된 불완전한 마하-젠더 간섭계(MZI)의 특성화를 바탕으로 연구를 구축합니다. 이들은 입력, 출력 및 빔 스플리터 자체 내부의 손실 요인을 포함하여 손실이 있는 간섭계 $T_\eta$를 모델링합니다.
3. 특이값 분해 (SVD): 수많은 외부 환경 모드를 추적하는 데 따른 계산적 폭발을 피하기 위해, 저자들은 SVD를 사용합니다. 불완전한 간섭계 메쉬로부터 발생하는 비유니터리 연산자 $A$는 $A = U_2 \sqrt{H} U_1$로 분해되며, 여기서 $U_1$과 $U_2$는 유니터리 행렬이고 $\sqrt{H}$는 소산 계수(singular values)를 포함하는 대각 행렬입니다. 이는 모든 손실 채널을 $M$개의 외부 모드(여기서 $M$은 물리적 모드의 수)로 응축하여, 시스템 크기를 다항식적 증가가 아닌 관리 가능한 $2M \times 2M$ 공분산 행리로 줄여줍니다.
4. 공분산 행렬 유도: 저자들은 이러한 손실을 겪는 가우시안 보존 샘플링(GBS) 과정의 공분산 행렬 $\sigma_S$에 대한 해석적 표현을 유도합니다. 이 행렬은 전역 심플렉틱 변환으로부터 환경 모드들을 부분적으로 트레이스 아웃(trace out)함으로써 얻어집니다.
5. 수치 시뮬레이션: 저자들은 GBS를 시뮬레이션하기 위해 가우시안 백엔드(Python 라이브러리 PNN)를 사용하는 파이프라인을 구현합니다. 이들은 다음을 비교합니다:
   • 이상적인 GBS 결과.
   • 유도된 비균일 손실 행렬 $\sqrt{H}$를 사용하는 현실적인 GBS 결과.
   • 균일한 손실($\sqrt{H} \propto I$)을 가정하는 GBS 결과.
     또한, 간섭계 메쉬의 Clements 및 Reck 분해 방식을 테스트합니다.

주요 기여

• 해석적 형식론: 본 논문은 SVD를 통해 유도된, 불완전한 간섭계로 수행되는 GBS 과정의 공분산 행렬에 대한 폐쇄형 표현식(식 24)을 제공합니다. 이를 통해 포크 공간을 거치지 않고도 노이즈가 있는 GBS 및 일반적인 다인(dyne) 기반 알고리즘을 클래식 백엔드에서 시뮬레이션할 수 있습니다.
• 효율적인 모드 감소: SVD를 활용함으로써, 저자들은 손실을 모델링하는 데 필요한 보조 모드의 수를 회로 깊이에 대한 다항식 의존성에서 선형 의존성($2M$)으로 줄여 계산 비용을 크게 낮추었습니다.
• 균일성 가정의 검증: 본 연구는 균일한 광학 손실을 가정하는 것의 타당성을 수치적으로 평가합니다. 연구 결과, 이러한 가정은 간섭계 분해 방식(Clements 대 Reck)과 모드의 수에 따라 달라지므로 보편적으로 참이 아님을 입증했습니다.
• 장치 특화적 특성화: 이 프레임워크는 개별 MZI의 특정 손실률과 같은 특정 하드웨어의 결함에 관한 실험 데이터를 시뮬레이션 파이프라인에 직접 통합할 수 있게 해줍니다.

결과

• 분포 변화: 수치 시뮬레이션 결과, 광학적 손실은 진공 상태의 인구(population)를 증가시키고, 이상적인 스퀴즈드 상태 입력에서는 금지되었던 홀수 광자 수를 검출할 수 있게 하여, 이상적인 경우에 비해 확률 분포를 크게 변화시킵 direction합니다.
• 분해 방식에 따른 의존성: 연구는 손실 행렬 $\sqrt{H}$가 항등 행렬(identity matrix)로부터 벗어나는 정도(즉, 비균일성)가 아키텍처에 따라 다르게 나타남을 보여줍니다:
  • Reck 방식(루프 기반 시간 빈 인코딩과 동등)의 경우, 항등 행렬로부터의 거리는 모드 수 $M$이 증가함에 따라 증가합니다. 저자들은 이 경향을 멱법칙(power law)과 로그 함수에 피팅하였으며, 멱법칙이 약간 더 잘 맞는다는 것을 발견했습니다.
  • Clements 방식(공간 인코딩)의 경우, 거리는 진동하며 $M=5$ 이후 묶인 범위(bounded range)로 수렴하는데, 이는 균일성으로부터의 편차가 더 안정적임을 시사합니다.
• Hafnian 발산: 행렬의 하프니안(Hafnian) 샘플링을 시뮬레이션할 때, 저자들은 현실적인 비균일 손실 모델이 균일 손실 근사보다 이상적인 경우로부터 더 빠르게 발산하는 것을 관찰했습니다. 이러한 차이는 광학적 손실이 증가할수록(투과율 $\eta$가 낮아질수록) 그리고 모드 수 $M$이 커질수록 더욱 뚜렷해집니다.

의의 및 주장
저자들은 본 연구가 향arian 특정 장치 및 아키텍처의 미래적인 특성화를 위한 "청사진"을 제공한다고 주장합니다. 제안된 프레임워크는 노이즈가 있는 중규모 양자(NISQ) 광학 플랫폼에서 양자 알고리즘을 현실적으로 시뮬레이션할 수 있게 합니다. 비균일한 손실을 정확하게 모델링함으로써, 이 방법은 균일한 손실 분포를 가정함으로써 발생할 수 있는 알고리즘 결과에 대한 잘못된 예측을 방지하는 것을 목표로 합니다. 저자들은 이 작업이 오류 완화 기술과 결려 현실적인 노이즈 모델을 통합하고, 잠재적으로 연속 변수 컴퓨팅을 위한 미래의 오류 정정 코드를 개발하기 위한 기초적인 단계라고 명시했습니다. 그들은 이 연구가 새로운 오류 정정 코드를 제안하는 것이 아니라, 노이즈 환경을 더 잘 이해하고 특성화하기 위한 시뮬레이션 도구임을 분명히 하고 있습니다.