反応拡散系

反応拡散系（はんのうかくさんけい、英: reaction-diffusion system）とは、空間に分布された一種あるいは複数種の物質の濃度が、物質がお互いに変化し合うような局所的な化学反応と、空間全体に物質が広がる拡散の、二つのプロセスの影響によって変化する様子を数理モデル化したものである。

反応拡散系は、化学の分野において自然な形で応用されるものである。しかし、化学的ではない動力学過程を表現する上でも、反応拡散系は応用される。例えば、生物学や地質学、物理学や生態学において、そのような応用例は見られる。数学的に言うと、反応拡散系は半線形放物型偏微分方程式の形を取るものである。

一般的には次のように記述される。

${\displaystyle \partial _{t}{\boldsymbol {q}}={\underline {\underline {\boldsymbol {D}}}}\,\nabla ^{2}{\boldsymbol {q}}+{\boldsymbol {R}}({\boldsymbol {q}}).}$

ここでベクトル q(x,t) の各成分はある物質の濃度を表し、${\displaystyle {\underline {\underline {D}}}}$ は拡散係数からなる対角行列を表し、R はすべての局所的な反応を表す。反応拡散方程式の解は、進行波の形成や波に似た現象、あるいは帯や六角形のような自己組織パターン（英語版）、あるいは散逸ソリトン（英語版）のようなより複雑な構造を含む、幅広い範囲の挙動を見せるものである。

最も簡単な反応拡散方程式は、空間一次元における単一物質の濃度 u に関するもので、

${\displaystyle \partial _{t}u=D\partial _{x}^{2}u+R(u)}$

と記述され、これはKPP方程式（Kolmogorov-Petrovsky-Piskounov の略）とも呼ばれる^[1]。反応項が無い場合、方程式は純粋な拡散過程のみを表す。そのような方程式はフィックの第二法則に関係するものである。反応項が R(u) = u(1-u) である場合、生物学的な人口の広がりを表現するために元々用いられた、フィッシャーの方程式が得られる^[2]。R(u) = u(1 − u²) である場合、レイリー＝ベナール対流を表すためのニューウェル＝ホワイトヘッド＝シーゲル方程式が得られる^[3][4]。R(u) = u(1 − u)(u − α) および 0 < α < 1 である場合には、燃焼理論に現れるより一般的なゼルドビッチ方程式が得られる^[5]。そしてその特別な退化的な例は R(u) = u² − u³ の場合に得られ、その方程式もまたゼルドビッチ方程式と呼ばれる^[6]。

一成分の系のダイナミクスは、ある特定の制限に関するものである。なぜならば、その発展方程式は変分系

${\displaystyle \partial _{t}u=-{\frac {\delta {\mathfrak {L}}}{\delta u}}}$

としても書かれ、したがってこれは次式で与えられる「自由エネルギー」${\displaystyle {\mathfrak {L}}}$ の永続的な減少を意味するからである。

${\displaystyle {\mathfrak {L}}=\int \limits _{-\infty }^{\infty }\left[{\frac {D}{2}}(\partial _{x}u)^{2}-V(u)\right]{\text{d}}x}$

ここで V(u) は R(u)=dV(u)/du であるようなポテンシャルを表す。

一つ以上の定常同次解を備える系において、典型的な解は、その同次状態をつなぐ進行波として与えられる。そのような解は、その形状を変えずに一定の速度で移動し、u(x, t) = û(ξ) と記述される。ここで ξ = x − ct であり c はその進行波の速度を表す。ここで、進行波は一般的に安定な構造を備えるが、非単調な定常解（例えば、前進と反前進のペアで構成される局所化された領域）は不安定であることに注意することである。c = 0 の場合には、この記述内容には次のような簡単な証明が存在する^[7]：u₀(x) が定常解、u=u₀(x) + ũ(x, t) が無限小摂動解であるなら、線型安定性解析によって次の方程式が導かれる。

${\displaystyle \partial _{t}{\tilde {u}}=D\partial _{x}^{2}{\tilde {u}}-U(x){\tilde {u}},\quad U(x)=-R^{\prime }(u)|_{u=u_{0}(x)}.}$

この解 ũ = ψ(x)exp(−λt) に対し、シュレディンガー型の固有値問題

${\displaystyle {\hat {H}}\psi =\lambda \psi ,\qquad {\hat {H}}=-D\partial _{x}^{2}+U(x),}$

が得られる。ただし、その負の固有値が解の不安定性に帰結するものである。平行移動不変性により、ψ = ∂_xu₀(x) は固有値 λ = 0 に対応する中立的な固有関数であり、その他のすべての固有関数は、ゼロ解の数について単調に増加する実固有値の絶対値について、増加する結び目の数に従って分類される。固有関数 ψ = ∂_x u₀(x) は少なくとも一つのゼロ解を持ち、非単調な定常解については対応する固有値 λ = 0 は最小のものではなく、したがって不安定性を意味する。

