Spectral Gaps via Imaginary Time

기술 요약: 허수 시간(Imaginary Time)을 통한 스펙트럼 갭 측정

문제 정의
스펙트럼 갭은 $\Delta E = E_1 - E_0$ (바닥 상태와 첫 번째 들뜬 상태 사이의 에너지 차이)로 정의되며, 형광 파장, 물질의 엑시톤 에너지 전달부터 단열 실행 시간(adiabatic run times) 및 상태 전이 충실도에 이르기까지 다양한 물리적 맥락에서 중요한 매개변수입니다. 이를 결정하는 것은 계산적으로 매우 까다롭습니다. 최악의 경우, 국소 해밀토니안(local Hamiltonian)이 갭을 갖는지 여부를 결정하는 것은 결정 불가능(undecidable)하며, 낮은 에너지 스펙트럼을 추정하는 것은 일반적으로 NP-hard 문제입니다. 기존의 수치적 방법들, 예를 들어 페르미온 양자 몬테카를로(부호 문제(sign problem)로 인해 어려움을 겪음)나 텐서 네트워크 방법(결합 차원(bond dimension)에 의해 제한됨)은 시스템 크기, 기저 크기, 측정 오버헤드 사이에 상당한 트레이드오프를 마주합니다. 허수 시간 진화(ITE)는 바닥 상태를 찾기 위한 강력한 고전적 도구이지만, 전체 상태 밀도를 분해하거나 양자 위상 추정(quantum phase estimation)을 수행하지 않고 특정 스펙트럼 갭을 추출하는 것은 여전히 구체적인 과제로 남아 있습니다.

방법론
저자들은 중첩된 교환자(nested commutators)의 기댓값 비를 통해 스펙트럼 갭 $\Delta E$를 직접 추정하는 방법을 제안합니다. 이는 허수 시간 궤적을 따라 평가된 두 기댓값의 비율을 이용합니다.

1. 이론적 프레임워크:
   유한 차원 해밀토니안 $H$와 스펙트럼 분해 $H = \sum E_n \Pi_n$가 있고, 바닥 상태($\Pi_0$)와 첫 번째 들뜬 상태($\Pi_1$) 부분 공간 모두에 지지 집합(support)을 가진 초기 상태 $|\phi_0\rangle$가 있다고 가정할 때, 시스템은 허수 시간 $\tau$ 동안 $|\psi(\tau)\rangle = e^{-\tau H}|\phi_0\rangle$로 진화합니다.
   본 방법은 바닥 상태와 첫 번째 들뜬 상태를 결합하는 "조정(coordinating)" 관측량 $O$를 활용합니다. 저자들은 $M$차 중첩 교환자 $[H, O]_M$을 정의합니다.
   핵심 결과(정리 1)는 $\tau \to \infty$일 때, 기댓값의 비가 에너지 갭의 제곱으로 수렴함을 보여줍니다:
   $$\frac{\langle [H, O]_{M+2} \rangle_\tau}{\langle [H, O]_M \rangle_\tau} = (E_0 - E_1)^2 + \mathcal{O}(e^{-\tau(E_2 - E_1)})$$
   이 추정량은 개별 고유벡터를 사용하는 대신 투영 형식(projector formulations)에 의존하기 때문에, 높은 에너지 준위의 축퇴(degeneracy)에 대해 불가지론적(agnostic)입니다.

2. 로그 추정량(Logarithmic Estimators):
   논문은 또한 $\ln|\langle [H, O]_M \rangle_\tau|$와 $\tau$ 사이의 기울기로 $E_0 + E_1$을 추출하고, 비율 피팅의 잔차(residual)로부터 두 번째 갭 $E_2 - E_1$을 추출하는 로그 추정량을 유도하여, 작은 수의 수치적 나눗셈이 불안정할 때 사용할 수 있는 대안을 제공합니다.

