Gray-Scott 散逸構造 / Dissipative structures

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

How to use the reaction–diffusion simulator, and an introduction to the non-equilibrium thermodynamics and dissipative-structure theory behind it.

使い方 / How to use

• プリセット / Presets：Mitosis・Coral・Worms など、Pearsonの分類にある代表的なパターンにワンクリックで切り替えられます。 / Switch with one click to representative patterns from Pearson's classification, such as Mitosis, Coral and Worms.
• Feed rate F / Kill rate k：反応の根本パラメータ。スライダーで自由に動かせます。境界線のすぐ近くで劇的にパターンが変わることがあります。 / The fundamental parameters of the reaction. Move them freely with the sliders; the pattern can change dramatically right next to a boundary line.
• キャンバス上のドラッグ / Drag on the canvas：クリックしたまま動かすと、その場所に新しい"種"（v成分）を追加できます。 / Hold and move to add a new "seed" (component v) wherever you drag.
• 系を閉じる / Close the system：ONにするとFを強制的に0にし、外部からの物質流入を止めます。パターンが崩壊していく様子は、開放系であることが構造維持の必要条件であることの直接的な観察です。 / When on, F is forced to 0 and the inflow of material from outside stops. Watching the pattern collapse is a direct observation that being an open system is a necessary condition for keeping the structure.
• エントロピー / 散逸パネル / Entropy and dissipation panel：下記「HとDの意味」を参照。 / See "H and D" below.
• 🔬 全プリセットを自動計測 / Auto-measure all presets：8つのプリセットを自動で切り替えながら、定常状態のH・Dを計測してテーブル化・CSV出力します。 / Cycles through the 8 presets automatically, measures H and D in the steady state, and outputs a table and CSV.

Gray-Scott方程式とは / The Gray–Scott equations

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

This is a system in which the concentration fields of two chemicals, u (substrate) and v (catalyst), react while diffusing. The autocatalytic reaction term uv² (v catalytically increases v) is the driving force for spatial pattern formation, just as in Turing patterns.

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

F(1−u) is the rate at which u flows into the system from outside (the source term of an open system), and (F+k)v is the rate at which v is removed from the system (the sink term). These F and k are what make the Gray–Scott system a non-equilibrium open system; the simulator's "Close the system" toggle simply sets this F to 0.

なぜ「散逸構造」なのか / Why "dissipative structures"?

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

Ilya Prigogine, who won the 1977 Nobel Prize in Chemistry, developed the theory that in an open system far enough from thermodynamic equilibrium, nonlinear interactions can stably maintain self-organized structures (dissipative structures) while discarding entropy into the environment. Bénard convection cells and the Belousov–Zhabotinsky reaction, described below, are classic examples.

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

Gray–Scott patterns can be understood in the same framework. The pattern persists only while the system is open (F>0); when you close the system (F=0), the inflow of material stops and the pattern gradually breaks down, relaxing toward thermodynamic equilibrium (a uniform state). You can confirm this yourself by flipping the toggle in the simulator.

ベロウソフ・ジャボチンスキー反応との関係 / Relation to the Belousov–Zhabotinsky reaction

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

In fact, what directly inspired Prigogine to build the theory of dissipative structures was a real chemical reaction, the Belousov–Zhabotinsky (BZ) reaction. Spread thinly in a petri dish, it shows target patterns and rotating spiral waves as actual changes in color. Mathematically there is a reaction–diffusion model devoted to it, the "Oregonator," which belongs to the same "reaction term plus diffusion term" framework as Gray–Scott.

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

And if you run the BZ reaction in an ordinary flask (a batch reactor), the reagents are used up and the oscillations stop. That is exactly the same phenomenon as the "Close the system" toggle. To keep oscillations and patterns going indefinitely you need a device that keeps feeding in reagents, a CSTR (continuous stirred-tank reactor), which is the chemical counterpart of F.

Pearsonの分類とプリセット / Pearson's classification and the presets

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

