Low-Depth Initial-State Preparation for Ground-State Energy Estimation of Two-Dimensional Strongly Correlated Systems

기술 요약: 2차원 강상관 계의 바닥 상태 에너지 추정을 위한 저심도 초기 상태 준비

문제 정의
양자 위상 추정(QPE)과 같은 바닥 상태 에너지 추정용 양자 알고리즘은 실제 바닥 상태와 비영(non-negligible)의 충실도(fidelity)를 갖는 초기 상태를 필요로 한다. 바닥 상태 에너지를 성공적으로 추출할 확률은 이 충실도($F_0$)에 따라 결정된다. 만약 $F_0$가 시스템 크기에 따라 지수적으로 감소한다면, 필요한 반복 횟수($1/F_0$)로 인해 알고리즘은 실행 불가능해진다. 변분 양자 고유값 솔버(VQE)는 하이브리드 접근 방식을 제공하지만, 기울기(gradient) 계산 비용과 얕은 회로의 클래식 시뮬레이션 가능성 문제에 직면한다. 또한, 정확한 클래식 근사치(예: 행렬 곱 상태, MPS)를 양자 회로로 컴파일하는 과정은 과도한 회로 심도를 초래할 수 있으며, 이는 논리 큐비트와 매직 상태(magic state)가 제한적인 초기 결함 허용 양자 컴퓨터(FTQC)에 부담이 된다. 핵심 과제는 깊은 최적화 루프나 클래식의 난해함에 의존하지 않고, 2차원(2D) 허바드 모델과 같이 강상관 계를 위한 고충실도의 초기 상태를 준비할 수 있는 얕고 자원 효율적인 양자 회로를 설계하는 것이다.

방법론
저자들은 2D 하이젠베르크 모델과 반 채움(half-filled) 허바드 모델 사이의 슈리퍼-울프(Schrieffer–Wolff, SW) 변환을 통한 강결합 연결성을 활용하여 얕은 양자 회로를 구축하는 전략을 제안한다. 방법론은 세 가지 주요 단계로 진행된다:

1. 하이젠베르크 회로 구축:

   • 안사츠(Ansatz): 회로는 반강자성 하이젠베르크 모델의 $S=0$ (싱글렛) 부공간을 목표로 한다. 초기 상태는 격자의 완전 매칭(perfect matching) 상에 정의된 싱글렛 다이머(singlet dimer)들의 곱($| \Phi \rangle$)이다.
   • 역학(Dynamics): 상태는 지수화된 SWAP($e$SWAP) 게이트 층을 사용하여 진화한다. 이 게이트들은 SU(2) 대칭성을 준수하여 상태가 싱글렛 부공간 내에 머물도록 보장한다.
   • 구조: 격자 결합은 네 개의 서로 소인 부분 집합($B_{x0}, B_{x1}, B_{y0}, B_{y1}$)으로 분할된다. 회로는 여러 블록($N_b$)에 걸쳐 특정 순서($B_{y1} \to B_{y0} \to B_{x1} \to B_{x0}$)로 $e$SWAP 게이트를 적용한다. 이 구조는 하이젠베르크 바닥 상태의 공명 결합 가치(resonating-valence-bond, RVB) 물리학을 모사한다.
   • 페르미온 매핑: 스핀 상태는 보조 큐비트를 도입하는 국소 유니터리 임베딩($V$)을 통해 단일 점유 페르미온 부공간으로 매핑된다.

2. SW 게이트를 통한 허바드 확장:

   • 유한한 $U/t$에서의 허바드 모델에 존재하는 전하 요동(더블론 및 홀론)을 포착하기 위해, 하이젠베르크 회로는 전하 요동 게이트로 "드레싱(dressing)"된다.
   • 이 게이트들은 1차 SW 생성자($\Omega$)로부터 직접 유도되며, 두 사이트 유니터리($W_{\langle i,j \rangle}$)로의 스즈키-트로터(Suzuki–Trotter) 분해를 통해 근사된다.
   • 최종 변분 안사츠는 회전 각도 $\theta$와 $\phi$로 매개변수화된 하이젠베르크 진화($U_B$)와 전하 요동 블록($W_B$)을 결합한다.

