Floquet-ADAPT-VQE: A Quantum Algorithm to Simulate Non-Equilibrium Physics in Periodically Driven Systems

기술 요약: Floquet-ADAPT-VQE

문제 정의
주기적으로 구동되는(Floquet) 양자 시스템은 시공간 결정(time crystals)이나 Floquet 위상 구조와 같이 정적 시스템에서는 나타나지 않는 비평형 현상을 보여준다. 그러나 이러한 시스템을 시뮬레이션하는 것은 시간 진화(예: Trotterization)를 통해 발생하는 오차와 힐베르트 공간의 지수적 증가로 인해 평형 시스템을 시뮬레이션하는 것보다 훨씬 더 까다롭다. 디지털 양자 컴퓨터가 이러한 시뮬레이션을 위한 범용 플랫폼을 제공하지만, 현재의 접근 방식들은 다음과 같은 구체적인 한계에 직면해 있다:

1. 시간 영역 알고리즘은 종종 Trotterization과 반복적 양자 위상 추정(IQPE)을 포함하는 깊은 회로를 요구하며, 이는 짧은 시간 주기나 높은 주파수로 제한된다.
2. 주파수 영역 알고리즘(예: 변분 해밀토니안 안사츠)은 제한된 차원의 보조 공간에 의존하고, 임의의 큐비트 절단(truncation)에 대한 보조 행렬의 파울리 분해에 대한 일반적인 처방이 부족하며, 안사츠 초기화 문제로 인해 수렴 성능이 떨어진다.
3. 일반적 공백: 확장된 Floquet 해밀토니안에 대한 확립된 범용 파울리 문자열(Pauli-string) 분해법이 없으므로, Floquet-Hilbert 공간 형식 내에서 양자 알고리즘을 체계적으로 구현하는 데 어려움이 있다. 또한, 기존 방법들은 깊은 회로 없이 Floquet 상태에 대한 관측량의 기댓값을 계산하는 프레임워크가 부족하다.

방법론
저자들은 Floquet-Hilbert 공간 형식을 사용하여 주기적으로 구동되는 시스템을 시불변(time-independent) 형태로 매핑하도록 설계된 하이브리드 양자-고전 알고리즘인 Floquet-ADAPT-VQE를 제안한다.

1. Floquet-Hilbert 공간 매핑:

   • 시간 주기적 해밀토니안 $\hat{H}(t)$는 물리적 힐베르트 공간 $\mathcal{H}$와 푸리에 모드(Floquet zone)를 나타내는 보조 공간 $\mathcal{A}$의 텐서 곱 공간 $\mathcal{F} = \mathcal{A} \otimes \mathcal{H}$ 상에서 작용하는 유효한 시불변 확장 Floquet 해밀토니안 $\hat{\hat{H}}_F$로 변환된다.
   • 저자들은 보조 행렬(대각 및 비대각 성분)을 큐비트 기반 프로세서에 네이티브한 파울리 문자열로 분해하는 체계적인 방법을 제공한다. 이들은 해당 행렬들이 보조 큐비트 수 $n_a$에 대해 전체 $4^{n_a}$가 아닌 $O(2^{n_a})$로 스케일링되는 파울리 문자열의 부분 집합으로 분해될 수 있음을 입증한다.
   • 절단 후에도 해밀토니안이 에르미트(Hermitian) 성과 Floquet zone 간의 대칭성을 유지하도록 하기 위해, 비대칭 보정 항($A^{asy}$)을 도입하고 이에 대한 명시적인 파울리 분해를 제공한다.

