고온 전개를 이용한 동적 마츠바라 스핀 상관함수

초록

본 논문은 고온 전개(HTE)를 확장하여 마츠바라(허수 주파수) 스핀‑스핀 상관함수의 동적 버전을 계산하는 방법을 제시한다. Heisenberg 모델(J, S=½,1)에 대해 차수 12까지의 전개 계수를 그래프별로 정확히 구하고, 이를 이용해 임의의 격자와 사이트 쌍에 대한 정적·동적 전단위자극성을 효율적으로 얻는다. 또한 고주파 전개 계수를 실주파수 스펙트럼의 모멘트와 연결시켜, 전통적인 해석적 연속 없이 동적 구조인자를 재구성할 수 있음을 보인다.

상세 분석

본 연구는 기존 고온 전개(HTE)가 열역학량이나 동시 시간 상관함수에 국한됐던 한계를 넘어, 허수 시간(마츠바라) 영역에서의 스핀‑스핀 상관함수 G(iνₘ)까지 확장한다는 점에서 혁신적이다. 핵심 아이디어는 Heisenberg Hamiltonian을 비상호작용 부분 H₀=0와 상호작용 부분 V=J∑V_{ij} 로 분리하고, 표준 섭동 이론을 이용해 G(iνₘ)의 전개 계수를 c^{(n)}(iνₘ) 형태로 전개한다. 여기서 n차 전개는 n개의 결합(bond) 연산자 V_{b}와 두 스핀 연산자 Sᶻ_i(τ), Sᶻ_{i’}(τ’)의 시간 순서 적분을 포함한다. 허수 시간 적분은 (n+2) 차원 다중 적분이 되며, 이를 해결하기 위해 최근 개발된 “Kernel trick”을 적용해 모든 적분을 폐쇄 형태의 유리수로 변환한다.

그래프 이론을 도입해 결합들의 조합을 ‘그래프’(order‑n graph)로 묶고, 각 그래프에 대한 기여 c_g^{(n)}(iνₘ)를 미리 계산한다. 전체 격자에 대한 전개는 그래프의 삽입 수(e(L,i,i’,g))를 가중치로 하여 합산한다. 이 과정은 기존 HTE에서 사용되는 ‘링크드‑클러스터’ 기법과 동일하지만, 허수 시간 적분이 추가된 점이 차별점이다. 전개 차수가 12까지 확장된 이유는 그래프 수가 약 10⁶개에 달하지만, 각 그래프에 대한 연산 복잡도는 n에 비례하는 수준으로 제한돼 실용적인 계산이 가능했다.

전개 결과는 두 부분으로 나뉜다. 첫 번째는 x=J/T에 대한 다항식 p^{(0)}(x) 로, 이는 정적(νₘ=0)인 스핀 감수성 χ_{ii’}와 동일하다. 두 번째는 역마츠바라 주파수 1/νₘ의 짝수 차수 항 p^{(2r)}(x) 로, 이는 고주파 전개 계수이며 실시간 상관함수의 짧은 시간 전개와 직접 연결된다. 특히 p^{(2r)}(x)∝⟨