Functional matrix product state simulation of continuous variable quantum circuits

1. 연구 배경 및 문제 제기 (Problem)

• 연속 변수 (CV) 양자 컴퓨팅의 중요성: 광자나 진동 모드를 활용하여 큐비트를 인코딩하는 CV 양자 시스템은 GKP (Gottesman-Kitaev-Preskill) 상태나 고양이 (Cat) 상태와 같은 '보손 큐비트'를 통해 높은 내결함성 (fault-tolerance) 을 제공할 수 있어 차세대 양자 컴퓨팅 플랫폼으로 주목받고 있습니다.
• 비고스 (Non-Gaussian) 시스템의 시뮬레이션 난제: CV 시스템에서 범용 양자 계산을 위해서는 비고스 연산이 필수적이지만, 이를 수치적으로 시뮬레이션하는 것은 매우 어렵습니다.
  • 기존 방법들 (포크 공간 기반 텐서 네트워크, 가우스 근사 등) 은 비고스 상태가 포함된 회로를 시뮬레이션할 때 계산 비용이 기하급수적으로 증가하거나 정확도가 떨어지는 한계가 있습니다.
  • 특히, 높은 압축 (squeezing) 이나 비고스 입력 상태 (예: GKP 상태) 를 다룰 때 기존 기술들은 확장성 (scaling) 이 부족합니다.

2. 제안된 방법론 (Methodology)

저자들은 함수 행렬 곱 상태 (Functional Matrix Product State, FMPS) 라는 새로운 수치 기법을 제안했습니다. 이는 이산 변수 시스템의 MPS 를 연속 변수 시스템으로 자연스럽게 확장한 것입니다.

• 실공간 표현 (Real-space Representation):
  • CV 상태의 파동 함수 $f(q)$를 실공간 (position space) 에서 직접 표현합니다.
  • 연속 슈미트 분해 (Continuous Schmidt decomposition) 를 기반으로 파동 함수를 텐서 네트워크 형태로 분해합니다.
  • 수식적으로 $f(q) \approx \sum \gamma_1 \cdot \gamma_2 \cdots \gamma_m$ 형태로 표현되며, 여기서 $\gamma$는 연속 변수 $q$와 이산 인덱스 (bond index) 를 가진 함수입니다.
• 이산화 및 보정 (Discretisation & Bounding Box):
  • 무한한 실공간을 유한한 격자 (grid) 로 이산화합니다.
  • 각 게이트 연산 (변위, 압축, 위상 회전, 빔 스플리터 등) 에 따라 상태의 지지 영역 (support) 이 변하므로, 계산 영역 (bounding box) 을 동적으로 조정합니다.
  • 격자 점 사이의 값을 보간하기 위해 3 차 스플라인 (cubic splines) 을 사용하여 파동 함수의 매끄러움과 대칭성을 보존합니다.
• 오차 관리:
  • 경계 상자 절단 (Truncation to bounding box): 상태가 거의 0 이 되는 영역을 잘라냅니다.
  • 이산화 (Discretisation): 격자 수 $N$을 적절히 설정하여 오차를 최소화합니다.
  • MPS 차수 절단 (Bond dimension truncation): 특이값 (singular values) 의 가중치를 기반으로 결합 차수 (bond dimension) 를 제한하여 계산 효율성을 높입니다.
• 노이즈 모델링:
  • 광자 손실 (photon loss) 과 같은 노이즈를 시뮬레이션하기 위해, 선형 광학 소자와 손실 채널이 교환 가능 (commute) 한 성질을 이용합니다. 이를 통해 회로 끝에서 한 번에 노이즈 채널을 적용하여 계산 비용을 줄입니다.

3. 주요 기여 (Key Contributions)

1. FMPS 프레임워크 개발: CV 양자 회로를 실공간에서 효율적으로 시뮬레이션할 수 있는 새로운 텐서 네트워크 방법론을 정립했습니다.
2. 범용 연산자 및 측정 구현: 변위, 압축, 위상 회전, 빔 스플리터, 3 차 위상 게이트 (cubic phase gate) 등 범용 CV 양자 계산을 위한 모든 기본 연산자와 광자 수 측정, 동위상 측정 (homodyne), 이위상 측정 (heterodyne) 등을 FMPS 형식에 맞게 구현했습니다.
3. 비고스 상태에서의 성능 입증: GKP 상태, 고양이 상태 등 기존 방법으로 시뮬레이션하기 어려웠던 고도로 비고스적인 입력 상태를 가진 회로에서 기존 방법 (Strawberry Fields 라이브러리의 Bosonic backend) 보다 월등히 우수한 성능을 보였습니다.
4. 확장성 분석: 얕은 회로 (shallow circuits) 에서는 결합 차수가 포화되어 모드 수에 대해 부분 지수적 (sub-exponential) 으로 확장됨을 보였습니다.

4. 실험 결과 (Results)

저자들은 무작위 얕은 회로와 캐스케이드 (cascaded) 회로를 시뮬레이션하여 FMPS 의 성능을 검증했습니다.

• 캐스케이드 회로 (Cascaded Circuits):
  • 입력 상태가 GKP, 고양이, 압축 상태일 때, FMPS 는 모드 수 (modes) 가 증가함에 따라 계산 시간이 부분 지수적으로 증가했습니다.
  • 반면, 기존 방법 (Strawberry Fields) 은 비고스 입력 상태의 경우 모드 수에 대해 지수적으로 계산 시간이 증가하여 메모리 부족으로 시뮬레이션이 중단되었습니다.
• 넓은 회로 (Wide Circuits):
  • 1 차원 배열에 랜덤 빔 스플리터 레이어를 적용한 회로에서도 FMPS 는 GKP 상태와 같은 복잡한 입력에 대해 기존 방법보다 우월한 성능을 보였습니다.
• 노이즈 시뮬레이션:
  • 광자 손실을 포함한 시뮬레이션에서도 FMPS 는 확장성을 유지하며 동작했습니다. 노이즈를 추가하면 계산 시간이 증가하지만, 여전히 기존 방법보다 효율적이었습니다.
• 정확도:
  • 격자 수 ($N \approx 200$) 와 결합 차수를 적절히 설정하면, 해석적 해나 기존 시뮬레이션 결과와 매우 높은 일치도 (quadrature distance $\epsilon \sim 10^{-6}$) 를 보였습니다.

5. 의의 및 결론 (Significance)

• 기술적 돌파구: CV 양자 컴퓨팅, 특히 GKP 큐비트와 같은 비고스 자원을 다루는 소규모 양자 프로토콜의 특성을 분석하고 설계하는 데 필수적인 도구를 제공합니다.
• 효율성: 기존 방법론이 처리하지 못했던 고도로 비고스적인 상태와 얕은 회로를 정확하게 시뮬레이션할 수 있게 하여, 실험적 결과의 검증 및 새로운 양자 알고리즘 개발을 가속화할 것으로 기대됩니다.
• 미래 전망: 이 연구는 CV 양자 시스템의 오류 정정 및 확장 가능한 아키텍처 설계에 중요한 기초를 제공하며, 향후 더 복잡한 CV 양자 회로 시뮬레이션의 표준 기법으로 자리 잡을 가능성이 큽니다.

요약하자면, 이 논문은 실공간 기반의 함수 행렬 곱 상태 (FMPS) 를 도입하여 비고스 연속 변수 양자 회로의 시뮬레이션 난제를 해결하고, 기존 방법론보다 확장성과 정확도를 획기적으로 개선한 획기적인 연구입니다.