상호작용 퀀치 후 자발적 공간 텍스처 형성

초록

본 연구는 실시간 실공간 Gutzwiller 근사를 이용해 반강도 반채워진 Hubbard 모델에 대한 상호작용 퀀치를 조사한다. 동질적인 초기 상태는 미세한 잡음에도 불구하고 동적으로 불안정해지며, 약한 퀀치에서는 거의 균일한 이완이, 중간 퀀치에서는 Mott‑유사 도메인의 핵생성과 느린 공동성장이, 강한 퀀치에서는 전역 코히어런스 붕괴 후 부분적인 회복이 나타난다. 이는 공간적 자기조직화가 비평형 강상관 물질의 보편적 현상임을 보여준다.

상세 분석

논문은 시간‑의존 Gutzwiller 근사(TDGA)를 실공간에 구현한 Gutzwiller‑von‑Neumann Dynamics(GvND) 방식을 도입하여, 48×48 삼각 격자상의 반채워진 Hubbard 모델을 직접 시뮬레이션한다. 초기 상태는 비상호작용( U_i=0)이며, t=0에서 온사이트 상호작용을 급격히 U_f 로 바꾸는 퀀치를 가한다. 잡음은 10⁻⁵ t_nn 수준의 Anderson 전위로 도입해 대칭을 인위적으로 깨지 않으면서도 공간적 불안정성을 촉발한다.

GvND는 슬래버‑보존 변수를 Φ_i(t)와 준입자 밀도 행렬 ρ(t)를 동시에 진화시키며, Φ_i는 로컬 전자 구성을(빈, 단일점, 이중점) 나타내는 3×3 대각 행렬이다. ρ는 유효 준입자 해밀턴 H_qp(Φ)와의 Von‑Neumann 방정식으로, Φ는 슬래버‑보존 해밀턴 H_sb(Φ,ρ)로부터 얻은 방정식으로 각각 진화한다. 이 두 시스템 간의 피드백이 공간적 비동질성 증폭의 핵심 메커니즘이다.

시뮬레이션 결과는 세 단계로 구분된다. 첫 번째는 동질적인 TDGA와 일치하는 초기 코히어런트 진동 단계로, 이때 double occupancy D(t)는 주기 T와 평균값 D*를 가진 진동을 보이며, T는 U_f/W에 대해 뚜렷한 최대를 나타낸다. 두 번째 단계에서는 D(t)의 전역 진동이 급격히 소멸하고, 공간 표준편차 σ_D(t) 가 급증한다. 이는 로컬 진동 위상이 서로 탈동조화되면서 발생하는 ‘공간 탈동조화’ 현상이다. 세 번째 단계에서는 상호작용 강도에 따라 다른 양상을 보인다.

• 약한 퀀치(U_f≪U_c)에서는 σ_D가 일시적으로 증가했지만 곧 감소하고, D(x,y)의 분포는 거의 가우시안 형태를 유지한다. 시스템은 거의 균일한 준정상 상태에 도달한다.
• 중간 퀀치(U_f≈U_c)에서는 σ_D가 크게 상승하면서 Mott‑유사 저이중점 영역이 핵생성(nucleation)되고, 시간이 지남에 따라 이들 영역이 합쳐지며 공동성장(coarsening)한다. D_i의 확률 분포는 명확한 이중 피크(bimodal)를 나타내어, 저이중점(≈0)과 높은 이중점(≈D*)이 공존함을 보여준다. 이는 공간적 상전이와 유사한 동역학적 불안정성을 의미한다.
• 강한 퀀치(U_f≫U_c)에서는 초기 진동이 빠르게 소멸하고, σ_D가 급증한 뒤 점차 감소한다. 이때 전역 코히어런스가 부분적으로 회복되며, D(x,y)의 평균값은 초기 진동 단계와 유사한 값을 유지한다. 그러나 로컬 구조는 여전히 강한 변동성을 보이며, 전하 밀도 변동 σ_n은 지속적으로 감소한다. 이는 전하 재분배보다는 상관 효과에 의해 텍스처가 형성됨을 시사한다.

이러한 현상은 기존의 DMFT이나 클러스터 DMFT가 포착하지 못하는 ‘공간적 자기조직화’를 강조한다. 특히 중간 퀀치에서 관찰된 핵생성과 공동성장은 비열적(Non‑thermal) 고정점으로의 이행을 가속화하는 새로운 이완 채널을 제공한다. 저자들은 시스템 크기에 따른 스케일링 분석을 부록 E에 제시했으며, 공동성장 시간은 L^z 형태의 크기 의존성을 가질 가능성을 언급한다.

결론적으로, 실공간 TDGA는 강상관 시스템의 비평형 동역학에서 공간적 불안정성을 자연스럽게 포착하며, 동질적인 이론이 놓칠 수 있는 중요한 물리적 메커니즘—특히 텍스처 형성, 핵생성, 공동성장—을 밝힌다. 향후 연구는 양자 플럭투에이션과 공간 플럭투에이션을 동시에 포함하는 방법론 개발이 필요함을 제시한다.