Symmetry-Adapted Eigensolver from Real-Space Sampling of Tensor-Network-Bridged Quantum Circuits

기술 요약: 텐서 네트워크 브릿지 양자 회로의 실공간 샘플링으로부터 도출된 대칭 적응형 고유값 솔버(Symmetry-Adapted Eigensolver)

문제 정의
강상관 양자 다체 계(strongly correlated quantum many-body systems)를 고전 컴퓨터로 시뮬레이션하는 것은 힐베르트 공간의 지수적 증가로 인해 어려움을 겪고 있다. 행렬 곱 상태(Matrix Product States, MPS)와 같은 텐서 네트워크 방법은 효율적인 압축을 제공하지만, 고차원 또는 복잡한 시스템에서 얽힘(entanglement)을 정확하게 포착하려면 결합 차원(bond dimension)을 급격히 증가시켜야 하며, 이는 막대한 메모리 및 계산 비용을 초내로 이어진다. 반대로, 양자 회로를 병렬 서브 회로로 분할하는 양자 컴퓨팅 접근 방식은 상당한 얽힘을 가진 시스템을 충실하게 표현하는 데 어려움을 겪을 수 있다. 본 논문은 파동 함수의 구성(configuration) 샘플링을 변분 상태(variational state)의 정의로부터 분리함으로써 양자 상관관계를 효율적으로 표현하는 새로운 전략의 필요성을 다룬다.

방법론
저자들은 병렬 양자 회로와 MPS 백본(backbone)으로 구성된 복합 안사츠(composite ansatz)를 활용하되, 이 복합 객체를 변분 파동 함수가 아닌 엄격한 **샘플 생성기(sample generator)**로 취급하는 "대칭 적응형 고유값 솔버(Symmetry-Adapted Eigensolver)"를 제안한다.

1. 복합 생성기 ($\phi(a_{\text{real}})$):

   • 격자는 크기가 $b$인 $n_b$개의 블록으로 나뉜다.
   • 국소 얽힘(Local Entanglement): 각 블록 내에서는 양자 회로 $C_i$(또는 축약 밀도 행렬(RDM) 고유벡터로부터 유도된 등거리 사상(isometry))가 국소 얽힘을 포착한다.
   • 전역 상관관계(Global Correlations): MPS 백본 $B_i$가 이 블록들을 연결하여, 렌ormalized 기저 $\gamma_i$를 연결함으로써 실공간 진폭 $\phi(a_{\text{real}})$을 생성한다.
   • 샘플링: 실공간 구성 $a_{\text{real}}$은 조건부 축약 밀도 행렬을 통한 완벽한(직접) 샘플링을 사용하여 확률 분포 p(a_{\text{real}}) = |\phi(a_{\text{real}})|^^2로부터 직접 추출되며, 이는 마르코프 체인(Markov chain)의 필요성을 제거한다.

2. 대칭 적응형 변분 상태 ($\psi(a_{\text{symm}})$):

   • 실제 변분 상태는 전체 격자 대칭성(병진, 반사, 반전, 입자 수)으로부터 구축된 대칭 적응형 기저에서 정의된다.
   • 모듈러스 규칙(Modulus Rule): 대칭 진폭의 모듈러스는 전체 등가 클래스(equivalence class)의 실공간 구성들에 대해 누적된 총 확률 가중치의 제곱근으로 정의된다:
     $$|\psi(a_{\text{symm}})| = \sqrt{\sum_{a_{\text{real}} \in \{g a_{\text{repr}}\}} |\phi(a_{\text{real}})|^2}$$
   • 위상 규칙(Phase Rule): 대칭 진폭의 위상은 해당 클래스의 단일 대표 구성 $a_{\text{repr}}$으로부터만 상속된다:
     $$\arg \psi(a_{\text{symm}}) = \arg \phi(a_{\text{repr}})$$
   • 핵심 차이점: 전통적인 투영 방법에서는 대칭 관계에 있는 진폭들이 코히어런트하게 합쳐져 간섭(interfering)을 일으키는 것과 달리, 이 방법은 구성을 확률적으로 결합한다. 이 대칭 상태는 생성기의 군 투영(group projection)이 아니라, 클래스 전반의 집합이 MPS 백본이 부담할 필요가 없는 장거리, 대칭 강제 얽힘을 공급한다.