2. 알고리즘 설계 (Floquet-ADAPT-VQE):

   • 비용 함수 (Cost Function): $\hat{\hat{H}}_F$의 기댓값을 직접 최소화하는 대신(이는 절단 오차가 가장 큰 최하위 Floquet zone을 타겟팅하게 됨), 알고리즘은 비용 함수 $C_\lambda(\vec{\theta}) = \langle \psi(\vec{\theta}) | (\hat{\hat{H}}_F - \lambda \hat{\mathbb{I}})^2 | \psi(\vec{\theta}) \rangle$를 최소화한다. $\lambda$를 중앙 Floquet zone($-\Omega/2 < \lambda < \Omega/2$) 내로 선택함으로써, 알고리즘은 절단 오차가 최소화되는 상태를 타겟팅한다.
   • 초기 상태: 핵심적인 기여는 특정 보조 초기 상태 $|0\rangle$(대각 보조 행렬 $A_d$의 고유값이 0인 고유 상태)을 식별한 것이다. 저자들은 이 초기 상태의 제곱 Floquet 에너지가 보조 큐비트 수($n_a$) 및 구동 주파수($\Omega$)에 독립적임을 증명하고 수치적으로 입증하였다. 이는 초기 상태가 중앙 Floquet zone에 가깝게 유지되도록 보장하여, 무작위 초기 상태에 비해 수렴 성능을 크게 향상시킨다.
   • 연산자 풀 (Operator Pool): 그래디언트 측정 비용을 관리하기 위해 보조 공간에 대한 파울리 연산자와 물리적 공간에 대한 문제 특화 연산자(예: 2-local 연산자)를 결합한 "혼합 Floquet 풀"을 구성한다.

3. 관측량 및 역학 (Observables and Dynamics):

   • 이 프레임워크는 고정 깊이 양자 회로를 사용하여 시간 의존적 기댓값을 계산하고 시간 주기적 Floquet 상태 $|\phi(t)\rangle$를 준비할 수 있게 한다.
   • 보조 큐비트에 시간 의존적 $Z$-회전을 적용하고 특정 기저에서 측정함으로써, 원래의 힐베르트 공간에서의 Trotterization과 달리 회로 깊이가 진화 시간에 따라 증가하지 않고도 마이크로모션(micromotion) 역학을 추출할 수 있다.

주요 결과
저자들은 외부 자기장이 인가된 주기적 구동 XYZ 스핀 체인의 고전적 수치 시뮬레이션을 통해 알고리즘을 검증하였다:

• 수렴성: 시뮬레이션 결과, 제안된 초기 상태($|0\rangle \otimes |\phi_p\rangle$)는 $n_a$나 $\Omega$가 증가함에 따라 성능이 저하되는 무작위 초기 상태와 달리, 제곱 Floquet 에너지가 $n_a$ 또는 $\Omega$에 따라 증가하지 않음을 확인하였다.
• 정확도: 알고리즘은 다양한 구동 파라미터(다양한 구동 주파수 및 진폭 포함)에 대해 정확한 대각화(exact diagonalization, $\hat{U}_F$를 통한) 결과와 일치하는 준에너지(quasienergies) 및 Floquet 고유 상태를 성공적으로 계산하였다.
• 스펙트럼 재구성: 시프트 파라미터 $\lambda$를 조정함으로써, 알고리즘은 중앙 Floquet zone 내의 전체 준에너지 스펙트럼을 성공적으로 재구성하였다.
• 관측량: 알고리즘은 순 자화(net magnetization) 및 2점 상관 함수(two-point correlation functions)의 기댓값을 성공적으로 계산하였다.

의의 및 주장
본 논문은 다음과 같은 연구 문헌의 근본적인 공백을 해결한다고 주장한다:

1. 범용 프레임워크: 확장된 Floquet 해밀토니안의 체계적인 파울리 분해를 제공함으로써, Floquet-Hilbert 공간 내에서 NISQ 및 결함 허용(fault-tolerant) 양자 알고리즘 모두의 구현을 가능하게 한다.
2. 효율성: 준에너지를 계산하기 위해 Trotterization과 IQPE를 피하는 접근 방식을 제공하여, 시간 영역 방법보다 더 얕은 회로로 임의의 구동 주파수에서 시스템을 시뮬레이션할 수 있게 한다.
3. 강건성(Robustness): 보조 공간의 크기 및 구동 주파수와 독립적으로 수렴을 보장하는 특정 초기 상태 전략을 제공하며, 이는 기존 변분 접근법의 알려진 병목 현상이다.
4. 역학 접근성: 고정 깊이 회로를 사용하여 마이크로모션 역학 및 시간 주기적 상태에 접근할 수 있는 독특한 능력을 제공하며, 이는 스트로보스코픽(stroboscopic) 역학만을 접근하거나 깊은 시간 진화 회로를 요구하는 기존 방법들과 차별화된다.

저자들은 Floquet-ADAPT-VQE를 현재 및 미래의 양자 하드웨어에서 주기적으로 구동되는 다체 시스템의 시공간 결정 및 빛-물질 상호작용과 같은 비평형 현상을 연구하기 위한 기초적인 도구로 제시한다.