Restoring detailed balance in non-Hermitian Markov processes

本論文は、ディソン写像を用いることで、詳細釣合いが破れた非エルミート・マルコフ過程においてエントロピーの単調増加を回復させ、それによって広範な非平衡系への標準的な統計物理学のツールや推論手法の適用を可能にする、一般的かつ計算量的に扱いやすい手法を導入するものである。

要約

あなたは、ある人々が繰り広げる混沌とした「椅子取りゲーム」を見ていると想像してください。何人かは立っており（感染）、何人かは座っています（健康）。彼らは特定のルールに基づいて、動き回り、立ち上がり、座ります。科学の世界では、これは**確率過程（stochastic process）**と呼ばれます。つまり、偶然と不確実性に支配されたシステムです。

通常、科学者はこうしたシステムをエントロピーを見ることで理解しようとします。エントロピーとは、「乱雑さ」や「不確実性」の尺度だと考えてください。完璧で穏やかな世界では、乱雑さは時間の経過とともに、同じままか、あるいはより乱雑になる（綺麗になることはない）ことが期待されます。これが「時間の矢」です。

しかし、この論文は、多くの現実世界のシステム（病気の蔓延や噂の広まりなど）において、この「乱雑さ」は素直に振る舞わないことを説明しています。それは上がったり、下がったり、また上がったりします。これは非単調な動きです。このため、科学者が次に何が起こるかを予測するために標準的な数学的ツールを使うことは非常に困難になります。なぜなら、それらのツールは、物事が落ち着くまで常に乱雑になり続ける（あるいは同じままである）という概念に基づいているからです。

問題点：壊れた鏡

著者らは、このシステムを行列（ここではHと呼びます）という数学的対象を用いて記述しています。「扱いやすい」システムでは、この行列は対称であり、完璧な鏡のようです。鏡に映った像は、元の物体と一致します。しかし、こうした混沌とした不均衡なシステムでは、行列は**非エルミート（non-Hermitian）**になります。

自分の姿が歪んで見える不思議な鏡を見ているところを想像してみてください。反射した像は存在しますが、それは歪んでいます。この歪みのために、システムの「乱雑さ（エントロピー）」は奇妙な挙動を示し、本来増えるべき時に減少することさえあります。これは、科学者が世界を理解するために用いる標準的な物理法則を打ち砕くものです。

解決策：「ダイソン写像（Dyson Map）」

著者らは、ダイソン写像と呼ばれる巧妙な数学的トリックを紹介しています。これは、特別なメガネやレンズのようなものだと考えてください。

この歪んだ「不思議な鏡」のシステムに、この「メガネ」をかけると、魔法のようなことが起こります：

1. 歪みが消える： 歪んだ行列が、完璧な対称行列（エルミート行列）へと変換されます。
2. ルールが回復する： 突然、「乱雑さ（エントロピー）」が正常に振る舞い始めます。それは、システムが穏やかで安定した状態に達するまで、着実に、かつ滑らかに増加していきます。
3. 視点は正確なまま： たとえメガネを通した見え方が異なっていても、基礎となる情報は保持されています。以前と同じ確率や結果を計算することができます。ただし、最終的な答えに小さな「補正係数（メトリック）」を適用する必要があります。

比喩：ガタつくテーブルを直すこと

左右にぐらぐらと傾き続ける、ガタついたテーブルを想像してください。動きが予測できないため、その上にコーヒーカップを置くのは困難です。

ダイソン写像は、調整可能な脚を追加することに似ています。テーブル自体を動かすのではなく、それを見る「視点」を変えるのです。すると、突然、テーブルは完全に水平に見えます。コーヒーカップは静止しています。「ぐらつき」は消えました。

しかし、テーブルは依然として同じ床の上にあります。もし部屋のどこにテーブルがあるのかを正確に知りたい場合は、脚を追加したことを覚えておかなければなりません。この論文は、その「水平なテーブル」の視点を、元の「ガタついたテーブル」の現実へと翻訳するための数学を提供しています。これにより、情報が失われることなく、元の現実へと戻すことができます。

検証内容：疾患モデル

この手法が機能することを証明するために、著者らは病気の蔓延モデル（SISモデルと呼ばれる）にこの方法を適用しました。

• 修正前： 彼らは病気の広がりにおける「不確実性」を観察しました。病気が広がるにつれて不確実性は上昇し、その後落ち着くにつれて「低下」し、再び上昇しました。それは混乱した、ギザギザの線でした。
• 修正後： ダイソン写像を用いることで、彼らは同じプロセスを観察しました。不確実性は、まるで丘を転がり落ちるボールのように、滑らかに、かつ着実に上昇し、底に到達しました。

この滑らかで予測可能な挙動により、病気が消滅するか流行（エピデミック）になるかの「分岐点（相転移）」を容易に見つけることができました。元の、壊れた視点では、これを特定することは非常に困難でした。

大きな展望

この論文は、この「レンズ（ダイソン写像）」を用いることで、標準的な熱力学の法則を破るような、乱雑で予測不可能なシステムを、ルールに従うクリーンで秩序あるシステムへと変容させることができると主張しています。

