Worldline-Susceptibility Scheduling for Quantum Annealing Beyond Local-Adiabatic Evolution

기술 요약: 양자 어닐링을 위한 월드라인 자화율 스케줄링

문제 정의
양자 어닐링(Quantum Annealing, QA)의 성능은 전체 실행 시간을 진화 경로를 따라 어떻게 배분할지를 결정하는 어닐링 스케줄에 결정적으로 의존한다. 이론적으로 최적인 롤랜드-서프(Roland–Cerf) 국소 단열 스케줄은 순간적 스펙트럼 갭($\Delta(s)^{-2}$)에 반비례하도록 시간을 할당하지만, 대규모 최적화 문제에서는 실질적으로 불가능하다. 이러한 스케줄을 구축하려면 즉각적인 저에너지 스펙트럼에 대한 완전한 지식이 필요하며, 이는 지수적으로 거대한 해밀토니안의 반복적인 대각화를 요구한다. 결과적으로, 명시적인 스펙트럼 정보 없이도 임계 영역을 식별할 수 있는 계산 비용이 저렴한 대리 모델(surrogate)이 필요하다.

방법론
저자들은 시뮬레이션 양자 어닐링(Simulated Quantum Annealing, SQA) 과정에서 측정된 월드라인 자화율(worldline magnetization susceptibility, $\chi_m$)에 기반한 대리 스케줄링 프레임워크를 제안한다. 방법론은 다음과 같이 진행된다:

1. SQA 시뮬레이션: 양자 어닐링 과정을 스즈키-트로터(Suzuki–Trotter) 분해를 사용하여 동등한 고전 통계 역학 시스템으로 매핑한다. 이는 횡장 이징 모델(transverse-field Ising model)을 $M$개의 결합된 복제본(Trotter slices)을 가진 $(d+1)$ 차원 고전 시스템으로 변환한다.
2. 관측량 측정: 평형 상태의 구성(configurations)을 고전 몬테카를로 기법을 사용하여 샘플링한다. 이를 통해 월드라인 자화 $m$을 계산하고, 요동-소산 정리(fluctuation–dissipation theorem)를 통해 자화율 $\chi_m$을 도출한다: $\chi_m = NM(\langle m^2 \rangle - \langle |m| \rangle^2)$.
3. 대리 모델 구축: 측정된 자화율 $\chi_m(s)$는 역제곱 스펙트럼 갭의 대리 지표 역할을 한다. 시간 할당 가중치 $w(s) = \chi_m(s) + \chi_0$ (여기서 $\chi_0$는 정규화 상수)를 정의한다. 누적 시간 할당 함수 $\tau(s)$는 이 가중치를 적분하여 구축되며, 최종 스케줄 $s(t)$는 $\tau(s)$를 역전(inverting)하여 얻는다.
4. 검증: 오픈 소스 소프트웨어인 Qanneal로 구현된 이 프레임워크는 셔링턴-커크파트(Sherrington–Kirkler) 스핀 글래스 인스턴스에 대해 정확한 대각화(Exact Diagonalization, ED) 결과와 대리 스케줄을 검증한다. 성능은 선형, 정확한 롤랜드-서프, 그리고 제안된 대리 스케줄을 사용하여 시간 의존 슈뢰딩거 방정식을 풀었을 때 얻은 최종 바닥 상태 확률($P_{GS}$)을 통해 정량화된다.

주요 결과
본 연구는 월드라인 자화율 스케줄이 상당수의 SK 인스턴스에 대해 기존의 선형 어닐링 및, 역설적이게도 정확한 롤랜드-서프 스케줄보다 일관되게 우수한 성능을 보임을 입증한다.

• 임계 영역 식별: 월드라인 자화율의 피크는 상세한 갭 프로파일을 정량적으로 재현하지는 못하더라도, 스펙트럼 갭이 최소가 되는 광범위한 영역을 신뢰성 있게 식별한다.
• 우수한 유한 시간 성능: 대표적인 인스턴스($n=10$ 및 $n=12$)에 대해, 대리 스케줄은 다양한 총 어닐링 시간 동안 정확한 롤랜드-서프 스케줄보다 높은 바닥 상태 확률을 생성한다.
• 실패 모드 발견: 저자들은 정확한 국소 단열 스케줄의 성능 저하를 설명하는 두 가지 뚜렷한 유한 시간 실패 모드를 식별하였다:
  1. 경계-갭 트랩(Boundary-Gap Trap): 최소 갭이 어닐링의 끝($s \approx 1$)에서 발생할 때, 롤랜드-서프 스케줄은 횡장 구동기(transverse-field driver)가 사실상 사라진 영역에 대부분의 실행 시간을 할당한다. 이는 추가된 시간이 동적으로 효과를 발휘하지 못하게 만든다. 자화율 기반 스케줄은 양자 요동이 여전히 활발한 더 넓은 크로스오버 영역에 시간을 집중함으로써 이를 방지한다.
  2. 진동 불안정성(Oscillatory Instability): 최소 갭이 내부(interior)에 위치할 때, 정확한 스케줄은 갭 주변에 실행 시간을 너무 급격하게 집중시켜 코히어런트 다준위 간섭(coherent multilevel interference)을 유발한다. 이는 성공 확률이 총 어닐링 시간에 따라 비단조적이고 진동하는 의존성을 갖게 만드는데, 더 매끄러운 자화율 기반 스케줄은 이를 회피한다.
• 통계적 견고성: 무작위성 실현(disorder realizations)과 시스템 크기($n=10$ ~ $n=20$)에 걸친 대규모 수치 실험은 대리 스케줄의 우수성이 특정 인스턴스의 산물이 아닌 일반적인 특징임을 확인시켜 준다. 시스템 크기가 커짐에 따라 경계-갭 인스턴스의 빈도가 높아지며, 이는 정확한 롤랜드-서프 스케줄의 성능을 대리 스케줄에 비해 더욱 악화시킨다.

의의 및 주장
본 논문은 현실적인 유한 시간 조건 하에서는 스케줄링에 있어 **견고성(robustness)**이 즉각적인 스펙트럼 갭 프로파일을 정밀하게 재현하는 것보다 더 중요하다고 주장한다. 저자들의 주장은 다음과 같다:

• 저렴한 평형 관측량(구체적으로 SQA로부터의 월드라인 자화율)이 스펙트럼 갭 기반 전략의 효과적인 대리 지표 역할을 할 수 있다.
• 정확한 국소 단열 조건은 점근적 극한에서는 이론적으로 최적이지만, "경계-갭 트랩"과 "진동 불안정성"으로 인해 유한 시간 영역에서는 오히려 역효과를 낼 수 있다.
• 제안된 방법론은 정확한 스펙트럼 계산의 계산적 불가능성을 우회하여 스케줄 최적화를 위한 확장 가능한 경로를 제공하며, 이를 통해 정확한 대각화가 도달할 수 없는 규모의 문제에도 적용 가능하다.

이 작업은 Qanneal 프레임워크에 구현되어 SQA, 자화율 측정 및 스케줄 구축을 위한 통합 플랫폼을 제공함으로써, 양자 어닐링 역학에 대한 재현 가능한 연구를 촉진한다.