긴거리 경쟁 상호작용을 갖는 개방 양자 시스템 변분 접근법

초록

본 논문은 행렬곱 연산자(MPO)와 시간‑의존 변분 몬테카를로(t‑VMC)를 결합한 새로운 알고리즘을 제시한다. 이 방법은 1·2차원 스핀‑½ 격자에서 알제브라적으로 감소하는 장거리 경쟁 상호작용을 포함한 개방 양자 시스템의 비평형 동역학과 정상 상태를 효율적으로 시뮬레이션한다. N = 200까지 확장 가능한 이 접근법은 트로터 오류 없이 Lindblad 마스터 방정식을 변분적으로 풀어, 공간적으로 변조된 자기질서가 비평형 정상 상태에서 나타나는 현상을 밝혀낸다.

상세 분석

이 연구는 개방 양자 다체 문제를 해결하기 위해 두 가지 강력한 텐서 네트워크 기법을 융합한다. 먼저, 시스템의 밀도 행렬을 행렬곱 연산자(MPO) 형태로 표현함으로써, 고차원 리우빌리언 공간을 효율적으로 압축한다. MPO는 각 사이트의 물리 차원을 d와 결합 차원 χ(결합 차원)으로 구성된 텐서들의 사슬이며, 이를 벡터화하면 행렬곱 상태(MPS)와 동일한 구조가 된다.

시간 전개는 Lindblad 마스터 방정식 dρ/dt = ℒρ 로 기술되며, 변분 원칙(Dirac‑Frenkel) 을 적용해 변분 파라미터 a = {a_i}가 정의하는 매니폴드 M 위에서 정확한 동역학 ℒρ 를 정규 직교 투사한다. 이 과정에서 얻어지는 변분 방정식은 S·ȧ = f 형태이며, 여기서 S는 파라미터 공간의 메트릭 텐서(그람 행렬), f는 “변분 힘”이라 부르는 ℒρ와의 내적이다.

S와 f는 직접 계산하면 지수적으로 커지는 힐베르트 공간 때문에 불가능하므로, 저자들은 Monte Carlo 샘플링을 도입한다. 구체적으로, 현재 MPO가 정의하는 확률분포 p(x) ∝ |⟨x|ρ⟩|² 에서 샘플 {x}를 추출하고, 각 샘플에 대해 로그미분 연산자 Δ_i와 로컬 Lindblad 추정량 L_loc(x) = ⟨x|ℒ|ρ⟩/⟨x|ρ⟩ 를 계산한다. 이렇게 하면 S_ij = ⟨Δ_i†Δ_j⟩_p,  f_i = ⟨Δ_i† L_loc⟩_p 로 표현되어, 샘플 평균만으로 근사값을 얻을 수 있다.

특히, 기존 신경망 기반 VMC와 달리 여기서는 Δ_i와 L_loc을 직접 텐서 수축으로 평가한다. 이는 비국소적인 장거리 상호작용을 포함한 Lindblad 연산자를 n‑local 연산자들의 합으로 분해하고, 각 항을 MPO 텐서와 효율적으로 계약함으로써 구현된다. 결과적으로, 복잡한 1/r^α 형태의 상호작용도 높은 정확도로 시뮬레이션 가능하다.

시간 전진은 Heun 2차 방법과 적응형 스텝 사이징을 사용해 안정성을 확보하고, 신호‑대‑노이즈 비율을 정규화해 수치적 발산을 방지한다. 중요한 점은 이 변분 스킴이 Trotter 분해를 전혀 사용하지 않으므로, 장시간 전개에서 누적되는 Trotter 오차가 없다는 것이다. 따라서 정상 상태에 대한 정밀한 추정이 가능해진다.

벤치마크에서는 1차원 이방성 Heisenberg 체인(N = 200, χ = 20)과 2차원 사각 격자(N = 4×4)에서 t‑VMC+MPO 결과를 전통적인 t‑MPS(두 번째 차수 Suzuki‑Trotter)와 비교하였다. 자기화, 근접·다음 근접 스핀‑스핀 상관함수, 그리고 정상 상태 자기화까지 전반적으로 뛰어난 일치를 보였으며, 특히 정상 상태 수렴 속도와 정확도에서 t‑VMC+MPO가 우수함을 확인했다.

마지막으로, 장거리 경쟁 상호작용(예: J₁ < 0, α₁ = 3인 반데르발스형, J₂ > 0, α₂ = 6인 반강제형)을 포함한 1·2차원 모델을 조사하였다. 비평형 정상 상태에서 스핀 밀도 파동이 형성되어, 공간적으로 변조된 자기질서(스트라이프드 패턴)가 나타났다. 이는 기존 폐쇄계에서 보고된 프러스트레이션 현상이 개방계에서도 유지될 수 있음을 시사한다.

전반적으로, 이 논문은 MPO 기반 변분 Monte Carlo가 장거리 상호작용과 외부 구동·소산을 동시에 다루는 개방 양자 다체 문제에 대해 확장성, 정확성, 구현 용이성 측면에서 새로운 표준이 될 가능성을 보여준다.