これは、システムの物理的な現実を変えるものではありません。単に、それを記述するために使用される数学的な言語を変えるだけです。これにより、研究者は、これまであまりに混沌としていて扱えなかった問題（例えば、病気がどのように広がるか、あるいは噂がどのように群衆の中を伝わるかなど）を、まず数学を「真っ直ぐに直す」ことで、強力で標準的なツールを用いて解決できるようになるのです。

技術概要

技術要約：非エルミート・マルコフ過程における詳細釣合の回復

問題提起
非平衡状態にある確率過程、特に離散状態を持つマルコフ過程によってモデル化されるプロセスでは、詳細釣合を破る非対称な遷移率が頻繁に含まれる。これらの非エルミート的な確率生成子（$H$）は、エントロピーが一時的な増加の後に減少するという、非単調なエントロピー生成を引き起こす。この挙動は、エントロピー最大化や、不確実性の絶え間ない増大という仮定に依存する標準的な熱力学的解釈や推論手法を根底から揺るがすものである。さらに、これらのシステムにおける時間可逆性の欠如は、低エントロピー構成を符号化している可能性のある平衡状態の物理的解釈を困難にする。テンソルネットワーク法は、高次元の相互作用系における「次元の呪い」に対処するが、非エルミート動力学とそれに伴う非単調なエントロピーという根本的な問題は解決しない。

手法
著者らは、非エルミート的な確率生成子を、実固有値スペクトルを保持したままエルミート演算子へと変換するために、ダイソン写像（$\eta$）の使用を提案している。核心となる手法は以下の通りである：

1. ダイソン変換： マスター方程式 $\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. 計量演算子： この変換により、元のシステムへの寄与を再重み付けする計量演算子 $\Omega = \eta^\dagger \eta$ が導入される。この演算子は、変換された座標系における測定を元の確率ベクトルへと結びつけ、統計的平均 $\langle O \rangle$ が $\langle \phi | \tilde{\Omega} \Xi' O' | \phi \rangle$ を通じて回収されることを保証する。
3. エントロピーの回復： 変換されたシステムにおいて、レニー・エントロピー $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$ を交換子を用いて展開することにより、$\Lambda$ を解くための一連の代数方程式（係数 $G_{ijk}$ を含む）を導出する。このアプローチは、交換子の下での閉じた構造を仮定しているが、これはすべての非自明な相互作用問題において成立するとは限らない。
   • 数値的手法： 解析的な閉鎖が困難な一般的なケースに対しては、反復的な数値アルゴリズムが提案されている。この手法は $\Lambda$ を代数内の回転軸として扱い、生成子の反エルミート成分（$\Delta H$）を許容誤差内に抑えるまで、反復的な回転を行う。

主な結果
本手法は、静的な接触グラフ上の感受性・感染・感受性（SIS）流行病モデルに適用された：

• 単調なエントロピー増大： 変換されたシステムでは、平衡に達するまでレニー・エントロピー $S'$ が単調に増加するが、元のシステムのエントロピーは一時的な極大値を示した後に減少する。
• 統計的等価性： 変換にもかかわらず、変換されたシステムで計算された統計的平均（感染ノードの平均および標準偏差）は、計量演算子 $\tilde{\Omega}$ が正しく適用されている限り、元のシステムの計算結果と高い精度（コルモゴロフ・スミルノフ距離 $< 10^{-3}$）で一致した。
• 相転移の検出： 両方のシステム（元のシステムと変換されたシステム）におけるレニー・エントロピーの極大値は、感染率と回復率の比率（$\beta/\gamma$）において一致している。これは、変換されたシステムが、時間可逆性を回復しているにもかかわらず、元の非エルミート系に見られる相転移のシグネチャを保持していることを裏付けている。
• 構造の保存： 変換されたシステムにおける定常分布（$|\phi\rangle$）は、元の分布（$|P\rangle$）の構造的指数を保持している。これは、ダイソン写像の対角成分が変換を支配し、非対角成分がエントロピーを増大させるための混合として機能することを示唆している。

意義と主張
本論文は、広範なマルコフ系において詳細釣合を回復させるための、一般的かつ計算可能な手法を提供すると主張している。非エルミート的な生成子をエルミート的なものに変換することで、このアプローチは、厳密に正のエントロピー生成によって特徴付けられる「馴染み深い物理的解釈」をこれらのシステムに付与する。これにより、平衡状態を決定するためのエントロピー最大化などの標準的な統計物理学のツールを適用することが可能になる。

著者らは、この形式体系が、実数かつ非負のスペクトルを持つシステムを通じて、確率過程と量子動力学（ウィック回転 $t \to i\tau$ による）との間の接続を確立していることを強調している。この類似性は、量子シミュレータが、ノイズやゆらぎに支配される困難な古典的確率問題を解くために利用できる可能性を示唆している。しかし、著者らは以下のような限界についても謙虚に述べている。現在の手法は実数スペクトルを必要とし、複素固有値（振動挙動）を持つシステムには直接対処できないこと、また、解析的な構築は、すべての相互作用問題に存在するとは限らない特定の代数的閉鎖に依存していることである。本手法は、非単調なエントロピーの複雑さを、計量演算子 $\tilde{\Omega}$ を介した統計的平均の計算の複雑さと実質的にトレードオフしている。