Closest Accessible Symmetry reduction: a tool for Hamiltonian interpolation analysis

기술 요약: 해밀토니안 보간 분석을 위한 가장 가까운 접근 가능한 대칭성 축소 (Closest Accessible Symmetry Reduction)

문제 정의
해밀토니안 보간(Hamiltonian interpolations)의 분석, 특히 양자 단열 계산(AQC) 및 양자 상전이 맥락에서의 분석은 일반적으로 보간 매개변수 $s$를 이산화하고 순간적인 대상들을 연구하는 방식(예: 정확한 대각화, 텐서 네트워크 또는 기하학적 접근법 활용)에 의존한다. 이는 볼록 보간 $H(s) = (1-s)A + sB$에서 전체 스펙트럼 정보가 이론적으로 초기($A$) 및 최종($B$) 해밀토니안의 결합 구조에 의해 결정됨에도 불구하고, 중간 단계의 해밀토니안 구조를 점 단위로 분석하게 되는 개념적 불균형을 초래한다. 기존의 슈리퍼-울프(Schrieffer-Wolff, SW) 변환과 같은 방법들은 직접적인 회전이 저에너지 섹터들을 연결하는 섭동 영역(perturbative regimes)에 국한되어 있으며, 이러한 회전이 정의되지 않거나 큰 오차를 유발하는 임계 영역을 가로지르는 보간에서는 실패한다. 본 논문은 시간 이산화나 섭동 가정에 과도하게 의존하지 않고, $A$와 $B$의 결합 구조로부터 보간에 대한 전역적 스펙트럼 정보를 추출할 수 있는 방법론의 필요성을 다룬다.

방법론: 가장 가까운 접근 가능한 대칭성 (CAS)
저자들은 **가장 가까운 접근 가능한 대칭성(Closest Accessible Symmetry, CAS)**에 기반한 프레임워크를 도입하여, 힐베르트 공간을 약하게 결합된 유사 고유 공간(pseudo-eigenspaces)으로 재귀적으로 분해한다.

1. 핵심 개념: 정확한 대칭성을 찾는 대신(이는 계산적으로 불가능함), 제한되고 다루기 쉬운 인증 가능한 반사(힐베르트 공간의 이분할) 제품군 내에서 "가장 가까운" 대칭성을 탐색한다.
2. 최적화: 해밀토니안 $H$에 대하여, 비대각 결합 노름(off-diagonal coupling norm)을 최소화하는 투영 연산자 $P_0$ (반사 $T = 2P_0 - 1$을 정의함)를 찾는다:
   $$\min_{P_0 \in \mathcal{A}} \| P_0 H (1-P_0) + (1-P_0) H P_0 \|_2$$
   여기서 $\mathcal{A}$는 문제의 구조로부터 도출된 접근 가능한 투영 연산자 제품군이다.
3. 접근 가능한 집합의 구성: 접근 가능한 집합 $\mathcal{A}$는 $A$와 $B$에 존재하는 파울리 항들의 동역학적 리 대수(Dynamical Lie Algebra, DLA) 분석을 사용하여 구성된다. 이 방법은 특정 항들이 제거된 수정된 해밀토니안과 교환 법칙이 성립하는 단일 파울리 문자열(또는 그들의 곱)로 표현 가능한 대칭성으로 탐색 범위를 제한한다. 이는 $Z_2$ 대칭성을 유도하기 위해 $A$와 $B$에서 제거해야 할 최소한의 항들을 효과적으로 식별한다.
4. 재귀적 분해: 이 과정은 재귀적으로 적용된다. 각 단계 $r$에서 해밀토니안은 식별된 CAS의 섹터들로 투영된다. 이를 통해 다음을 얻는다:
   • 유사 고유값 ($\mu_k$): 투영된 해밀토니안의 대각 블록들.
   • 잔여 결합 ($\nu$): 투영 과정에서 버려진 명시적인 비대각 항들.
   • 재귀는 $2 \times 2$ 행렬이 남을 때까지 계속되며, 이들은 $s$의 함수로서 정확하게 대각화된다.