進行波の速度 c を決定するために、移動座標系を考え、定常解を探すことが出来る。

${\displaystyle D\partial _{\xi }^{2}{\hat {u}}(\xi )+c\partial _{\xi }{\hat {u}}(\xi )+R({\hat {u}}(\xi ))=0.}$

この方程式は、位置 û、時間 ξ、力 R、減衰係数 c に対する質量 D の動きに対する機械的な類似性を備えるものである。

空間一次元からより高次の空間次元に議論を移しても、依然として有効となる内容は数多く存在する。平らな、あるいは曲がった進行波は典型的な構造で、曲がった波の局所速度が局所曲率半径に依存するに従い、新たな効果が生じるものである（このことは極座標系を考えることで分かる）。この現象は、いわゆる曲率駆動不安定性を導く^[8]。

二成分の系は、一成分の系と比較してより幅広い現象を許すものである。アラン・チューリングによって初めて提唱されたある重要なアイデアに、局所的な系においては安定であっても拡散の存在する状況では不安定となる状態というものがある^[9]。拡散は一般的には安定化効果と関連するものであるので、一聴するとこのアイデアは直感に反するもののようでもある。

しかしながら、線型化安定性解析によって、一般的な二成分系

${\displaystyle \left({\begin{array}{c}\partial _{t}u\\\partial _{t}v\end{array}}\right)=\left({\begin{array}{cc}D_{u}&0\\0&D_{v}\end{array}}\right)\left({\begin{array}{c}\partial _{xx}u\\\partial _{xx}v\end{array}}\right)+\left({\begin{array}{c}F(u,v)\\G(u,v)\end{array}}\right)}$

を線型化するとき、定常同次解の平面波摂動

${\displaystyle {\tilde {\boldsymbol {q}}}_{\boldsymbol {k}}({\boldsymbol {x}},t)=\left({\begin{array}{c}{\tilde {u}}(t)\\{\tilde {v}}(t)\end{array}}\right)e^{i{\boldsymbol {k}}\cdot {\boldsymbol {x}}}}$

は次を満たすことが分かる。

${\displaystyle \left({\begin{array}{c}\partial _{t}{\tilde {u}}_{\boldsymbol {k}}(t)\\\partial _{t}{\tilde {v}}_{\boldsymbol {k}}(t)\end{array}}\right)=-k^{2}\left({\begin{array}{c}D_{u}{\tilde {u}}_{\boldsymbol {k}}(t)\\D_{v}{\tilde {v}}_{\boldsymbol {k}}(t)\end{array}}\right)+{\boldsymbol {R}}^{\prime }\left({\begin{array}{c}{\tilde {u}}_{\boldsymbol {k}}(t)\\{\tilde {v}}_{\boldsymbol {k}}(t)\end{array}}\right).}$

チューリングのアイデアは、反応函数のヤコビアン R' の符号によって特徴付けられた系の四つの同値類においてのみ、理解されるものである。特に、有限の波ベクトル k が最も不安定なものであると仮定されたとき、そのヤコビアンは符号

${\displaystyle \left({\begin{array}{cc}+&-\\+&-\end{array}}\right),\quad \left({\begin{array}{cc}+&+\\-&-\end{array}}\right),\quad \left({\begin{array}{cc}-&+\\-&+\end{array}}\right),\quad \left({\begin{array}{cc}-&-\\+&+\end{array}}\right)}$

を備えるものでなければならない。この系の類は、その第一の描写にちなみ活性因子・抑制因子系（activator-inhibitor system）と呼ばれる。すなわち、基底状態の近くではある成分は両成分の生産を促進するが、一方で別の成分はそれらの成長を阻害している。その最も傑出した代表例は、フィッツフュー＝南雲方程式

${\displaystyle {\begin{aligned}\partial _{t}u&=d_{u}^{2}\,\nabla ^{2}u+f(u)-\sigma v,\\\tau \partial _{t}v&=d_{v}^{2}\,\nabla ^{2}v+u-v\end{aligned}}}$

である。ここで ƒ(u) = λu − u³ − κ は活動電位がどのように神経を移動するかを表している ^[10][11]。また d_u、d_v、τ、σ および λ は正定数である。

活性因子・抑制因子系にパラメータの変化が施されたとき、均質な基底状態が安定であるような条件から、それが線型不安定であるような条件へと移ることがある。その対応する分岐は、支配的な波数 k = 0 を備える大域的な振動均質状態へのホップ分岐であるか、支配的な有限の波数を備える大域的なパターン状態へのチューリング分岐のいずれかであり得る。空間二次元における後者の分岐は、通常、ストライプや六角形のパターンを導くものである。