Ever since Pearson (1993) sorted the Gray–Scott F–k parameter space into about twelve classes labeled with Greek letters, it has been explored thoroughly around the world for more than 30 years (thousands of coordinates are published in MROB's archive). The presets in this simulator are representative coordinates from that work.

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

Even with the same (F, k), the behavior can change with the numerical discretization (grid resolution, Laplacian weights, time step). This is a widely known property in the field and not an implementation bug.

HとD：エントロピーと散逸速度の指標 / H and D: indicators of entropy and dissipation

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

This simulator computes two kinds of statistics in real time. It is important not to confuse the two.

H：パターンのShannonエントロピー / H: Shannon entropy of the pattern

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

The values of the v field are put into a 32-bin histogram, H is computed as H = −Σp·log(p) and normalized by log(32), giving an index of spatial "disorder." When structure (spots, spirals) forms, H actually goes down. This is the same picture as Schrödinger's idea that life "feeds on negative entropy": local order (low entropy) is achieved by pushing a larger increase of entropy onto the environment.

D：散逸(流束)速度の近似 / D: an approximation of the dissipation (flux) rate

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

D ≈ mean[F(1−u) + kv]. It is the size of the material flux through the system, and serves as an approximate index of the rate of entropy production in a non-equilibrium steady state. Turning on "Close the system" makes F=0, and D drops sharply.

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

In our experiments, H (pattern complexity) and D (dissipation rate) tended to rank the presets in roughly the same order (Waves > Chaos > Holes > Spots > Coral > Solitons≈Worms > Mitosis). However, at root both reflect "what fraction of the domain is reacting," so they are not independent and should be interpreted with care.

Maximum Entropy Production Principle (MEPP) との関連 / Relation to the Maximum Entropy Production Principle (MEPP)

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

The hypothesis (MEPP) that "in an open system far from equilibrium, the dissipative-structure mode with the largest rate of entropy production is the one that is stabilized and selected" is still debated in non-equilibrium thermodynamics. Jeremy England's "dissipative adaptation" theory is a modern development that derives a related idea from more rigorous stochastic thermodynamics. Note, however, that MEPP is not a proven law: counterexamples and limits of applicability have been pointed out, and it remains a hypothesis.

膜・区画化という視点 / Membranes and compartmentalization

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

The concentration fields of Gray–Scott and Lenia are continuous, and their boundaries are only gradations of the same substance. In contrast, autopoiesis theory (Maturana & Varela, 1974) takes the essence of life to be "a system that separates itself from the outside with a membrane, while its own internal metabolic network keeps producing that membrane." It was also one of the earliest models in ALife.

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

Furthermore, Zwicker, Seyboldt, Weber, Hyman and Jülicher (2016) showed in simulation that simply combining phase separation (liquid–liquid phase separation, with sharp Cahn–Hilliard-type interfaces) and non-equilibrium chemical reactions makes droplets grow and divide repeatedly, even without a membrane ("membraneless protocells"). This suggests that the "discontinuity" life needs may be achieved not through a special material such as a lipid membrane, but through the more general combination of sharp interfaces created by phase separation and non-equilibrium chemical reactions.

参考文献 / References

• J. E. Pearson, "Complex Patterns in a Simple System," Science, 1993.
• Reaction-Diffusion by the Gray-Scott Model: Pearson's Parameterization (MROB)
• I. Prigogine, Nobel Prize in Chemistry 1977 lecture, "Time, Structure and Fluctuations."
• Belousov-Zhabotinsky reaction — Scholarpedia
• R. Swenson (1989); G. Dewar (2003) — Maximum Entropy Production Principle
• J. England, "Statistical physics of self-replication," 2013（dissipative adaptation）
• H. Maturana & F. Varela, "Autopoiesis: the organization of living systems," 1974.
• D. Zwicker et al., "Growth and Division of Active Droplets: A Model for Protocells," arXiv:1603.01571, 2016.

シミュレーターで試す / Try the simulator →

hamzWorkerz（ハムズワーカーズ） ・ プライバシーポリシー / Privacy