Restoring detailed balance in non-Hermitian Markov processes

기술 요약: 비가역적(Non-Hermitian) 마르코프 과정에서의 상세 균형 복원

문제 정의
평형 상태에서 벗어난 확률 과정, 특히 이산 상태를 가진 마르코프 과정으로 모델링되는 경우, 상세 균형(detailed balance)을 깨뜨리는 비대칭 전이율을 빈번하게 포함한다. 이러한 비가역적(non-Hermitian) 확률 생성자($H$)는 엔트로피가 일시적 증가 이후 감소할 수 있는 비단조적 엔트로피 생성을 초래한다. 이러한 거동은 엔트로피 극대화나 불확정성의 지속적 증가 가설에 의존하는 표준 열역학적 해석 및 추론 기법을 저해한다. 또한, 시스템의 시간 가역성 결여는 저엔트로피 구성을 인코딩할 수 있는 평형 상태에 대한 물리적 해석을 복잡하게 만든다. 텐서 네트워크 방법론이 고차원 상호작용 시스템의 "차원의 저주"를 해결하긴 하지만, 비가역적 역학 및 그로 인한 비단조적 엔트로피라는 근본적인 문제는 해결하지 못한다.

방법론
저자들은 비가역적 확률 생성자를 실수 고유값 스펙트럼을 보존하면서 에르미트(Hermitian) 연산자로 변환하기 위해 **다이슨 맵(Dyson maps, $\eta$)**을 사용하는 것을 제안한다. 핵심 방법론은 다음과 같다:

1. 다이슨 변환(Dyson Transformation): 마스터 방정식 $\partial_t |P(t)\rangle = -H|P(t)\rangle$에 가역적 변환을 적용하여, 변환된 방정식 $\partial_t |\phi(t)\rangle = -\tilde{H}|\phi(t)\rangle$를 얻는다. 여기서 $|\phi(t)\rangle = \eta|P(t)\rangle$이고 $\tilde{H} = \eta H \eta^{-1}$이다.
2. 메트릭 연산자(Metric Operator): 이 변환은 원래 시스템의 기여도를 재가중하는 메트릭 연산자 $\Omega = \eta^\dagger \eta$를 도입한다. 이 연산자는 변환된 좌표계에서의 측정을 원래의 확률 벡터와 연결하여, 통계적 평균 $\langle O \rangle$를 $\langle \phi | \tilde{\Omega} \Xi' O' | \phi \rangle$를 통해 회복할 수 있도록 보장한다.
3. 엔트로피 복원: 변환된 시스템에서는 레니 엔트로피(Rényi entropy) $S' = -\ln |\phi(t)|^2$가 시간에 따라 단조 증가함이 증명된다($\frac{1}{2}\frac{d}{dt}S'(t) \geq 0$). 이는 원래 시스템에서 제곱 노름 $|P(t)|^2$(및 엔트로피)가 국소적 최솟값을 가질 수 있는 것과 대조적이다.
4. 계산적 구축:
   • 해석적 접근법: 저자들은 폐쇄된 연산자 대수 내에서 에르미트 생성자 $\Lambda = \ln \eta$를 찾고자 한다. $H = H^\dagger$ 조건을 교환자(commutator)를 통해 확장함으로써, $\Lambda$를 구하기 위한 일련의 대수 방정식(계수 $G_{ijk}$를 포함)을 유도한다. 이 접근법은 교환자 하에서의 폐쇄적 구조를 가정하지만, 이는 모든 비자명한 상호작용 문제에서 성립하지 않을 수도 있다.
   • 수치적 접근법: 해석적 폐쇄가 어려운 일반적인 경우를 위해 반복적 수치 알고리즘이 제안된다. 이 방법은 $\Lambda$를 대수 내의 회전축으로 취급하여, 생성자의 반-에르미트(anti-Hermitian) 성분($\Delta H$)이 허용 오차 범위 내로 줄어들 때까지 반복적인 회전을 수행한다.

주요 결과
본 방법론은 정적 접촉 그래프 상의 감염-회복(SIS) 역학 모델에 적용되었다:

• 단조적 엔트로피 성장: 변환된 시스템에서는 레니 엔트로피 $S'$가 평형에 도달할 때까지 단조 증가하는 반면, 원래 시스템의 엔트로피는 최대치를 찍은 후 감소하는 양상을 보인다.
• 통계적 등가성: 변환에도 불구하고, 메트릭 연산자가 올바르게 적용된다면 변환된 시스템에서 계산된 통계적 평균(감염된 정점의 평균 및 표준 편차)은 원래 시스템의 결과와 높은 정밀도로 일치한다(콜모고로프-스미르노프 거리 $< 10^{-3}$).
• 상전이 탐지: 두 시스템(원래 시스템과 변환된 시스템)의 레니 엔트로피 최대치는 감염률과 회복률의 비율($\beta/\gamma$)의 동일한 임계점에서 일치한다. 이는 변환된 시스템이 시간 가역성을 복원했음에도 불구하고, 원래의 비가역적 시스템에서 발견되는 상전이의 특징을 그대로 유지하고 있음을 확인시켜 준다.
• 구조 보존: 변환된 시스템의 정상 상태 분포($|\phi\rangle$)는 원래 분포($|P\rangle$)의 구조적 지표를 보존한다. 이는 다이슨 맵의 대각 성분이 변환을 지배하는 반면, 비대각 성분은 엔트로피를 증가시키기 위해 성분들을 혼합하는 역할을 함을 시사한다.

의의 및 주장
본 논문은 광범위한 마르코프 계에서 상세 균형을 복원할 수 있는 일반적이고 계산 가능한 방법을 제공한다고 주장한다. 비가역적 생성자를 에르미트 생성자로 변환함으로써, 이 접근법은 해당 시스템에 엄격하게 양수인 엔트로피 생성을 특징으로 하는 "친숙한 물리적 해석"을 부여한다. 이를 통해 엔트로피 극대화와 같은 표준 통계 물리학 도구를 사용하여 평형 상태를 결정할 수 있다.

저자들은 이 형식론이 실수 및 비음수 스펙트럼을 가진 시스템을 통해 확률 과정과 양자 역학(윅 회전 $t \to i\tau$를 통해) 사이의 연결을 확립한다고 강조한다. 이러한 유사성은 노이즈와 변동이 지배적인 어려운 고전적 확률 문제를 해결하기 위해 양자 시뮬레이터를 활용할 수 있다는 가능성을 시사한다. 그러나 저자들은 다음과 같은 한계점을 조심스럽게 언급하였다: 현재의 접근법은 실수 스펙트럼을 요구하며, 복소 고유값(진동 거동)을 가진 시스템을 직접 다룰 수 없다. 또한, 해석적 구축은 모든 상호작용 문제에서 존재하지 않을 수 있는 특정 대수적 폐쇄성에 의존한다. 이 방법은 본질적으로 비단조적 엔트로피의 복잡성을 메트릭 연서 $\tilde{\Omega}$를 통한 통계적 평균 계산의 복잡성으로 치환하는 것이다.