3. 최적화 및 기울기(Gradients):

   • 에너지 기대값과 기울기는 $a_{\text{real}}$을 샘플링하고 이를 $a_{\text{symm}}$으로 매핑함으로써 평가된다.
   • 기울기 계산은 모듈러스와 위상 규칙을 효과적인 진폭으로 융합하여, 미분이 생성기 $\phi$에 작용하는 동시에 국소 에너지는 대칭 기저에서 평가되도록 한다.

주요 기여

• 샘플링과 정의의 분리: 본 연구는 파동 함수를 실공간 생성기로부터 샘플링하지만 대칭 적응형 기저에서 정의함으로써, 생성된 구성들이 재가중치(reweighting) 없이 타겟 대칭 분포의 편향되지 않은 샘플이 되도록 하는 프레임워크를 도입했다.
• 확률 기반 집합: 진폭 대신 확률을 집합함으로써 대칭 클래스 내의 간섭 문제를 피하고, 장거리 얽힘을 변분 상태에 자연스럽게 통합하여 매우 작은 MPS 결합 차원을 사용할 수 있게 한다.
• 등거리 초기화(Isometric Initialization): 양자 회로는 정확히 풀 수 있는 더 작은 시스템으로부터 얻은 축약 밀도 행렬(RDM)의 지배적인 고유벡터들로 구성된 등거리 사상으로 구현된다. 이는 깊고 매개변수가 많은 변분 양자 회로를 사용하는 것을 우회한다.
• 직접 샘플링: 이 방법은 생성기의 직접 샘플링을 활용하여 MCMC 방법과 관련된 자기상관 시간(autocorrelation times)을 피한다.

결과
이 방법은 1차원 및 2차원에서의 좌절된(frustrated) 스핀-1/2 $J_1-J_2$ 반강자성 하이젠베르크 모델에 대해 병진, 반사, 반전 및 입자 수 대칭성을 포함하여 벤치마킹되었다.

• 1D 시스템: 64-사이트 주기적 체인에 대해, 결합 차원 $D$ (최대 $D=6$)에 대한 선형 외삽 후 에너지 당 절대 오차 $10^{-5}$를 달로 달성했다. 이는 표준적인 접근 방식에 비해 현저히 적은 변분 매개변수를 사용하여 달성된 결과이다.
• 2D 시스템: $6 \times 6$ 토러스에 대해, $D \le 6$에서 유사한 정확도($10^{-5}$ 오차)를 얻었다.
• 얽힘: 대칭 적응형 기저가 장거리 얽힘을 운반한다는 것이 입증되었다. $D=1$에서, 대칭 상태는 정확한 바닥 상태의 이분 얽힘 엔트로피(bipartite entanglement entropy)를 재현한 반면, 생성기를 직접 실공간 파동 함수로 사용했을 때는 분할 영역에 대해 곱(product) 상태로 남았다.

의의 및 주장
본 논문은 이 접근 방식이 실공간 양자 샘플링, 대칭 제약 측정, 그리고 텐서 네트워크 재구성이 시너지를 내는 일관된 프레임워크를 구축한다고 주장한다.

• 효율성: 주요 의의는 DMRG 수준의 정밀도(오차 $\sim 10^{-5}$)에 일반적인 DMRG가 필요로 하는 결합 차원(종종 수천 단위)보다 수십 배 작은 결합 차원($D \le 6$)으로 도달할 수 있다는 점이다. 이는 대칭 적응형 기저가 얽힘의 부담을 짊어짐으로써 MPS 백본이 나머지 부분만을 표현하도록 하기 때문이다.
• 하드웨어 호환성: 이 아키텍처는 분산 양자 하드웨어에 적합하다. 블록들은 클래식하게 연결되며(양자 프로세서 간의 양자 링크 불필요), 생성 단계는 전체 시스템 크기($N$)와 독립적인 큐비트 수를 사용하여 양자 하드웨어에서 직접 실행 가능하다.
• 한계: 저자들은 생성 단계가 하드웨어 준비가 되어 있음에도 불구하고, 에너지 추정치가 현재 진폭 쿼리(진폭의 비율)를 요구하며, 이는 양자 프로세서에서 직접 측정 가능한 기대값 형태가 아니라는 점을 언급했다. 이 추정치를 하드웨어 구현을 위해 재구성하는 것이 향후 과제로 식별되었다.

본 연구는 일반적인 부호 문제(sign problem)를 해결하거나 모든 텐서 네트워크 방법을 대체한다고 주장하는 것이 아니라, 대칭성이 지배적인 특징인 강상관 시스템을 위한 특정적이고 효율적인 전략을 제시한다.