Variational Real-Time Dynamics on Reduced Operator Manifolds

기술 요약: 축소된 연산자 매니폴드 상의 변분 실시간 역학

문제 정의
상관 관계가 있는 양자 다체계(quantum many-body systems)를 정확하게 실시간 시뮬레이션하는 것은 고전 알고리즘과 근미래 양자 하드웨어 모두에 있어 중요한 과제로 남아 있습니다. 텐서 네트워크 및 신경 양자 상태(neural quantum states)와 같은 고전적 방법들은 얽힘 성장(entanglement growth) 또는 확률적 샘플링 비용과 관련된 병목 현상에 직면해 있습니다. 양자 하드웨어 상에서 표준적인 접근 방식인 트로터화(Trotterization)는 시뮬레이션 시간 및 정확도에 따라 불리하게 스케일링되는 회로 깊이를 요구합니다. 변분 양자 실시간 진화(VQRTE)는 매개변수화된 회로 안사츠(ansatz)로 역학을 제한함으로써 절충안을 제시하지만, 맥라클란 변분 원리(McLachlan variational principle)에 기반한 기존 구현 방식은 양자 기하 텐서(QGT)의 추정을 필요로 합니다. 이 과정은 파라미터 수의 제곱에 비례하는 상당한 측정 오버헤드를 발생시키며, 노이즈 및 불량 조건(ill-conditioning)에 민가합니다. 또한, 기존 방법들은 시뮬레이션 정확도와 필요한 측정 자원 사이의 트레이드오프를 조절하는 데 어려움을 겪는 경우가 많습니다.

방법론: 연산자 투영 변분 양자 실시간 진화 (OVQRTE)
저자들은 매개변수화된 양자 회로를 업데이트할 때, 전체 상태 진화를 안사츠 접공간(tangent space)에 투영하는 대신 선택된 다항식 크기의 관측량 부분 집합에 대해 에렌페스트 방정식(Ehrenfest equations)을 강제하는 방법인 연산자 투영 변분 양자 실시간 진화(OVQRTE)를 소개합니다.

• 핵심 원리: OVQRTE는 전체 상태 역학을 구하는 대신, 일련의 기대값 $S = \{\hat{O}_\alpha\}$ (예: 국소 파울리 문자열 또는 축소 밀도 행렬 연산자)을 추적합니다. 파라미터 업데이트 $\dot{\theta}$는 이 관측량들의 변분 변화율과 정확한 에렌페스트 변화율 사이의 잔차를 최소화함으로써 결정됩니다:
  $$\dot{\theta} = \arg \min_{v} \sum_{\alpha} \left| \dot{x}^{var}_{\alpha} - i\langle [\hat{H}, \hat{O}_\alpha] \rangle \right|^2$$
  이는 선형 시스템 $A\dot{\theta} = C$를 산출하며, 여기서 $A$는 관측량 자코비안(Jacobian)의 그람 행렬(Gram matrix)이고 $C$는 교환자 기대값들을 포함합니다.
• 자원 효율성: 결정적으로, OVQRTE는 QGT를 명시적으로 구축하는 것을 피합니다. 행렬 $A$는 오직 기대값과 그 미분값들만으로 구성되며, 이는 큐비트 단위 가환(QWC) 그룹을 사용하여 동시에 측정될 수 있습니다. 이는 측정 비용을 (QGT 기반 방식의 파라미터 수의 제곱에서) 선형 또는 동시 섭동 확률적 근사(SPSA)를 사용할 경우 파라미터 수와 무관한 수준으로 감소시킵니다.
• 연산자 풀 축소: 하드웨어 노이즈와 회로 깊이를 더욱 완화하기 위해, 저자들은 "큐비트 들뜸 기반"(Qubit-Excitation-Based, QEB) 변형 방식을 제안합니다. 이 접근법은 페르미온-큐비트 매핑(예: Jordan-Wigner 변환)을 단순화하여, 비국소적 패리티 문자열(Pauli-Z 체인)을 항등 연산자로 대체하고 끝단 연산자만을 유지합니다. 이는 측정의 국소성과 파울리 문자열의 가중치를 크게 줄여 측정의 견고성을 높입니다.
• QSCI를 통한 바닥 상태 추출: 논문은 OVQRTE가 양자 선택적 구성 상호작용(Quantum-Selected Configuration Interaction, QSCI)의 생성기로 기능할 수 있음을 보여줍니다. 시간 진화 궤적에서 계산 기저 상태들을 샘플링함으로써, 압축되고 물리적으로 유의미한 부분 공간을 구축합니다. 바닥 상태는 이 부분 공간 내에서 해밀토니안을 대각화함으로써 근사됩니다. 하드웨어에서의 노이즈로 인한 대칭성 위반을 처리하기 위해, 저자들은 자기 일관적 구성 회복(S-CORE) 절차를 사용하며, 반 채워진 안데르손 불순물 모델(AIM)의 경우 입자-홀(PH) 대칭 완성을 적용합니다.

