[論文レビュー] Boolean Irreducibility and Phase Transitions in NK-Kauffman Networks

Federico Zertuche|arXiv (Cornell University)|Aug 1, 2012

Nonlinear Dynamics and Pattern Formation被引用数 1

ひとこと要約

この論文は、構造的複雑性を捉えるブール関数の分類であるブール非可約性を組み込んだことで、NK-カウフマンネットワークの平均場解析を再考し、修正された遷移曲線を導出する。主な結果は、$ p_c = 1/2 $ における臨界結合数 $ K_c \approx 2.621 $ へのシフトであり、Derrida らの元々の $ K_c = 2 $ を置き換えるもので、数値シミュレーションが修正された閾値を支持している。

ABSTRACT

In a series of articles published in 1986 Derrida, and his colleagues studied two mean field treatments (the quenched and the annealed) for extit{NK}-Kauffman Networks. Their main results lead to a phase transition curve $ K_c \, 2 \, p_c \left( 1 - p_c ight) = 1 $ ($ 0 < p_c < 1 $) for the critical average connectivity $ K_c $ in terms of the bias $ p_c $ of extracting a $1$ for the output of the automata. Values of $ K $ bigger than $ K_c $ correspond to the so-called chaotic phase; while $ K < K_c $, to an ordered phase. In~[F. Zertuche, {\it On the robustness of NK-Kauffman networks against changes in their connections and Boolean functions}. J.~Math.~Phys. {\bf 50} (2009) 043513], a new classification for the Boolean functions, called {\it Boolean irreducibility} permitted the study of new phenomena of extit{NK}-Kauffman Networks. In the present work we study, once again the mean field treatment for extit{NK}-Kauffman Networks, correcting it for {\it Boolean irreducibility}. A shifted phase transition curve is found. In particular, for $ p_c = 1 / 2 $ the predicted value $ K_c = 2 $ by Derrida {\it et al.} changes to $ K_c = 2.62140224613 \dots $ We support our results with numerical simulations.

研究の動機と目的

• ブール非可約性の概念を統合することで、NK-カウフマンネットワークの平均場処理を再表現すること。
• Derrida らが導出した古典的遷移曲線を、ブール関数の構造的複雑性を考慮することで修正すること。
• 特に対称的バイアス $ p_c = 1/2 $ におけるブール非可約性が臨界結合数 $ K_c $ に与える影響を調査すること。
• 数値シミュレーションを用いて、修正された遷移閾値を検証すること。

提案手法

• クエンチドおよびアンナールド処理を拡張し、ブール非可約性を含む平均場アプローチを採用する。
• 低次元の成分に簡略化できない関数を区別する、ブール関数の非可約性に基づく分類を導入する。
• 平均バイアス $ p_c $ と非可約性補正済み結合数を組み込んだことで、再び遷移条件を導出する。
• 非可約性を補正した修正された遷移曲線を導出：$ K_c \cdot 2 \cdot p_c (1 - p_c) = 1 $ であり、これにより $ K_c $ が上昇する。
• 特に $ p_c = 1/2 $ における $ K_c $ のシフトを確認するため、理論的予測を数値シミュレーションで検証する。

実験結果

リサーチクエスチョン

• RQ1 ブール非可約性は、NK-カウフマンネットワークの平均場解析にどのように影響を与えるか？
• RQ2 ブール非可約性を段階的遷移モデルに組み込むと、再評価された臨界結合数 $ K_c $ はどのようになるか？
• RQ3 $ p_c = 1/2 $ において、遷移閾値は顕著にシフトするか？もしそうなら、その大きさは？
• RQ4 数値シミュレーションは、非可約性によって生じる高い $ K_c $ の理論的予測を確認できるか？

主な発見

• ブール非可約性を平均場処理に組み込むことで、古典的 Derrida らのモデルとは異なったシフトされた遷移曲線が得られる。
• $ p_c = 1/2 $ における再評価された臨界結合数は $ K_c = 2.62140224613\dots $ であり、元の $ K_c = 2 $ より顕著に高い。
• $ K_c $ のシフトは、ブール非可約性が捉える構造的複雑性に起因し、混沌沌とした状態に達するための有効なネットワーク結合数を増加させる。
• 数値シミュレーションは理論的予測を支持し、混沌沌とした状態の閾値が上昇していることを確認した。
• 結果から、バイアスだけでなく、ブール関数の構造そのものが、NKネットワークの力学的状態を決定づける重要な要因であることが示された。