• 亜臨界チューリング分岐：フィッツフュー＝南雲型の二成分反応拡散系におけるノイズの多い初期状態からの六角形パターンの形成。
• 
  t = 0のノイズの多い初期状態。
• 
  t = 10の系状態。
• 
  t = 100のほとんど収束した状態。

フィッツフュー＝南雲の例に対し、そのチューリング分岐およびホップ分岐のための線型安定領域の境界を作る中立安定曲線は、次式で与えられる。

${\displaystyle {\begin{aligned}q_{\text{n}}^{H}(k):&{}\quad {\frac {1}{\tau }}+(d_{u}^{2}+{\frac {1}{\tau }}d_{v}^{2})k^{2}&=f^{\prime }(u_{h}),\\[6pt]q_{\text{n}}^{T}(k):&{}\quad {\frac {\kappa }{1+d_{v}^{2}k^{2}}}+d_{u}^{2}k^{2}&=f^{\prime }(u_{h}).\end{aligned}}}$

分岐が亜臨界であるなら、基底状態とパターンが共存するようなヒステリシスな領域において、しばしば局所的な構造（散逸ソリトン（英語版））が観測される。その他、頻繁に現れる構造としては、パルス列（周期進行波としても知られる）、螺旋波、ターゲットパターンがある。それら三つの解のタイプは、局所的なダイナミクスが安定なリミットサイクルを備えるような二成分（あるいはより多くの成分）の反応拡散方程式の、一般的な構造である^[12]。

• フィッツフュー＝南雲型の二成分の拡散反応系に現れる他のパターン。
• 
  回転する螺旋。
• 
  ターゲットパターン。
• 
  定常的な局所化されたパルス（散逸ソリトン）。

様々な系に対して、二つよりも多い成分の反応拡散方程式が提唱されている。例えば、ベロウソフ・ジャボチンスキー反応のモデル^[13]、血液凝固のモデル^[14]あるいは平面の気体放電（英語版）系^[15]などが挙げられる。

より多くの成分を含む系では、一成分あるいは二成分の系では起こり得ないさまざまな現象（例えば、大域的フィードバックのない空間多次元における安定ランニングパルスなど）が起こることが知られている^[16]。扱う系の性質に依存して起こり得る現象についての導入と系統的な概要については、^[17]で与えられている。

近年、反応拡散系はパターン形成（英語版）に対する基本的なモデルとして、多くの関心を集めている。上述のパターン（進行波、スパイラル、ターゲット、六角形、ストライプ、散逸ソリトン）は様々なタイプの反応拡散系において見られる。しかしそこには多くの矛盾、例えば、局所反応項においてそれらが見られるなど、が存在する。反応拡散過程は、生物学における形態形成と関連する過程に対する本質的な基盤であることも、述べられている^[18] 。そして、それは動物の毛皮や、皮膚の色素沈着に対しても関連付けられている^[19][20]。

反応拡散方程式のさらなる応用例は、生態の侵入^[21]や感染症の拡がり^[22]、腫瘍の成長^[23][24][25]や、傷の治癒^[26]などに見られる。反応拡散系が関心を集める別の理由には、それらが非線型偏微分方程式であるにもかかわらず、解析的な扱いがたびたび可能となることが挙げられる^{[7][8][27][28][29]}。

化学における反応拡散系のよく管理された実験には、現在三つの方法があることが知られている。第一に、ゲル型リアクター^[30]あるいはキャピラリーチューブ^[31]が用いられること。第二に、キャタリティック表面上の温度パルスが調べられること^[32][33]。第三に、反応拡散系を用いて神経パルスの進行がモデル化されること^[10][34]、である。

それらの一般的な例とは別に、適切な環境下ではプラズマ^[35]や半導体^[36]のような電気輸送系も、反応拡散の手法によって表現することが出来ることが判明している。それらの系に対して、パターン形成に関する様々な実験が行われている。

• Java applet 反応拡散のシミュレーション
• Another applet グレイ＝スコットの反応拡散
• Java applet いくつかの種の蛇に見られるパターン形成をシミュレートするための反応拡散
• Gallery 反応拡散の画像と動画
• TexRD software 反応拡散を基盤とするランダム・テクスチャー・ジェネレーター
• Reaction-Diffusion by the Gray-Scott Model: Pearson's parameterization グレイ＝スコット反応拡散のパラメータ空間の可視図
• A Thesis 反応拡散パターンと分野の概観についての論文
• ReDiLab - Reaction Diffusion Laboratory ベロウソフ＝ジャボチンスキー、グレイ＝スコット、ウィラモフスキー＝レスラーおよびフィッツフュー＝南雲モデルをシミュレートした Flash & GPU based application（完全なソースコード付き）