주요 결과

• Heisenberg 모델 벤치마크: 1차원 하이젠베르크 모델 시뮬레이션 결과, OVQRTE는 표준 VQRTE에 비해 측정 비용을 실질적으로 줄여, 시간 단계당 요구되는 회로 수를 1~2 차수(order) 낮추었습니다. 상태 벡터 시뮬레이션에서 OVQRTE는 VQRTE와 대등한 수치 정밀도를 달acia했습니다. 유한 샷 샘플링(하드웨어 노이즈 시뮬레이션) 환경에서 OVQRTE는 우수한 견고성을 보였으며, 동일한 샷 예산(shot budget)에 대해 VQRTE보다 낮은 상태 충실도(infidelity)와 더 정확한 자화 역학을 유지했습니다.
• 하드웨어 상의 Anderson Impurity Model (AIM): 이 알고리즘은 단일 밴드 안데르손 불순물 모델을 시뮬레이션하기 위해 최대 24개의 큐비트를 사용하여 IQM Emerald 초전도 프로세서에서 벤치마킹되었습니다. 연구에는 장치의 네이티브 연결성에 SWAP 게이트 없이 매핑된 Local Unitary Cluster Jastrow (LUCJ) 안사츠가 사용되었습니다. 반복적인 호출 과정에서 발생하는 하드웨어 및 샘플링 노이즈의 누적으로 인한 통제되지 않은 피드백 불안정성을 피하기 위해, 저자들은 전체 실시간 궤적을 고충실도로 재현하는 대신, 고전적 gGut 솔루션으로부터 얻은 수렴된 AIM 인스턴스에 대해 OVQRTE+QSCI를 벤치마킹하는 데 집중했습니다.
• QSCI 통합: OVQRTE로 생성된 상태들은 QSCI를 위한 부분 공간을 구축하는 데 사용되었습니다. 결과적으로 도출된 바닥 상태 에너지 추정치와 상태 밀도는 정확한 고전적 참조값과 비교했을 때 높은 정확도를 보였습니다. 이는 전체 실시간 궤적이 높은 충실도로 재현되지 않더라도, 진화된 상태들이 정보가 풍부한 저에너지 부분 공간을 점유하고 있다면, OVQRTE가 지표면(ground-state) 특성 추출을 위한 유효한 도구임을 입증합니다.

의의 및 주장
본 논문은 OVQRTE가 근미래 양자 하드웨어에서 상관 관계가 있는 응집 물질계를 조사하기 위한 유망한 접근법임을 확립합니다. 주요 의의는 다음과 같습니다:

1. 체계적 트레이드오프: 이는 연산자 집합 $S$의 크기와 선택을 통해 정확도와 측정 비용 사이의 트레이드오프를 조절할 수 있는 제어 가능한 노브(knob)를 제공하여 자원 효율적인 시뮬레이션을 가능하게 합니다.
2. 노이즈 탄력성: QGT를 피하고 축소된 연산자 집합(QEB)을 활용함으로써, 표준 변분 실시간 진화가 겪는 노이즈와 불량 조건 문제에 덜 민감합니다.
3. 역학 이상의 효용성: 본 연구는 전체 역학 궤적을 정확하게 재현하는 것이 다운스트림 애플리케이션에 반드시 필수적인 것은 아님을 강조합니다. OVQRTE는 바닥 상태 추정(QSCI를 통해)을 위한 정보가 풍부한 부분 공간을 효과적으로 생성하며, 이는 현재 하드웨어 내의 gGut과 같은 임베딩 이론 내에서 불순물 문제를 해결하기 위한 실행 가능한 도구임을 보여줍니다.

저자들은 OVQRTE가 기존 변분 방식의 자원 병목 현상을 우회하여 축소된 연산자 매니폴드를 활용함으로써, 상관 계통을 시뮬레이션할 수 있는 실질적인 경로를 제공한다고 결론짓습니다.