High-temperature series expansion of the dynamic Matsubara spin correlator

기술 요약: 동적 마츠바라 스핀 상관 함수의 고온 급수 전개

문제 정의
고온 급수 전개(High-temperature series expansion, HTE)는 특히 좌절된 상호작용(frustrated interactions)을 가진 양자 스핀 시스템에서 열역학적 양과 등시각(equal-time) 스핀 상관 함수를 계산하는 표준적인 도구이다. 그러나 기존의 HTE는 정적인 등시각 관측량에 국한되어 있다. 스핀 액체와 같은 평형 양자 스핀 시스템의 집단 현상을 이해하기 위해서는 스핀 역학에 접근하는 것이 필수적이다. 동적 스핀-스핀 상관 함수 $G_{ii'}(t)$와 그 푸리에 변환인 동적 구조 인자(dynamical structure factor)는 여기 스펙트럼과 준입자의 안정성에 대한 풍부한 정보를 담고 있다. 이러한 값들은 exact diagonalization이나 텐서 네트워크와 같은 방법으로 계산될 수 있지만, 이러한 접근법들은 열역학적 극한이나 고차원 문제에서 어려움을 겪는 경우가 많다. 본 논문은 HTE를 허수 주파수에서의 동적 마츠바라 스핀-스핀 상관 함수 $G_{ii'}(i\nu_m)$로 확장함으로써 이 간극을 해결하고자 한다.

방법론: 동적 HTE (Dyn-HTE)
저자들은 무차원 결합 상수 $x = J/T$와 역 마츠바라 주파수 $1/\nu_m$의 거듭제곱에 대한 마츠바라 상관 함수의 전개를 계산하는 "동적 HTE"(Dyn-HTE)라고 명명된 방법을 개발하였다.

1. 형식론: 이 전개는 단일 결합 상수 $J$와 스핀 길이 $S \in \{1/2, 1\}$를 갖는 하이젠베르크 모델에 대해 유도되었다. 상관 함수는 다음과 같이 전개된다:
   $$T G_{ii'}(i\nu_m) = \sum_{n=0}^{n_{max}} (-x)^n c^{(n)}_{ii'}(i\nu_m)$$
   여기서 계수 $c^{(n)}_{ii'}$는 다시 $1/\nu_m$의 거듭제곱으로 전개된다. 최종 형태는 유리수 계수를 갖는 무차원 다항식 $p^{(2r)}_{ii'}(x)$를 포함한다.