3. 충실도 추정 및 벤치마킹:

   • 매개변수 전이: 작은 $4 \times 4$ 격자에서 최적화된 매개변수는 추가 최적화 없이 더 큰 격자($6 \times 6$ ~ $10 \times 10$)로 전이된다.
   • 시뮬레이션: 회로 상태는 믿음 전파(belief propagation, BP) 가이드 단순 업데이트를 사용하는 투영 얽힘 쌍 상태(PEPS)를 사용하여 시뮬레이션된다.
   • 충실도 경계: 대규모 시스템에 대한 정확한 바닥 상태를 알 수 없으므로, 저자들은 바닥 상태 충실도의 하한을 추정한다. 이는 다음을 통해 이루어진다:
     • 회로 상태와 변분 몬테카를로(VMC) 참조 상태 사이의 충실도($F_{ref}$) 계산.
     • 참조 에너지, 정확한 바닥 상태 에너지(Stochastic Series Expansion으로부터 도출), 그리고 싱글렛 들뜸 간극($\Delta_0$)을 사용하여 VMC 상태의 정확한 바닥 상태에 대한 충실도($F_{ref,0}$)를 구속.
     • 푸비니-스터디(Fubini–Fubini) 각도에 대한 삼각형 부등식을 적용하여 $F_0$에 대한 엄격한 하한을 도출.

주요 기여 및 결과

• 효과적인 매개변수 전이: 하이젠베르크 모델의 경우, $4 \times 4$ 격자에서 최적화된 매개변수는 재최적화 없이 최대 $10 \times 10$ 격자에 대해서도 견고한 충실도 하한을 제공한다. 구체적으로 $N_b=3$일 때, 추정된 바닥 상태 충실도 하한은 $10 \times 10$ 격자에서도 양수($\approx 0.028$)를 유지하며, 이는 초기 상태가 진폭 증폭(amplitude amplification)에 적합함을 시사한다.
• 허바드 충실도 개선: $4 \times 4$ 격자에 대한 정확한 계산은 SW 유도 전하 요동 게이트를 포함하는 것이 허바드 바닥 상태에 대한 충실도를 크게 향상시킨다는 것을 보여준다. $N_b=4$일 때, 충실도는 하이젠베르크만 있을 때의 $\approx 0.44$에서 최적화된 전하 각도를 사용할 때 $\approx 0.89$로 증가한다.
• 자원 효율성: 저자들은 $10 \times 10$ 격자에 대한 상태 준비 비용을 추정한다. 회전당 오차를 $10^{-5}$로 가정하는 Clifford+$T$ 합성을 적용할 때, 총 비용은 약 $10^5$ $T$ 게이트 수준이다. 이는 상태 준비 비용을 초기 FTQC의 "메가퀽프(megaquop)" 영역 내에 위치시키며, 후속 바닥 상태 에너지 추정 알고리즘의 비용보다 현저히 낮다.
• 확장성: 회로 심도는 격자 결합 수($|B|$)에 따라 선형적으로 스케일링되며, 매개변수 공유 전략을 통해 작은 시스템의 최적화를 기반으로 큰 시스템을 위한 회로를 구축할 수 있다.

의의 및 주장
본 논문은 강상관 계의 바닥 상태 에너지 추정을 위한 저비용 초기 상태 준비의 실행 가능한 경로를 제시한다고 주장한다. 하이젠베르크 모델과 허바드 모델 사이의 물리적 연결을 활용함으로써, 저자들은 얕은 회로가 정확한 클래식 대각화가 닿지 않는 규모의 격자에서도 비영의 바닥 상태 충실도를 달erm할 수 있음을 입증했다. 이러한 결과는 얕은 회로 심도와 $T$-게이트 수를 최소화하는 것이 중요한 초기 FTQC 시대에 이 회로들이 매력적인 후보임을 시사한다. 이 연구는 매개변수가 전이 가능하고 회로 심도가 얕게 유지된다면, 작은 격자 시뮬레이션이 클래식하게 난해한 더 큰 시스템을 위한 양자 회로 설계를 효과적으로 가이드할 수 있다는 가설을 뒷받침한다. 저자들은 자신들의 경계값이 엄격하지만 실제 충실도는 더 높을 수 있음을 명시하며, 향-후 연구에서는 더 깊은 해밀토니안 변분 안사츠나 더 정확한 시뮬레이션 기법(예: full-update PEPS)을 통해 이러한 추정치를 더욱 정교화할 수 있다고 언급했다.