散逸構造と確率熱力学の接続——自己組織化理論はどこまで統合できるか

散逸構造理論と確率熱力学——二つの非平衡理論が交差する地点

非平衡系を記述する理論は、20世紀後半に二つの大きな流れとして発展してきた。一方は、プリゴジンとその共同研究者たちが構築した散逸構造理論であり、パターン形成・自己組織化・分岐・安定性という「構造の出現」を巨視的連続体の言葉で記述する。もう一方は、1990年代以降に確立した確率熱力学であり、揺らぎが本質的なメソスコピック・小系において、単一経路に沿った仕事・熱・エントロピー生成を確率法則として扱う。

この二つは長らく別個の理論として発展してきたが、近年の数理的研究はその橋渡しを急速に進めている。本記事では、両者の理論的基盤を整理したうえで、「散逸構造はメソスコピック確率過程の大数極限として現れる決定論的構造である」という統合的視点から、四つの有力な接続枠組みを解説する。

---

散逸構造理論の理論的基盤——局所平衡からGlansdorff–Prigogine規準まで

エントロピー生成率の双線形表式と線形応答

古典的散逸構造理論の出発点は局所平衡の仮定にある。非平衡系であっても、局所的には平衡状態と同じ熱力学変数を割り当てられると仮定することで、エントロピー生成率は流束 $J_\rho$Jρ​ と一般化力 $X_\rho$Xρ​（温度勾配・化学ポテンシャル勾配・親和力など）の積和$$P = \frac{d_i S}{dt} = \sum_\rho J_\rho X_\rho \ge 0$$P=dtdi​S​=ρ∑​Jρ​Xρ​≥0

として表現される。平衡近傍では $J_\rho = \sum_{\rho’} L_{\rho\rho’} X_{\rho’}$Jρ​=∑ρ′​Lρρ′​Xρ′​ というOnsagerの線形応答が成立し、相反関係 $L_{\rho\rho’} = L_{\rho’\rho}$Lρρ′​=Lρ′ρ​ が従う。

この線形枠組みの中で最小エントロピー生成定理が導かれるが、プリゴジン自身が明示したように、その有効性は厳密に線形で現象係数を定数とみなせる近傍平衡域にのみ限定される点は重要である。

遠平衡領域とGlansdorff–Prigogine過剰エントロピー生成

散逸構造理論の真の射程は、ここからさらに遠平衡の構造形成に踏み込んだことにある。Glansdorff–Prigogineは、非平衡定常状態からの偏差 $\delta J_\rho$δJρ​、$\delta X_\rho$δXρ​ を用いた過剰エントロピー生成の概念を導入し、$$\sum_\rho \delta J_\rho \, \delta X_\rho \ge 0$$ρ∑​δJρ​δXρ​≥0

が十分条件として満たされるとき、$\delta^2 S$δ2S がLyapunov関数として機能し安定性が従うことを示した。これが一般進化規準の古典的形式である。

プリゴジンとNicolis、プリゴジンとLefeverは、化学反応と拡散の組み合わせが均一定常状態をHopf分岐やTuring型不安定性を通じて時空間秩序へ移行させることを示した。Brusselatorはその最小模型として、化学振動・波・パターン形成の理論的試金石となった。散逸構造は孤立系から見れば「秩序の増大」に見えるが、開放系全体ではより多くの散逸を伴う自己組織化であり、エネルギー散逸と構造形成は不可分である。

---

現代確率熱力学の枠組み——経路エントロピーと普遍的揺らぎ関係

単一経路に沿うエントロピー生成とSeifertの積分揺らぎ定理

確率熱力学では、系の状態はマルコフ跳躍過程やLangevin過程として記述される。Seifertの定式化では、単一経路 $\Gamma$Γ に沿う全エントロピー生成を$$\Delta s_{\mathrm{tot}}[\Gamma] = \ln \frac{\mathbb{P}[\Gamma]}{\mathbb{P}^\dagger[\Gamma^\dagger]}$$Δstot​[Γ]=lnP†[Γ†]P[Γ]​

と前進・逆過程の経路確率比として定義し、そこから$$\bigl\langle e^{-\Delta s_{\mathrm{tot}}} \bigr\rangle = 1$$⟨e−Δstot​⟩=1

という積分揺らぎ定理が導かれる。第二法則の平均不等式 $\langle \Delta s_{\mathrm{tot}} \rangle \ge 0$⟨Δstot​⟩≥0 はこの系から自動的に回収される。

Jarzynski等式・Crooksの定理・熱力学的不確かさ関係

この経路熱力学から複数の普遍的関係が派生する。

Jarzynski等式は、初期が平衡分布であれば遠平衡駆動でも$$\bigl\langle e^{-\beta W} \bigr\rangle = e^{-\beta \Delta F}$$⟨e−βW⟩=e−βΔF

を与え、Crooksの定理は前進・逆過程の仕事分布に対して$$\frac{P_F(W)}{P_R(-W)} = e^{\beta(W – \Delta F)}$$PR​(−W)PF​(W)​=eβ(W−ΔF)

