High-order splitting of non-unitary operators on quantum computers

논문 요약: 양자 컴퓨터를 위한 비유니터리 연산자의 고차 분할 기법

1. 문제 제기 (Problem)

• 비유니터리 동역학의 중요성: 물리 과정의 대부분은 마찰, 점성, 확산, 경계 유출 등의 메커니즘으로 인해 소산 (dissipation) 과 비가역성을 수반합니다. 이는 고전 과학 계산뿐만 아니라 양자 시스템의 환경 상호작용에서도 필수적입니다.
• 양자 시뮬레이션의 한계: 양자 컴퓨터는 본질적으로 가역적인 유니터리 (unitary) 동역학 (국소 해밀토니안의 시간 진화) 을 구현하도록 설계되어 있습니다. 반면, 소산적인 동역학은 비유니터리 연산자로 표현되며, 이를 양자 프로세서에서 직접 구현하는 것은 어렵습니다.
• 기존 방법의 결함:
  • 고차 정확도를 얻기 위한 기존 곱 공식 (product formulas, 예: Strang splitting 등) 은 2 차 이상의 고차 근사를 위해 **음의 시간 계수 (negative time steps)**를 포함합니다.
  • 유니터리 동역학에서는 역시간 진화가 가능하지만, 소산적인 (비유니터리) 동역학에서는 역시간 진화가 수치적으로 불안정하며, 이를 유니터리 연산자로 블록 인코딩 (block-encoding) 하려면 재스케일링이 필요합니다.
  • 이로 인해 기존에는 소산적 동역학 시뮬레이션이 2 차 이하의 정확도로만 제한되었습니다.

2. 방법론 (Methodology)

저자들은 복소수 계수 (complex-coefficient) 곱 공식을 도입하여 고차 정확도를 유지하면서 안정적인 양자 알고리즘을 제안합니다.

• 연산자 분해:
  • 선형 시스템 $\frac{d\vec{\psi}}{dt} = M\vec{\psi}$에서 생성자 $M$을 에르미트 성분 ($H_1$, 소산적) 과 반에르미트 성분 ($iH_2$, 유니터리) 으로 분해합니다 ($M = H_1 + iH_2$).
  • $H_1$은 확산/감쇠를, $H_2$는 파동 전파/이류를 나타냅니다.
• 복소수 계수 곱 공식 적용:
  • 소산적 진화 ($e^{H_1 t}$) 에 대해 실수부가 양수인 복소수 계수 ($a_j$, $\Re(a_j) > 0$) 를 사용합니다. 이는 허수 시간 진화를 의미하며, 이는 유니터리 연산자로 변환 가능합니다.
  • 유니터리 진화 ($e^{iH_2 t}$) 에 대해 양수인 실수 계수 ($b_j > 0$) 를 사용합니다.
  • 각 단계 $S_j$는 다음과 같이 정의됩니다:
    $$S_j(\Delta t) = e^{iH_2 b_j \Delta t} e^{iH_1 \Im(a_j) \Delta t} e^{H_1 \Re(a_j) \Delta t}$$
  • 여기서 $e^{H_1 \Re(a_j) \Delta t}$는 감쇠 (수축) 를, $e^{iH_1 \Im(a_j) \Delta t}$는 추가적인 유니터리 위상 회전을 나타냅니다.
• 구체적 공식:
  • 4 차 공식 (Castella et al.): 5 단계로 구성되며, 계수 $a_j$와 $b_j$가 모두 양의 실수부를 가집니다.
  • 6 차 공식 (Bernier et al.): 16 단계로 구성되며, 대칭 켤레 (symmetric-conjugate) 구조를 가집니다.
• 회로 구현:
  • 파동 방정식 (감쇠된 파동) 을 예시로, 푸리에 공간에서 이산화된 모드를 다룹니다.
  • 유니터리 부분 (파동 전파) 은 제어된 $R_Y$ 게이트로, 소산 부분 (감쇠) 은 보조 큐비트 (ancilla) 와 후선택 (postselection) 을 통해 구현됩니다.
  • 허수 시간 진화 부분은 위상 게이트 (Phase gate) 로 변환되어 구현됩니다.