2. 그래프 기반 평가: 계산은 격자 본드(bond)에 대한 합을 $n$개의 에지(edge)를 가진 위상적으로 구별되는 그래프(격자 조각) $g^{(n)}$에 대한 합으로 재구성하는 그래프 기반 접근 방식에 의존한다.

   • 임베딩(Embedding): 특정 격자에 대한 각 그래프의 기여도는 그래프 이론 알고리즘과 공간군 대칭성을 사용하여 계산된 임베딩 인자 $e(L, i, i', g^{(n)})$에 의해 결정된다.
   • 재귀적 차감(Recursive Subtractions): 연결된 기여만을 분리하기 위해(진공 및 비연결 다이어그램 제거), 저자들은 "연결된 행렬식 공식"(fermionic diagrammatic Monte Carlo에서 채택됨)에 기반한 재귀적 차감 기법을 사용한다. 이는 하위 그래프 $g^{(k)} \subset g^{(n)}$의 기여를 차감하는 과정을 포함한다.

3. 시간 적분을 위한 커널 트릭(Kernel Trick): 핵심적인 기술적 혁신은 $n$차 항에서 요구되는 $(n+2)$중 허수 시간 적분을 평가하는 것이다. 저자들은 수치 적분 대신, [참조 22]에 기반하여 닫힌 형태의 해석적 표현을 제공하는 "커널 트릭"을 활용한다.

   • $H_0 = 0$(비상호작용 스핀)인 언퍼터브드 해밀토니안의 경우, 고유 상태는 자명한 곱 상태(product states)이다.
   • 시간 순서 적분은 외부 연산자의 위치와 마츠바라 주파수 인덱스 $m$에 의존하는 "커널 함수" $\tilde{K}_{n+2}$를 사용하여 정확하게 표현된다.
   • 이를 통해 확장 계수들을 유리수로 정확하게 계산할 수 있어, 확률적 오차를 피할 수 있다.

4. 알고리즘 최적화: 구현 과정에서는 계산 비용을 줄이기 위해 대칭성(예: 순환 성질, 스핀-플레이버 대칭성)을 활용한다. 저자들은 스핀 길이 $S=1/2$ 및 $S=1$에 대해 모든 가능한 그래프에 대한 정확한 유리수 계수를 $n_{max} = 12$ (약 $1.27 \times 10^6$개의 그래프)까지 사전 계산하였다.

주요 기여

• HTE의 확장: 본 연구는 HTE 프레임워크를 정적인 등시각 상관 함수에서 동적 마츠바라 스핀-스핀 상관 함수로 성공적으로 확장하였다.
• 정확한 유리수 계수: 저자들은 모든 그래프에 대해 $J/T$의 12차 항까지 유리수 형태의 정확한 확장 계수를 제공하는 알고리즘을 제시하였다.
• 공개 저장소: 임의의 격자 임베딩을 위한 도구와 사전 계산된 계수들을 저장소에 공개하여, 임의의 사이트 쌍 또는 파동 벡터에 대한 정적인 운동량 분해 감수율(momentum-resolved susceptibilities)을 계산할 수 있도록 하였다.
• 커널 트릭의 적용: 스핀 시스템에 대한 커널 트릭의 적용은 윅의 정리(Wick's theorem)가 적용되지 않는 스핀 연산자에서 어려웠던 고차 허수 시간 적분의 해석적 해결을 가능하게 하였다.

결과 및 검증

• 벤치마킹: 본 방법은 반강자성 $S=1/2$ 하이젠베르크 사슬과 삼각 격자 모델에 대해 테스트되었다.
• QMC와의 비교: $S=1/2$ 사슬에 대해, Dyn-HTE로부터 유도된 정적 감수율을 오차 제어가 가능한 양자 몬테카를로(QMC) 시뮬레이션과 비교하였다. 결과는 우수한 일치를 보였으며, 이는 확장 계수의 타당성을 입증한다.
• 삼각 격자: 좌절된 삼각 격자 모델에 대한 정적 감수율이 보고되었으며, 이는 좌절된 시스템에 대한 방법론의 적용 가능성을 보여준다.
• 근사 테스트: 저자들은 운동량 의존성에 대한 "재규격화된 평균장 형태(renormalized mean-field form)" 파라미터화를 벤치마킹하기 위해 정적 감수율을 사용하였으며, 그 정확성을 입증하였다.

의의 및 주장
본 논문은 Dyn-HTE를 다른 방법들이 실패하거나 특징 없는 결과를 낼 수 있는 좌절된 시스템의 고온 영역에서 스핀 역학을 연구하기 위한 강력한 도구로 자리매김한다.

• 해석적 연속(Analytical Continuation): 주요한 의의 중 하나는 마츠바라 상관 함수의 고주파 확장 계수를 실주파 스펙트럼 함수의 주파수 모멘트(frequency moments)와 동일시함으로써, 해석적 연속의 난해한 문제를 우회할 수 있다는 점이다. 이를 통해 표준적인 모멘트 기반 기술(예: 계속 분수법)을 사용하여 실주파 동적 구조 인자를 재구성할 수 있다. 이는 동반 서신(companion letter)에 상세히 기술되어 있다.
• 효율성: 저자들은 커널 트릭 덕분에 동적 확장이 고차원 적분을 포함함에도 불구하고, 일반적인 등시각 상관 함수를 위한 HTE에 비해 계산 비용이 (최악의 경우에도) $n$배 정도로 완만하게 증가한다고 언급하였다.
• 범위: 이 방법은 평형 양자 스핀 시스템을 위한 일반적인 도구로 제시되나, 현재는 단일 결합 하이젠베르크 모델 및 특정 스핀 길이에 제한되어 있다. 다만, 이 형식론은 자기장이나 다양한 결합 등을 포함한 향후 확장이 가능하다.

저자들은 HTE가 기본적으로 고온 도구임을 인정하며 저온 거동에 대한 즉각적인 물리적 해석에는 신중한 태도를 유지하면서도, 다른 방법으로는 얻기 어려운 스펙트럼 함수의 정확한 제약 조건(모멘트)을 제공한다는 점을 강조한다.