Emergent Fibonacci Structure from Local Growth Rules

技術要約：局所的な成長規則からのフィボナッチ構造の創発

問題提起
植物の葉序から準結晶に至るまで、螺旋状の形態は自然界に遍在している。既存の理論的枠組み（連続的な反応拡散系、生物力学、および幾何学的モデル）はこれらのパターンを記述しているが、離散的なフィボナッチ数の出現については完全には解明されていない。具体的には、なぜフィボナッチ数列が、生物学的な指示として明示的にエンコードされたり、あるいは大域的なエネルギー最小化によって課されたりすることなく、これほど一貫して現れるのかという点が不明である。従来のアプローチは、記号的な置換や幾何学的なインフレーションを通じてフィボナッチ構造を生成しており、これは数列をモデルに組み込んでいる（構築している）のであって、導出しているわけではない。本論文は、純粋に局所的な離散的成長規則と最小限のエネルギー制約から、フィボナッチ構造がどのように創発し得るかという問いに取り組むものである。

手法
著者は、以下の3つの物理的な動機付けに基づいた、二次元の離散細胞成長モデルを提案している：

1. 成長報酬： 面積の増加はエネルギー的に有利である。
2. 回転ペナルティ： 成長方向を変更することは、新しい面に沿って成長を支えるために再活性化が必要となる不活性細胞の数に比例したエネルギーコストを課す。
3. 制約された逐次的拡張： 成長は各ステップにおいて、単一の外向きの辺に制限される。

このモデルは、生物学的な細胞形成の離散的な性質を反映するため、完全に整数内で動作する。成長セグメントは、成長方向の長さ $g_i$ と垂直方向の長さ $p_i$ のペア $(g_i, p_i)$ によって定義される。幾何学的には、右回転の際に漸化式 $p_{i+1} = g_i + p_{i-1}$ を強制する。

回転するか直進するかという決定は、エネルギー不等式 $E = -\lambda V + \rho Q$ によって決定される。ここで、$V$ は面積利得、$Q$ は回転ペナルティである。直進を続ける場合と回転する場合のエネルギーを比較することで、本モデルは垂直方向の数列 $\{p_i\}$ に関する一般化された漸化式を導出する：
$$p_{i+1} = \lceil p_i + (\kappa - 1)p_{i-1} \rceil + p_{i-1}$$
ここで、$\kappa = \rho/\lambda$ は、回転ペナルティの重みと成長報酬の重みの比である。モデルの解析には、短期的挙動を測定するための連分数と、長期的挙動を測定するための漸近的成長比が用いられる。

主な貢献

• 賦課なしでの創発： 本論文は、面積増加のエネルギー的利益が回転のコストと均衡する臨界点（$\kappa = 1$）において、フィボナッチ漸化式（$p_{i+1} = p_i + p_{i-1}$）が正確に現れることを示している。これは、フィボナッチ数列を事前に定義したり、大域的な最適化を用いたりすることなく実現される。
• パラメータ化されたスパイラル・タイリングの族： 本モデルは、$\kappa$ によってパラメータ化された、連続的な整数列および対応する長方形螺旋タイリングを生成する。
• 特殊な数列の特定： 解析により、特定の整数値 $\kappa$ において既知の整数列を回収する：
  • $\kappa = 1$：フィボナッチ数。
  • $\kappa = 2$：ヤコブスタール数（$p_i$ について）および2の冪乗（$g_i$ について）。
  • $\kappa = 3$：OEIS A006130。
  • 非整数 $\kappa$：天井関数によるステップ依存の丸め処理が導入されているため、ホラダム族には適合しない、真に新しい整数列を生成する。
• 過渡的なレオナルド構造： 臨界点（$\kappa = 1 + \epsilon$）の近傍において、モデルは新たな漸近的成長率へと移行する前に、過渡的なレオナルド数（OEIS A00 ever1595）の数列を生成する。

結果

• 臨界点と相転移： $\kappa = 1$ において、モデルはフィボナッチ数列を正確に再現し、漸近的成長比 $r(1) = \phi \approx 1.618$ を持つ。
• 漸近的挙動： 一般的な $\kappa$ に対して、漸近的成長比は $r(\kappa) = (1 + \sqrt{1 + 4\kappa})/2$ で与えられる。$\kappa$ が1から逸脱すると、成長率は変化し、数列の記号的構造もシフトする。
• 連分数指標： 初期の項を重く評価する連分数指標は、$\kappa = 1$ で不連続な低下を示す。$\kappa < 1$ の場合、数列はより長い期間、フィボナッチ的であり続ける。$\kappa > 1$ の場合、数列は急速に異なる記号的構造（初期にはレオナルド数）へと移行し、指標値の急激な変化を引き起こす。
• 数値的検証： シミュレーションにより、$\kappa = 1$ のとき、数列がフィボナッチ数を正確に再現することが確認された。$\kappa = 2$ のとき、比は2に収束する。非整数 $\kappa$ の場合、理論的な $r(\kappa)$ への収束はより遅いが、明白である。

意義と主張
本論文は、有限の範囲を持つ生物学的システムにおけるフィボナッチ螺旋の普及は、**大域的な最適性ではなく、堅牢性（ロバストネス）**によって説明できると主張している。モデルの初期段階の挙動は、広範なパラメータ空間（具体的には $\kappa \leq 1$）においてフィボナッチ的であり続けるため、有限の生物はフィボナッチ・パターンを示すために正確な臨界点を占める必要はない。本モデルは、局所的なエネルギーバランスが不可欠な選択メカニズムであることを示唆しており、これにより、生物学的にエンコードされた指示なしに、単純な物理原則からフィボナッチ構造が自然に創発することを可能にしている。

本研究は、ドゥアディ＝クーデルのパラダイムおよび関連する離散的アプローチと整合しているが、これを簡略化された整数ベースのモデルを通じて達成している。フィボナッチ・タイリングを、より広い螺旋形態の族の中で、エネルギー的に特権的かつ構造的に堅牢なケースとして特定している。これは、自然界におけるこれらのパターンの普遍性に、潜在的な説明を与えるものである。