Gray-Scott 散逸構造

反応拡散系シミュレーターの使い方と、背景にある非平衡熱力学・散逸構造理論の解説。

使い方

Gray-Scott方程式とは

∂u/∂t = Du∇²u − uv² + F(1−u)
∂v/∂t = Dv∇²v + uv² − (F+k)v

u（基質）とv（触媒）という2つの化学物質の濃度場が、拡散しながら反応する系です。uv²という自己触媒的な反応項（vがvを触媒的に増やす）が、Turingパターンと同じ空間パターン形成の駆動力になります。

F(1−u)は外部から系にuが流入する速さ（開放系のソース項）、(F+k)vは系からvが排出される速さ（シンク項）です。このFとkがGray-Scott系を非平衡開放系たらしめている核心で、シミュレーターの「系を閉じる」トグルはこのFを0にするだけの操作です。

なぜ「散逸構造」なのか

1977年ノーベル化学賞のIlya Prigogineは、「熱力学的平衡から十分離れた開放系では、非線形な相互作用がエントロピーを環境に捨てながら自己組織化した構造（散逸構造）を安定に維持できる」ことを理論化しました。Bénard対流セルや、後述するベロウソフ・ジャボチンスキー反応が代表例です。

Gray-Scottのパターンも同じ枠組みで理解できます。開放系(F>0)のときだけパターンが維持され、系を閉じる(F=0)と外部からの物質流入が止まり、パターンは徐々に崩れて熱力学的平衡（一様な状態）に緩和していきます。これはシミュレーター上で実際にトグルを切り替えて確認できます。

ベロウソフ・ジャボチンスキー反応との関係

実は、Prigogineが散逸構造理論を作る直接のきっかけになったのが、実際の化学反応であるベロウソフ・ジャボチンスキー反応 (BZ反応)です。シャーレに薄く広げると、ターゲットパターンや回転する渦巻き（スパイラル波）が実際に色の変化として現れます。数学的には「Oregonator」というBZ反応専用の反応拡散モデルがあり、Gray-Scottと同じ「反応項＋拡散項」の枠組みに属します。

そして普通のフラスコ（バッチ反応器）でBZ反応をやると、試薬が消費され尽くして振動が止まります——これはまさに「系を閉じる」トグルと同じ現象です。振動・パターンを永続的に維持するには、試薬を絶えず供給し続けるCSTR（連続撹拌槽型反応器）という装置が必要で、これが化学の世界におけるFの実体です。

Pearsonの分類とプリセット

Gray-ScottのF・kパラメータ空間は、Pearson (1993) がギリシャ文字で約12クラスに分類して以来、30年以上にわたって世界中で調べ尽くされています（MROBのアーカイブに数千点の座標が公開）。このシミュレーターのプリセットはその代表的な座標です。

プリセットFk特徴
Mitosis0.03670.0649自己複製スポット（ただし種の形状に敏感）
Coral0.05450.0620珊瑚状の分岐パターン
Worms0.05800.0650虫状・迷路状パターン
Solitons0.03000.0620孤立した安定パターン
Spots0.02500.0600安定した斑点
Waves0.01400.0410波状・空間全体を占拠
Holes0.03900.0580穴状パターン
Chaos0.02600.0510渦巻き・脈動する時空カオス

同じ(F,k)の値でも、数値計算の離散化（格子解像度・ラプラシアンの重み・時間刻み幅）によって挙動が変わることがあります。これは分野で広く知られた性質で、実装のバグではありません。

HとD：エントロピーと散逸速度の指標

このシミュレーターは2種類の統計量をリアルタイムで計算しています。この2つを混同しないことが重要です。

H：パターンのShannonエントロピー

v場の値を32ビンのヒストグラムに分け、H = −Σp·log(p)で計算し、log(32)で正規化した空間的な「乱雑さ」の指標です。構造（斑点・渦巻き）が形成されるとHはむしろ下がります——これはSchrödingerの「負のエントロピーを食べて生きる」という話と同じ構図で、局所的な秩序（低エントロピー）は環境側により大きなエントロピー増加を押し付けることで成立します。

D：散逸(流束)速度の近似

D ≈ 平均[F(1−u) + kv]。系を流れる物質流束の大きさで、非平衡定常状態でのエントロピー生成速度の近似指標です。「系を閉じる」をONにするとF=0になりDは急落します。

実験的には、H（パターンの複雑さ）とD（散逸速度）はプリセット間でおおむね同じ順位で動く傾向が見られました（Waves > Chaos > Holes > Spots > Coral > Solitons≈Worms > Mitosis）。ただしこの2つは根っこでは「ドメイン中でどれだけの割合が反応中か」を反映しているという点で独立ではなく、解釈には注意が必要です。

Maximum Entropy Production Principle (MEPP) との関連

「非平衡から遠く離れた開放系では、可能な散逸構造のモードのうち、最もエントロピー生成速度が大きいものが安定化・選択される」という仮説（MEPP）が、非平衡熱力学で今も議論されています。Jeremy Englandの"dissipative adaptation"（散逸適応）理論は、これをより厳密な確率論的熱力学から導いた現代的な発展形です。ただしMEPPは証明された法則ではなく、反例や適用限界が指摘されている、あくまで仮説であることには注意が必要です。

膜・区画化という視点

Gray-ScottやLeniaの濃度場は連続的で、境界も同じ物質の濃淡にすぎません。これに対し、生命の本質を「膜で自己を外界から区切り、その膜自体を内部の代謝ネットワークが作り続けるシステム」とする立場がオートポイエーシス理論（Maturana & Varela, 1974）で、ALifeの最初期のモデルの一つでもあります。

さらに、Zwicker, Seyboldt, Weber, Hyman, Jülicher (2016)は、相分離（液−液相分離、Cahn-Hilliard型のシャープな界面）と非平衡化学反応を組み合わせるだけで、膜がなくても液滴が成長・分裂を繰り返すことをシミュレーションで示しました（"membraneless protocells"）。生命に必要な"不連続性"は、脂質膜という特別な物質ではなく、相分離が作る鋭い界面＋非平衡の化学反応という、より一般的な組み合わせで実現できる可能性を示唆しています。

参考文献