2차원 허버드 모델의 스핀·전하 전도성, 실시간 클러스터 확장으로 밝히다

초록

본 논문은 수치 연결 클러스터 전개(NLCE)와 그 동적 확장(dNLCE)을 이용해 2차원 정사각 격자 퍼미-허버드 모델의 스핀 및 전하 광학 전도성을 실시간 전류‑전류 상관함수로부터 직접 계산한다. 제한된 시간 구간을 보완하기 위해 피팅과 드루드 가중치(D = Drude weight)를 도입하고, 1차원 DMRG 결과와 비교해 정확성을 검증하였다. 강한 상호작용·반반점 채움에서 스핀 전도도의 DC값이 실험과 일치하도록 드루드 항이 핵심 역할을 함을 확인하였다.

상세 분석

이 연구는 2차원 퍼미‑허버드 모델(FHM)의 동적 응답을 정확히 구하기 어려운 문제를 해결하고자, 기존의 수치 연결 클러스터 전개(NLCE)를 시간 의존 형태로 확장한 dNLCE 방식을 도입하였다. NLCE는 무한 시스템의 광범위한 온도·화학 퍼텐셜 그리드에서 클러스터별 정확한 해를 이용해 열역학적 극한을 직접 계산할 수 있는 장점이 있다. 여기서는 전류 연산자를 포함한 4점 상관함수를 계산해야 하므로 클러스터 차수를 8차까지 확장했으며, 1차원에서는 10차까지 수행했다.

시간 영역에서 상관함수는 짧은 시간(τ≈1–2 /t)까지 수렴하지만, 긴 시간으로 갈수록 급격히 발산한다. 이를 보완하기 위해 저자들은
(y(\tau)=2D/β + A e^{-B\tau} + \sum_{n=1}^{N} C_n \cos(E_n\tau - F_n))
형태의 피팅 함수를 제안하고, 드루드 가중치 D를 자유 파라미터로 포함시켰다. 피팅 파라미터는 실제 데이터와의 최소 제곱 차이뿐 아니라, 세 가지 Kubo‑식(식 4‑6)으로부터 얻은 전도도 스펙트럼이 서로 일치하도록 추가 제약을 두었다. 이렇게 얻은 피팅 함수는 실시간 상관함수의 장기 행동을 효과적으로 추정하고, 푸리에 변환을 통해 실주파수 영역의 전도도 σ(ω)를 얻는다.

전도도는
(\sigma(\omega)=2\pi D\delta(\omega)+\sigma_{\text{reg}}(\omega))
형태로 분리되며, Drude 항은 DC 전도도의 핵심 기여를 의미한다. 특히 스핀 전도도의 경우, 낮은 온도·강한 상호작용 영역에서 D가 비제로임을 확인했으며, 이는 기존 QMC·DMFT 계산이 놓친 부분을 메워준다. 전하 전도도는 반반점에서 금속성 행동을 보이지만, 강한 결합(U/t≈8)에서는 Drude 가중치가 거의 사라져 절연성에 가까운 스펙트럼을 나타낸다.

계산 결과는 1차원 DMRG와 정량적으로 일치함을 보였으며, 2차원에서는 최근 광학 격자 실험 데이터와 비교해 전하 DC 전도도가 실험값과 잘 맞는다. 특히 스핀 전도도의 경우, 실험에서 관측된 이론‑실험 격차가 Drude 항을 포함함으로써 해소됨을 보여준다. 이는 강한 결합·반반점에서 스핀 전하가 각각 다른 메커니즘으로 전도성을 유지한다는 물리적 통찰을 제공한다.

마지막으로, 저자들은 dNLCE가 시간‑주파수 변환 과정에서 발생하는 불확실성을 피팅과 합성 규칙을 통해 정량화했으며, 향후 더 큰 클러스터 차수와 고성능 컴퓨팅을 결합하면 저온·저에너지 스케일까지 접근 가능할 것으로 전망한다.