を与える。さらに**熱力学的不確かさ関係（TUR）**は、定常マルコフカレントの精度がエントロピー生成によって制約されることを示し、$$\frac{\mathrm{Var}(J_T)}{\langle J_T \rangle^2} \cdot \Sigma_T \ge 2$$⟨JT​⟩2Var(JT​)​⋅ΣT​≥2

という形で表される。Crooksの定理はJarzynski等式を含意しており、これらは互いに階層的な関係にある。

これらの関係式は、局所詳細釣り合いとマルコフ性を根本仮定とする。実験面では、光ピンセットで操作したRNAヘアピンの折り畳み・解きほぐしにおける不可逆仕事分布からCrooksの定理が検証されており、確率熱力学の中心関係式が測定可能な法則として確立している。

---

両理論の比較——何が同じで何が違うのか

基本自由度と根本仮定の違い

両理論の根本的な差異は、「何を基本対象とするか」にある。古典散逸構造理論は濃度場・温度場・流束場の平均的進化を扱い、揺らぎは主に安定性や分岐枝選択の補助変数として現れる。確率熱力学では揺らぎこそが中心であり、第二法則さえ経路確率比の形で表現される。

観点	古典散逸構造	現代確率熱力学
基本自由度	連続体場・巨視的流束	単一軌道・状態確率・確率流
根本仮定	局所平衡・構成方程式の閉鎖	マルコフ性・局所詳細釣り合い
不可逆性の表現	エントロピー生成率・一般進化規準	経路エントロピー生成・揺らぎ定理
強い適用領域	パターン形成・分岐・安定性	小系・有限時間駆動・カレント揺らぎ

過剰エントロピー生成とexcess/housekeeping分解の接点

重要なのは、古典理論の「過剰エントロピー生成」と確率熱力学のexcess/housekeeping分解が、別々に生まれたにもかかわらず明確に接続されつつある点である。Maes–Netočný（2015）は、Glansdorff–Prigogineの分解がMarkov過程ではHatano–Sasa型のアプローチとして再現されること、さらに非線形拡散では過剰エントロピー生成率が局所自由エネルギー汎関数の時間微分になることを示した。

Yoshimuraらは一般非線形力学系に対し、定常状態に依存しないexcess/housekeeping分解を提案し、化学反応ネットワーク・極限周期・カオスにも適用可能であることを示した。ここに両理論の最も明白な接点がある。

---

四つの統合枠組み——散逸構造に確率熱力学を組み込む方法論

案A｜熱力学整合的な確率化された散逸構造

古典的なBrusselatorは不可逆反応を含む決定論的最小模型として書かれるが、確率熱力学に接続するには逆反応を含むchemostattedな可逆ネットワークへの埋め込みが必要である。

系サイズ $V$V を導入して $n = Vc + \sqrt{V}\,\xi$n=Vc+V​ξ と展開すると、$V$Vの次数で巨視的反応速度方程式 $\dot{c} = F(c)$c˙=F(c) が現れ、$\sqrt{V}$V​の次数で線形雑音近似が得られる。古典散逸構造の決定論的方程式そのものがマルコフ跳躍過程の大数極限として現れるのであり、同時に局所詳細釣り合いがあれば揺らぎ定理が構造形成モデルに直接付与される。

利点は、古典理論の分岐・パターン形成と確率熱力学の経路関係式を同じ模型に載せられる点にある。ただし、もとの古典模型が効果的不可逆反応で書かれている場合、熱力学的に整合的な埋め込みが別途必要になる。

案B｜確率場としての散逸構造

空間パターンや波を扱うには、反応拡散型の確率場——SPDE（確率偏微分方程式）や場のFokker–Planck記述——が必要になる。Nardiniらは場の経路確率からエントロピー生成を定義し、空間局所・波数空間の両方で不可逆性を分解できることを示した。

Brusselatorを例にとれば、確率反応拡散系の各Fourierモード $k$k に対して線形化が可能であり、$\mathrm{Re}\,\lambda_{\max}(k_*) = 0$Reλmax​(k∗​)=0 を満たす $k_* \neq 0$k∗​=0 でTuring不安定性、$\mathrm{Re}\,\lambda_{\max}(0) = 0$Reλmax​(0)=0 でHopf不安定性が起こる。これは古典散逸構造理論そのものだが、同時に経路重みの比から場のエントロピー生成を定義できるため、**どのモードがどれだけ不可逆性を担っているか**が測定可能になる。

案C｜大偏差・準ポテンシャル・マクロ確率熱力学

数理的に最も強固な統合案は、大偏差原理を介した接続である。マルコフ跳躍過程のマクロ極限において確率密度が $p(c,t) \asymp e^{-VI(c,t)}$p(c,t)≍e−VI(c,t) の形をとり、そこから決定論的力学とマクロ揺らぎ理論が同時に導出される。

