실시간 반복 스킴을 이용한 동적 평균장 이론과 양자 시뮬레이션 프레임워크

초록

본 논문은 실시간(retarded) 그린함수를 직접 다루는 DMFT 자기일관성 반복 알고리즘을 제안한다. 베타트 격자 반강체 모델을 대상으로, 5개의 배치 사이트를 가진 1차원 체인으로 임피던스 함수를 이산화하고, 정확대각화(ED)와 실시간 전파를 이용해 임피던스 함수를 반복적으로 피팅한다. 제한된 시간 해상도와 최소 배치 크기에도 불구하고 금속‑절연 전이와 허버드 밴드 형성을 정확히 재현한다. 이 접근법은 Trotter화 또는 VQE와 같은 근거리 양자 컴퓨팅 기법에 바로 적용 가능하도록 설계되었다.

상세 분석

이 연구는 기존 DMFT가 주로 복소수 주파수축 또는 가상시간(Imaginary‑time)에서 수행되는 점을 근본적으로 전환한다. 실시간(retarded) 그린함수를 직접 다루면 양자 하드웨어가 제공하는 유니터리 연산과 자연스럽게 매핑될 수 있다. 저자들은 베타트 격자에서 무한 차원 한계의 self‑energy가 국소적이라는 DMFT의 핵심 가정을 유지하면서, 임피던스 함수 Δ(t)를 실시간 영역에서 정의하고, Δ_R(t) = t*² G_R^loc(t) 라는 간단한 자기일관성 관계를 도출한다.

임피던스 함수를 이산화하기 위해 5개의 배치 사이트를 갖는 1차원 체인을 사용했으며, 이는 기존 두‑사이트 DMFT(두 개의 파라미터만 조정)보다 훨씬 풍부한 스펙트럼 해상도를 제공한다. 체인 파라미터 {t_i, ε_i}는 매 반복마다 Δ(t)와 Δ_R^chain(t;P) 사이의 L2 오차를 최소화하는 비선형 피팅 절차를 통해 업데이트된다. 이때 정확대각화(ED)를 이용해 배치 그린함수 g₁₁(t)를 계산하고, 실시간 전파 시 발생하는 수치적 발산을 억제하기 위해 지수 감쇠 η를 도입한다.

양자 하드웨어 구현을 염두에 둔 설계는 두 단계로 나뉜다. 첫 번째는 전체 체인(H_imp+H_hyb+H_bath)의 기저 상태를 VQE 등으로 준비하고, 두 번째는 그 상태에서 실시간 연산자를 Trotter화하거나 변분법으로 전파한다. 현재는 노이즈 없는 시뮬레이션을 위해 ED를 사용했지만, 알고리즘 자체는 양자 회로 깊이와 큐비트 수를 최소화하도록 설계되어 있다.

수렴성 측면에서 저자들은 선형 믹싱 파라미터 α를 0.8 정도로 설정해 Δ(t)의 진동을 억제하고, 10~15 회 반복 내에 자기일관성을 달성한다. 결과적으로 U/t*가 0에서 8까지 변화할 때, 금속 상태에서는 피크가 Fermi 레벨에 집중되고, 강한 상호작용에서는 피크가 사라지며 허버드 밴드가 ±U/2 근처에 나타난다. 이는 기존 두‑사이트 DMFT가 제공하지 못했던 정밀한 스펙트럼 구조를 실시간 프레임워크에서도 재현할 수 있음을 증명한다.

이 논문의 주요 기여는 (1) 실시간 그린함수 기반 DMFT 반복 스킴을 제시, (2) 최소 배치 크기로도 허버드 밴드와 Mott 전이를 포착, (3) 양자 컴퓨팅 친화적인 알고리즘 구조를 설계하여 향후 NISQ 디바이스에 바로 적용 가능하도록 한 점이다. 또한, 시간‑도메인 피팅을 통해 복소수 주파수 변환 과정에서 발생하는 수치적 불안정을 회피한다는 점도 주목할 만하다. 향후 연구에서는 배치 사이트 수를 늘려 연속적인 하이브리드화 함수를 더 정밀히 근사하고, 실제 양자 하드웨어에서 Trotter 오류와 디코히런스 효과를 포함한 실험적 검증이 필요하다.