Quantum simulation of classical RCL circuits via complexification and Schrödingerization

기술 요약: 복소화 및 슈뢰딩거화를 통한 고전적 RCL 회로의 양자 시뮬레이션

문제 정의
동역학적 과정의 시뮬레이션은 물리학의 근간을 이루지만, 고전적 컴퓨터는 규모가 커짐에 따라 양자 시스템을 시뮬레이션하는 데 있어 본질적인 어려움에 직면한다. 반대로, 양자 컴퓨터는 전통적으로 고전적 기계가 수행하던 계산을 가속화할 수 있는 잠재력을 제공한다. 핵심적인 과제는 양자 역학(복소 힐베르트 공간에서 정의됨)과 고전 역학(실수 벡터 공간에서 정의됨) 사이의 엄격한 수학적 동등성을 확립하는 것이다. "탈복소화(decomplexification)" 절차는 양자 시스템을 고전적 유사체(예: 전기 공진 회로)로 매핑할 수 있게 해주지만, 그 역 과정인 **복소화(complexification)**는 고전적 미분 방정식을 양자 하드웨어에 적합한 등가 양자 형식으로 들어 올리기 위해 필요하다. 저자들은 특히 RCL(저항-인덕터-커패시터) 회로를 지배하는 구체적인 운동 방정식과 같이, 고전적 미분 방정식을 양자 형식으로 변환하기 위한 효율적인 방법을 구축해야 할 필요성을 다룬다.

방법론
본 논문은 고전 및 양자 역학 사이를 연결하기 위해 두 가지 뚜렷한 방법론적 접근 방식을 채택한다:

1. 대각화를 통한 복소화:
   저자들은 먼저 표준적인 복소화 절차를 재검토한다. 실수의 2차 미분 방정식 시스템(예: 감쇠 진동자 또는 결합된 PT-대칭 회로)에서 시작하여, 이 시스템을 1차 행렬 형태인 $\dot{\mathbf{x}} = A\mathbf{x}$로 변환한다. 실상태 행렬 $A$의 고유값(종종 복소 공액임)을 분석함으로써, 저자들은 역학을 실수 공간 $\mathbb{R}^{2n}$에서 복소 공간 $\mathbb{C}^n$으로 확장한다. 이는 대각화 행렬 $S$를 찾아 시스템을 분리된 복소 방정식으로 변환하는 과정을 포함한다. 여기에 $i$를 곱하면 슈뢰딩거 유형의 방정식 $i\partial_t z = H z$를 얻게 되며, 여기서 $H$는 일반적인 경우 비에르미트(non-Hermitian) 해밀토니안이다. 저자들은 이를 단일 감쇠 진동자와 한 쌍의 PT-대칭 결합 진동자에 대해 명시적으로 입증하고, 그에 대응하는 해밀토니안을 유도하며, 유사 변환(similarity transformations)을 통해 표준 PT-대칭 형태를 회복할 수 있음을 보여준다.

2. 슈뢰딩거화 (제안된 양자 알고리즘):
   대각화 방법이 대규모 시스템에서는 대수적으로 까다로워질 수 있음을 인지하여, 저자들은 "슈뢰딩거화(Schrödingerization)" 절차[14]를 채택한다. 이 방법은 일반적인 선형 미분 방정식 $\partial_t |u(t)\rangle = A |u(t)\rangle$를 보조 변수 $p$를 도입하여 직접적인 유니타리 슈뢰딩거 진화로 매핑한다.

   • 분해: 생성자 행렬 $A$를 에르미트 부분($H_1$)과 반-에르미트 부분($H_2$)으로 분해한다: $A = H_1 + iH_2$.
   • 워프 변환(Warp Transformation): 확장된 힐베르트 공간 $\mathcal{H} = \mathcal{H}_I \otimes \mathcal{H}_p$ 내에서 새로운 상태 $|\omega(p, t)\rangle = e^{-p}|u(t)\rangle$를 정의한다. 여기서 $\mathcal{H}_p$는 보조 자유도를 나타낸다.
   • 해밀토니안 구성: 진화는 에르미트 해밀토니안 $H_{Sch} = H_1 \otimes P_p - H_2 \otimes I_p$에 의해 제어된다. 여기서 $P_p = -i\partial_p$는 운동량 연산자로 작용한다.
   • 양자 구현: 양자 컴퓨터에서 이를 구현하기 위해, 보조 변수 $p$를 유한 격자(finite grid) 상에 이산화하고 이를 양자 레지스터에 인코딩한다. 운동량 연산자 $P_p$는 양자 푸리에 변환(QFT)을 사용하여 구현된다. 진화는 운동량 기저(여기서 $P_p$는 대각 행렬임)에서 수행된 후 QFT를 통해 위치 기저로 다시 변환된다.
   • 회복(Recovery): 특정 회복 영역 내의 격자점 $p_j$를 선택하고 역 워프 인자 $e^{p_j}$를 적용함으로써 원래의 고전적 해를 회복한다.