反応ネットワークのHamiltonian $H(c,\pi) = \sum_\rho r_\rho(c)(e^{\Delta_\rho \cdot \pi} – 1)$H(c,π)=∑ρ​rρ​(c)(eΔρ​⋅π−1) に対して、WKB展開からHamilton–Jacobi方程式が得られ、$\pi = 0$π=0 を代入すると古典散逸構造の決定論的方程式が回収される。定常準ポテンシャル $I_{\mathrm{ss}}(c)$Iss​(c) は条件の良い場合にマクロ力学のLyapunov汎関数になり、稀な遷移率 $\kappa_{\alpha\to\beta} \asymp e^{-V\Delta I_{\alpha\beta}}$κα→β​≍e−VΔIαβ​ によって複数の散逸構造間の枝選択やノイズ誘起遷移が定量化できる。

Maes–Netočný が示したように、近傍平衡の非線形拡散では古典的なLyapunov関数は大偏差の準ポテンシャルと接続される。これが「最も完成に近い統合の核心」と言えるゆえんである。

案D｜作用素論・幾何学的接続

第四の案は、安定性解析と揺らぎ統計を作用素のスペクトルとして統一するものである。古典散逸構造理論では均一定常状態の安定性がJacobianや分散演算子の固有値問題として記述され、確率熱力学では仕事やカレントの生成関数がtilted generator $\mathcal{L}_\lambda$Lλ​ の最大固有値で決まる。$$\psi(\lambda) = \lim_{T\to\infty} \frac{1}{T} \log \bigl\langle e^{-\lambda A_T} \bigr\rangle$$ψ(λ)=T→∞lim​T1​log⟨e−λAT​⟩

この $\psi(\lambda)$ψ(λ) はGallavotti–Cohen型対称性 $\psi(\lambda) = \psi(1-\lambda)$ψ(λ)=ψ(1−λ) を持ち、さらに決定論的TUR$$\dot{\Sigma} \cdot \frac{\mathcal{D}_{\mathrm{sc}}}{|\dot{c}|^2} \ge 2$$Σ˙⋅∣c˙∣2Dsc​​≥2

が成立する。分岐のスペクトルと揺らぎのスペクトルを、同一の作用素族として並べて扱うことができるというのがこの案の本質的な強みである。

Yoshimura–Itoは決定論的化学反応ネットワークに対してもメソスコピック揺らぎとの接続から決定論的TURを導いており、2024年には流体方程式に対するhousekeeping/excess分解と相対揺らぎ由来のEPR下界も構成されている。

---

実験的検証の現状と今後の課題

現代確率熱力学の主要関係式は小規模実験で強く検証されているが、古典的な散逸構造実験系そのものに対して、Glansdorff–Prigogine型の安定性判定・Crooks/Jarzynski型の経路関係・TUR・モード別不可逆性を同一装置・同一データ系列で統合的に検証した標準的実験例は、現時点ではまだ確立していない。これは既存研究が、小系での経路熱力学実験と、場のパターン形成に対する理論・数値研究に分かれて発展してきたためだと推測される。

今後の実証研究において重要な観測量として、エントロピー生成の分解（excess/housekeeping）、フラックス・カレント揺らぎ（TUR検証）、相関関数・応答関数・構造因子 $S(k,\omega)$S(k,ω)（揺動散逸関係の破れとして観測可能）、そして分岐枝間の遷移率と待ち時間分布（準ポテンシャル障壁 $\Delta I$ΔI による制約）が挙げられる。

空間分解した濃度場・速度場・応答関数の計測が可能な実験系では、散逸構造と確率熱力学の統合検証は十分に射程に入っている。2025年のNagayamaらによる反応拡散系の幾何学的熱力学研究は、Fisher–KPP波やBrusselatorの対称性変化で数値検証を行っており、実験への橋渡しが具体化しつつある。

---

まとめ——統合の現在地と残された問い

本記事の要点を以下に整理する。

• 散逸構造理論は局所平衡・双線形エントロピー生成・Glansdorff–Prigogine安定性を核心とし、遠平衡でのパターン形成・分岐・自己組織化を扱う
• 確率熱力学は局所詳細釣り合いとマルコフ性を基礎に、揺らぎ定理・Jarzynski等式・Crooksの定理・TURを単一経路レベルで定式化する
• 両者はexcess/housekeeping分解と大偏差準ポテンシャルを介して接続され、ここで古典的Lyapunov関数と現代的自由エネルギー汎関数・稀少経路の作用が出会う
• 四つの統合枠組みのうち、案C（大偏差・マクロ確率熱力学）を核心に、案B（確率場）と案D（作用素論・幾何学的分解）を組み合わせる構成が現時点で最も有望と判断される
• 2025年のTomé–de Oliveira論文はGlansdorff–Prigogine原理を確率熱力学から導出するという方向性を明示しており、統合は比喩ではなく具体的な数理的研究課題として進行中である

「散逸構造は、確率熱力学を土台とするメソスコピック開放系の大数極限で現れる決定論的構造であり、その安定性・枝選択・遷移・精度は、揺らぎ定理・過剰エントロピー生成・大偏差準ポテンシャル・TURによって補足される」——これが現時点での最も厳密な統合命題である。