3. 양자 구현 (ITQDE):
   유니터리 진화가 기본이지만 허수 시간 진화는 기본이 아닌 양자 하드웨어에서 이를 구현하기 위해, 저자들은 허수 시간 양자 동역학 에뮬레이션(ITQDE)을 활용합니다. 이들은 허드바드-스트라토노비치 변환(Hubbard-Stratonovich transformation)을 사용하여 유클리드 커널 $e^{-\tau H^2}$를 실시간 유니터리의 가우시안 평균으로 표현합니다. 이 커널은 가우스-에르미트 쿼드러처(Gauss-Hermite quadrature)를 사용하여 이산화되며, 진화를 유니터리의 가중 합으로 근사합니다.
   이후 기댓값들은 $\langle \phi_0 | U_j^\dagger P_\alpha U_k | \phi_0 \rangle$ 형태의 오버랩의 실수 부분들의 가중 합으로 축소됩니다 (여기서 $P_\alpha$는 파울리 문자열). 이러한 오버랩은 표준 하다마드 테스트(Hadamard test)를 사용하여 측정할 수 있습니다. 저자들은 커널의 대칭성 덕분에 이러한 오버랩의 실수 부분만 필요하다고 언급했습니다.

주요 결과
저자들은 두 가지 표준 모델에 대해 QuSpin 라이브러리를 통한 클래식 시뮬레이션을 통해 검증했습니다:

• 1D 횡장 이징 모델 (TFIM): 시스템 크기 $L=4$에 대해, 추정된 갭의 상대 오차는 $\tau$에 따라 지수적으로 감소하다가($\sim 10^{-8}$에 도달), 첫 번째 들뜬 상태가 수치적 노이즈 아래로 억제됨에 따라 매우 큰 $\tau$에서 다시 증가했습니다.
• 1D 페르미-허드바드 모델: 반 채워진 체인($L=4$)에 대해 유사한 오차의 지수적 감소가 관찰되었으며(상대 오차 약 $\sim 10^{-4}$), 이 방법은 비자명한 스핀 및 전하 구조를 가진 상호작용하는 페르미온 환경을 성공적으로 처리했습니다.
• 양자 시뮬레이션: Qiskit의 Aer 백엔드와 $10^6$ 샷을 사용하여 TFIM에 대한 개념 증명 시뮬레이션이 수행되었습니다. 이 방법은 ITQDE 프레임워크와 호환되었으나, 달성된 최소 오차는 $\sim 10^{-1}$이었습니다. 저자들은 이를 "대역폭 트레이드오프(bandwidth tradeoff)" 때문이라고 설명했습니다. 즉, 높은 상태를 억제하기 위해 필요한 가우시안 필터가 $\tau$와 쿼드러처 차수 $K$가 세심하게 조정되지 않을 경우 관련 오버랩을 샷 노이즈 바닥(shot-noise floor)으로 몰아넣기 때문입니다.

의의 및 주장
본 논문은 스펙트럼 추정을 위한 "겸손하지만 진정한 구조적 이점"을 제공한다고 주장합니다.

• 단순성 및 범위: 이 방법은 개별 고유 상태를 분해하거나, 양자 위상 추정을 수행하거나, 전체 상태 밀도를 재구성하지 않고도 낮은 에너지 갭을 추정하는 방법을 제공합니다. 이는 스펙트럼 문제를 하나의 매개변수 흐름(one-parameter flow)을 따른 완화(relaxation) 문제로 전환합니다.
• 강건성: 이 추정량은 높은 에너지 준위의 축퇴에 대해 강건하며, 클래식 희소 선형 대수와 양자 합-항(sum-of-terms) 해밀토니안 모두에 적용 가능합니다.
• 겸손함: 저자들은 이 접근 방식이 로컬 해밀토니안 문제의 최악의 복잡성을 우회하지 않는다는 점을 명시적으로 밝힙니다. 이는 상세한 상태 밀도가 필요한 경우 전체 스펙트럼 재구성이나 위상 추정을 대체하기 위한 것이 아닙니다.
• 향-후 전망: 이 연구는 개념 증명 역할을 합니다. 저자들은 이 방법이 양자 하드웨어에서 경쟁력을 갖추기 위해서는 매개변수($\tau, K$)를 최적화하기 위한 이전 연구의 전체 대역폭 분석을 도입하고, 더 나은 조정 관측량을 설계하며, 더 측정 효율적인 오버랩 추정량을 사용해야 한다고 언급했습니다. 이 논문은 갭 중심의 관점을 완화율과 갭이 주요 관심사인 미분 방정식을 풀기 위한 유효 해밀토니안을 사용하는 광범위한 프로그램과 결합되는 도구로 포지셔닝합니다.