Tensor-Network-Based Unraveling of Non-Markovian Dynamics in Large Spin Chains via the Influence Martingale Approach

기술 요약: 영향 마틴게일(Influence Martingale) 접근법을 통한 대규모 스핀 사슬의 비마르코프적 역학의 텐서 네트워크 기반 언래블링(Unraveling)

문제 정의
개방형 양자 계(open quantum system)의 역학을 고전적으로 시뮬레이션하는 것은 시스템 크기에 따른 힐베르트 공간의 지수적 증가로 인해 상당한 난관에 봉착해 있다. 행렬 곱 상태(Matrix Product States, MPS)와 같은 텐서 네트워크 방법은 닫힌 계(closed systems)와 마르코프적 개방형 계의 스케일링 문제는 성공적으로 해결해 왔으나, 비마르코프적(non-Markovian) 역학을 시뮬레이션하는 것은 여전히 도전적인 과제이다. 기억 효과(memory effects)와 정보 역류(information backflow)를 특징으로 하는 비마르코프적 영역은 종종 시스템의 전체 이력을 추적하거나 메모리 커널을 포함하는 적분-미분 방정식을 풀어야 하며, 이는 막대한 계산 비용을 초래한다. 또한, 몬테카를로 파동 함수(Monte Carlo Wave-Function, MCWF) 방법과 같은 표준 확률적 언래블링 기술은 양자 점프(quantum jumps)를 확률적으로 해석하기 위해 양수의 붕괴율(decay rates)에 의존한다. 시간-켤레 없는(time-convolutionless, TCL) 마스터 방정식으로 기술되는 비마르코프적 영역에서는 붕괴율이 일시적으로 음수가 될 수 있으며, 이는 표준 언래블링을 무효화하는 "음의 확률" 문제를 야기한다. 기존의 비마르코프적 시스템을 위한 궤적 기반 접근법들은 소규모 시스템(예: 최대 13 큐비트)에 국한되어 왔으며, 열화(thermalization) 및 얽힘 확산(entanglement spreading)과 같은 집단 현상을 연구하는 데 필요한 다체 영역(50~100 큐비트)으로 확장하는 데 어려움을 겪고 있다.

방법론
저자들은 1차원 대규모 양자계에서 시변 마르코프적 역학과 진정한 비마르코프적 역학을 모두 처리할 수 있는 효율적인 알고리즘으로서 **텐서 점프 방법(Tensor Jump Method, TJM)**을 확장하는 방안을 제안한다. 이 접근법의 핵심은 세 가지 통합된 구성 요소로 이루어진다:

1. 텐서 네트워크 프레임워크: 이 방법은 순수 상태 궤적을 표현하기 위해 MPS를 활용하여 밀도 행렬의 지수적 스케일링을 회피한다. 시간 진화는 시간 의존 변분 원리(Time-Dependent Variational Principle, TDVP)에 의해 제어되며, 이는 정확도와 계산 비용 사이의 균형을 맞추기 위해 결합 차원(bond dimension)을 동적으로 조절한다. 알고리즘은 스즈키-트로터(Suzuki–Trotter) 분해를 사용하여 유니터리(결맞음) 진화와 소산적(비결맞음) 진화 단계를 분리한다.
2. 영향 마틴게일 형식론: 비마르코프적 TCL 마스터 방정식에서 발생하는 음의 붕괴율($\gamma_k(t) < 0$) 문제를 해결하기 위해, 저자들은 영향 마틴게일 접근법을 채택한다. 직접적인 언래블링 대신, 이들은 엄격하게 양수인 이동된 붕괴율 $r_k(t) = \gamma_k(t) + C_t$를 갖는 참조 앙상블을 도입한다 (여기서 $C_t$는 $r_k(t) > 0$을 보장하는 시간 의존적 이동값이다). 물리적 역학은 영향 마틴게일 $\mu_t$를 사용하여 참조 궤적을 재가중함으로써 복구된다. 마틴게일은 연속적인 드리프트(drift)와 이산적인 점프 기여를 포함하여 확률적으로 진화하며, 이를 통해 앙상블 평균 $\mathbb{E}_{ref}[\mu_t |\psi(t)\rangle\langle\psi(t)|]$은 정확한 비마르코프적 밀도 행렬 진화를 재현한다.
3. 영향 반경(Influence Radius): 대규모 시스템(최대 100 스핀)에서의 확장성을 확보하기 위해, 저자들은 **"영향 반경"**이라는 개념을 도입한다. 수치적 분석 결과, 비마르코프적 노이즈가 국소적 관측량에 미치는 영향은 노이즈 근원으로부터의 거리와 함께 감쇠한다는 것이 밝혀졌다. 따라서 마틴게일 보정(이동값 $C_t$ 및 그에 따른 가중치 $\mu_t$)은 시스템 전체가 아니라 관측량의 유한한 공간적 근방 내에 있는 붕 decay rate들에 대해서만 계산하면 된다. 이는 마틴게일 가중치 변동과 관련된 계산 오버헤드를 크게 줄여준다.