5. 오차 경계 및 하이브리드화: 방법론은 유사 갭(pseudo-gap) $\tilde{\Delta}$와 잔여 결합의 노름 $\|\nu\|$를 사용하여 실제 스펙트럼 갭 $\Delta$에 대한 엄격한 경계를 제공한다. 저자들은 하이브리드화 매개변수 $\chi_{kj} = \|\nu_{kj}\|_2 / \tilde{\Delta}_{kj}$를 정의한다.
   • $\chi_{kj} < 1/2$인 경우, 유사 레벨들은 약하게 하이브리드화되어 있으며, 바일의 부등식(Weyl's inequalities)과 페쉬바흐-슈어 맵(Feshbach-Schur map)을 통해 타이트한 오차 경계를 제공한다.
   • $\chi_{kj} \ge 1/2$인 경우, 강하게 하이브리드화된 상태들의 클러스터가 식별된다. 저자들은 기저 상태가 약하게 결합된 경우 또는 기저 상태가 강하게 하이브리드화된 클러스터의 일부인 경우(임계점을 나타냄)에 대한 구체적인 경계를 도출한다.

주요 기여

• 전역적 스펙트럼 표현: 이 방법은 약하게 상호작으는 섹터들 사이의 결합을 반복적으로 절단함으로써, 로컬한 설명들의 집합이 아닌 보간 전체에 대한 단일한 계층적 설명을 제공한다.
• 분석적 상전이 진단: 프레임워크는 유사 스펙트럼(pseudo-spectrum)을 통해 서로 다른 유형의 양자 상전이를 구분한다:
  • 1차 상전이: 두 특정 레벨 사이에서 유사 갭이 닫히며, 나머지 스펙트럼에 대한 하이브리드화가 낮게 나타나는 특징을 가진다.
  • 연속(임계) 상전이: 들뜬 상태의 타워(tower)가 기저 상태로 붕괴하며, 거대한 강한 하이브리드화 클러스터를 형성하는 특징을 가진다. 이러한 전이를 매개하는 CAS는 매우 비국소적(non-local)일 것으로 예상된다.
• 계산 가능성: 탐색을 "접근 가능한" 대칭성(DLA 및 파울리 문자열로부터 도출됨)으로 제한함으로써, 지수적인 복잡도를 가진 정확한 대칭성을 찾는 것을 피하면서도 근사 오차에 대한 제어를 유지한다.
• 고전적 휴리스틱과의 연결: CAS 브랜치를 통한 재귀적 경로는 최적화 문제(예: Ising 모델)에서 저에너지 상태를 찾기 위한 고전적 휴리스틱으로 해석될 수 있다.

결과
저자들은 횡장(transverse-field) 해밀토니안과 이징 모델 사이의 단열 보간을 사용하여 프레임워크를 검증하였으며, 특히 **좌절된 이징 링(frustrated Ising ring, $N=101$)**에 대해 테스트하였다.

• 1차 상전이: 표준 CAS 축소(최대 비대각 노름 최소화)는 약 3.5%의 상대 오차로 1차 상전이(섭동적 반교차)의 위치를 성공적으로 예측했다. 또한 전이 근처에서 기저 상태, 제1 들뜬 상태, 제2 들뜬 상태의 하이브리드화를 정확히 식별했다.
• 연속 상전이: 표준 축소 방식은 오차가 해당 영역에서 최대가 되기 때문에 임계점($s_c$)에 대한 정량적 추정치를 제공하는 데 실패했다. 그러나 보간의 중간 지점에 편향을 주는 **가중 비용 함수(weighted cost function)**를 도입함으로써, 저자들은 매우 정확한 임계점 예측($s_c^{pred} = 0.5125 \pm 0.0025$ vs. true $0.5135 \pm 0.0005$)을 달ей했다.
• 스케일링: 분석 결과, 관련 비대각 결합이 나타나는 재귀 깊이는 시스템 크기에 따라 $r \approx (N+1)/2$로 스케일링됨을 확인하였으며, 이는 상태 간의 해밍 거리(Hamming distance)와 관련된 섭동적 반교차에서의 갭의 지수적 억제와 일치한다.
• 유사 스펙트럼 충실도: 작은 시스템에 대한 수치 시뮬레이션은 유사 고유 공간이 (지지 집합의 이동이 발생하는 교차점 근처를 제외하면) 실제 고유 공간과 높은 충실도를 유지함을 보여주었다.

의의 및 주장
본 논문은 CAS 축소가 다체 양자계의 스펙트럼 재구성에 대한 새로운 관점을 제공한다고 주장한다.

• 이는 유사 스펙트럼과 매개되는 대칭성의 국소성을 통해 상전이의 성격(1차 vs 연속)을 식별하는 진단 도구 역할을 한다.
• 근사 및 그 불확실성에 대한 완전한 분석적 설명을 제공하여, 개념적 프레임워크를 위해 시간 이산화에 의존하지 않고도(비록 대규모 수치 효율성을 위해 이산화가 사용되기는 하지만) 상평도를 연구할 수 있게 한다.
• 저자들은 이 프레임워크가 단열 양자 계산과 고전적 휴리스틱 사이의 간극을 메운다고 주장하며, CAS 재귀 트리를 통한 "경로"가 최적화 문제를 해결하기 위한 일련의 고전적 알고리즘을 정의한다고 제안한다.
• 이 방법은 텐서 네트워크와 같은 기존 기술을 대체하기보다는 보완하는 방식으로 제시되며, 로컬 분석으로는 종종 가려지기 쉬운 $A$와 $B$의 결합 특성으로부터 전역적 구조적 통찰을 추출하는 방법을 제공한다.