3. 주요 기여 (Key Contributions)

1. 고차 소산 시뮬레이션의 가능성 입증: 음의 시간 계수 없이, 실수부가 양수인 복소수 계수를 사용하여 소산적 동역학을 고차 (4 차, 6 차) 로 정확하게 분할할 수 있음을 이론적으로 증명했습니다.
2. 직접 블록 인코딩: 스케일링 없이 소산 진화 연산자를 직접 블록 인코딩하여 고차 정확도를 달성하는 방법을 제시했습니다.
3. 실제 하드웨어 검증: IonQ Forte 1(트랩드 이온 양자 컴퓨터) 에서 4 차 및 6 차 분할 알고리즘을 실행하여, 잡음이 있는 환경에서도 고차 방법이 유효함을 검증했습니다.

4. 실험 결과 (Results)

• 시뮬레이션 (Statevector):
  • 감쇠된 파동 방정식을 128 개의 그리드 포인트로 이산화하여 시뮬레이션했습니다.
  • 4 차 및 6 차 분할 공식은 이론적인 오차 스케일링 ($O(\Delta t^4)$, $O(\Delta t^6)$) 을 따랐으며, 낮은 오차 허용 범위 (high precision) 에서 1 차, 2 차 방법보다 훨씬 적은 CNOT 게이트 수로 동일한 정확도를 달성했습니다.
• 하드웨어 실행 (IonQ Forte 1):
  • 1 차 및 2 차: 낮은 정확도를 보였으며, 2 차는 4 차와 유사한 정확도를 보였습니다.
  • 4 차: 회로 깊이가 증가했음에도 불구하고 1 차 및 2 차보다 더 높은 정확도를 보였습니다. 성공 확률도 이상적인 수축 진화와 더 가까웠습니다.
  • 6 차: 회로 깊이가 너무 깊어져 잡음 (noise) 의 영향이 지배적이 되어, 오히려 1 차 결과보다 정확도가 떨어졌습니다.
  • 결론: 4 차 분할이 잡음이 있는 NISQ(Noisy Intermediate-Scale Quantum) 장치에서 정확도와 회로 깊이의 최적 균형을 이룹니다.

5. 의의 및 결론 (Significance)

• 실용적 접근법: 고차 분할 기법은 소산적 동역학 시뮬레이션에 있어 정확도와 효율성을 동시에 제공하는 실용적인 해결책입니다. 특히 4 차 방법이 잡음 환경에서도 우월한 성능을 입증했습니다.
• 자원 효율성: 고차 방법은 높은 정밀도가 요구될 때 CNOT 게이트 수를 기하급수적으로 줄여줍니다. 게이트 복잡도는 $O(2^{p/2} t^{1+1/p} \epsilon^{-1/p})$로 스케일링되어, 기존 2 차 방법 ($O(t^{3/2})$) 보다 시간 의존성이 훨씬 낮습니다.
• 미래 전망:
  • 6 차 이상의 계수를 찾는 연구가 필요하며, $b_j > 0$ 조건을 완화하면 더 높은 차수의 공식 발견이 쉬울 수 있습니다.
  • 이 방법은 열역학, 유체 역학 (확산 - 대류), 화학 반응 등 다양한 소산적 물리 현상을 양자 컴퓨터로 시뮬레이션하는 데 적용될 수 있습니다.

요약하자면, 이 논문은 복소수 계수를 활용한 고차 연산자 분할 기법을 통해 양자 컴퓨터에서 소산적 (비유니터리) 동역학을 안정적이고 정확하게 시뮬레이션할 수 있음을 이론 및 실험을 통해 입증하였으며, 특히 4 차 방법이 현재의 잡음 있는 양자 하드웨어에서 가장 유망한 접근법임을 보여주었습니다.