주요 기여

• TJM의 확장: 본 논문은 시간 의존적 붕괴율(음수가 될 수 있는 붕괴율 포함)을 수용할 수 있도록 영향 마틴게일 형식론을 통합하여 텐서 점프 방법을 일반화하였다.
• 대규모 시스템으로의 확장성: 제안된 알고리즘은 최대 100개의 스핀 큐비트를 포함하는 1차원 횡장 이징 사슬(transverse-field Ising chains)에 대해 검증되었으며, 이는 표준 방법들을 사용한 비마르코프적 궤적 시뮬레이션으로는 이전에 접근 불가능했던 규모이다.
• 영향 반경 개념: "영향 반경"의 도입은 자원 효율적인 프레임워크를 제공한다. 이는 마틴게일 보정을 위한 국소적 범위를 정량화하여, 국소적 관측량에 대한 계산 비용이 전체 시스템 크기에 따라 반드시 증가하지 않도록 하는 확장 가능한 시뮬레이션을 가능하게 한다.
• 통합 프레임워크: 이 접근법은 근본적으로 다른 알고리즘으로 전환할 필요 없이, 시간-불균질 마르코프 과정(양의 시간 의존적 비율)과 비마르코프 역학(음의 비율)을 모두 시뮬레이션할 수 있는 단일 계산 프레임워크를 제공한다.

결과

• 벤치마킹: 저자들은 30개 사이트 사슬에 대해 수치적으로 정확한 행렬 곱 연산자(Matrix Product Operator, MPO) 시뮬레이션과 궤적 기반 결과를 벤치마킹하였다. 결과는 궤적 수($N_{traj}$)에 따른 체계적인 수렴을 보여주었으며, 특정 모델 파라미터에 대해 최대 결합 차원($\chi_{max}$)에 대한 의존성이 약하게 나타났다. 이는 상대적으로 작은 궤적 앙상블($N_{traj} = 100$)에서도 방법이 올바른 역학적 경향을 포착함을 나타낸다.
• 마틴게일 거동: 시뮬레이션을 통해 영향 마틴게일의 앙상블 평균 $\mathbb{E}[\mu_t]$가 1에 가깝게 유지됨을 확인하였으며, 이는 방법의 트레이스 보존(trace preservation)을 검증한다. 마틴게일은 물리적 붕괴율 $\gamma_k(t)$가 음수인 구간 동안 특이한 시간적 변동을 보였다.
• 스케일링 및 영향 반경: 비마르코프적 붕괴율의 크기를 변화시킴으로써, 저자들은 국소적 관측량의 오차가 시스템 전체의 통합된 이동된 점프 강도(integrated shifted jump intensity)에 따라 스케일링됨을 입증하였다. 100-사이트 사슬(비마르코프 노이즈가 경계 또는 무작위 사이트에 국부화된 경우)에 "영향 반경" 개념을 적용했을 때, 마틴게일 계산을 관련 국소 근방으로 제한함으로써 $N_{traj} = 100$만으로도 높은 정확도를 유지할 수 있음을 보여주었으며, 이는 국소 측정을 위한 계산 비용을 전체 시스템 크기로부터 효과적으로 디커플링(decoupling)시킨다.

의의
본 논문은 밀도 행렬 직접 진화 및 표준 확률적 방법의 한계를 극 넘어, 비마르코프적 다체 개방계 역학을 시뮬레이션하기 위한 확장 가능하고 효율적인 경로를 제공한다고 주장한다. 텐서 네트워크를 영향 마틴게일 형식론과 결려함으로써, 저자들은 기억 효과가 중요한 영역에서 집단적 다체 현상(엔탱글먼트 확산 및 소산적 전이 등)을 연구할 수 있게 하였다. 이 방법은 막대한 계산 자원을 소비하지 않고도 현실적인 노이즈 조건 하에서 대규모 양자계의 결맞음 해제(decoherence)와 소산을 특성화할 수 있는 탐색적 연구를 위한 도구로 제시된다. 이 프레임워크는 다양한 물리적 시나리오, 미시적 배스(bath) 모델, 그리고 상호작용하는 다체 해밀토니안에 대한 향후 조사를 위한 토대로 자리매김하고 있다.