주요 기여 및 결과

• 이론적 프레임워크: 본 논문은 탈복소화와 복소화 사이의 수학적 관계를 명확히 하며, 단일 감쇠 진동 및 PT-대칭 결합 시스템을 포함한 고전 RCL 회로에 대한 양자 등가 해밀토니안을 명시적으로 유도한다.
• 알고리즘 구현: 저자들은 슈뢰딩거화 방법을 양자 컴퓨터에 구현하기 위한 상세한 청사진을 제공한다. 여기에는 보조 변수의 이산화, 운동량 기저 해밀토니안의 구성, 그리고 구체적인 회로 시퀀스(상태 준비 $\to$ 역 QFT $\to$ $H_{mom}$ 하의 유니타리 진화 $\to$ QFT $\to$ 회복)가 포함된다.
• 수치적 검증: 이 방법은 단일 감쇠 조화 진동기($\omega=2, \gamma=1$)에 대해 $N_p=2048$개의 격자점(보조 큐비트 11개 + 시스템 큐비트 1개)을 사용하는 양자 시뮬레이터를 통해 수치적으로 테스트되었다.
  • 정확도: 재구성된 해 $u_{rec}(t)$는 대부분의 시뮬레이션 시간 동안 위치($x$) 및 속도($y$) 성분에 대해 정확한 고전적 해 $u_{exact}(t)$와 매우 잘 일치하였다.
  • 오차 분석: 정규화된 방향 오차는 초기 및 중간 진화 단계에서는 작게 유지되었다. 그러나 후기 시간에는 오차가 증가했다. 저자들은 이를 회복 과정 때문이라고 설명한다. $H_1$이 양의 고유값을 가지므로, 해는 큰 양의 $p$에서 회복되어야 한다. 이는 물리적 진폭이 이미 0에 가까워진 상태에서 작게 지수적으로 감소하는 보조 프로파일로 나누는 과정을 수반하며, 이 과정이 수치적 이산화 및 푸리에 오차를 증폭시킨다.

의의 및 주장
본 논문은 고전적 미분 방정식을 양자 하드웨어가 직접 다룰 수 있는 형태로 재구성하는 것이 고전 시스템을 효율적으로 시뮬레이션하는 데 있어 매우 중요하다는 점을 시사한다. 저자들은 직접적인 복소화 방법이 이론적으로는 타당하지만, 대수적 복잡성으로 인해 대규모 시스템에는 종종 부적합할 수 있다고 주장한다. 반면, 슈뢰딩거화 절차는 고전 시스템을 슈뢰딩거 방정식에 의해 제어되는 양자 역학으로 매핑하는 단순하고 일반화 가능한 경로를 제공한다.

연구는 결론적으로, QFT 기반의 슈뢰딩거화 접근 방식이 양자 컴퓨터에서 고전 미분 방정식을 풀기 위한 유망한 도구라고 밝힌다. 감쇠 진동기의 역학을 성공적으로 수치적으로 재현한 것은, 보조 운동량 변수와 회복 영역을 다룰 때 후기 시간의 수치적 민감도를 완화하기 위한 주의가 필요하다는 전제하에 이 방법의 실현 가능성을 입증한다. 본 논문은 모든 고전 시스템에 대한 일반적인 시뮬레이션 문제를 해결했다고 주장하는 것이 아니라, RCL 회로와 같은 선형 시스템에 대한 실행 가능한 프레임워크를 구축하는 것을 목적으로 한다.