Chaos is not an invariant of fractional-order systems: the clock shapes orbits, attractors and bifurcations

技術的要約：「カオスは分数次システムの不変量ではない：時計が軌道、アトラクター、および分岐を形作る」

問題提起
常微分方程式（整数次システム）の定性的理論において、ポアンカレ以来確立された基本的な原則は、ベクトル場に正の状態依存関数（「時計」）を乗じても、軌道の通過速度が変わるだけで、軌道自体の幾何学的形状は変わらないというものである。この再パラメータ化不変性は、相図、リミットサイクル、アトラクター、および分岐構造が、時間則とは独立してベクトル場のみに依存する性質であることを意味している。

しかし、分数次システム（特にCaputo微分によって定義されるもの）は、状態の過去の履歴を重み付けする非局所的なメモリカーネルを持つ。本論文は、この古典的な再パラメータ化不変性が分数次システムにおいても成立するかどうかを調査している。具体的には、もし分数次システム $D^\alpha x = f(x)$ が $D^\alpha x = g(x)f(x)$ （ここで $g(x)$ は正の状態依存関数）に変更された場合、それらの解は単に異なる速度で通過しているだけで、同じ幾何学的軌道を記述するのだろうか？という問いである。著者は、分数次 $\alpha \in (0, 1)$ において、この不変性は成立しないと断じている。

手法
本研究は、解析的な展開と数値的な検証を組み合わせて用いている：

1. 解析的導出: 著者は、Caputo微分方程式 $D^\alpha x = F(x)$ の解の短時間挙動を、等価なボルテラ積分方程式を用いて分析する。参照システム（$F=f$）と、時計を適用したシステム（$F=gf$）の両方について、$t^\alpha$のべき乗を用いたテイラー近似的な展開を行う。そして、変換された時間$s(t) = \int_0^t g(x(\tau)) d\tau$によって再パラメータ化された参照解の解$x(t)$ と比較することで、不一致項を分離する。
2. 普遍関数: 解析の結果、次数 $\alpha$ の関数である、偏差を定量化する2つの普遍関数 $c(\alpha)$ および $d(\alpha)$ が得られる。
3. 数値的検証: 理論的な予測を以下の4つのシステムでテストした：
   • 非線形振動子（係数 $c(\alpha)$ および $d(\alpha)$ の検証のため）。
   • 分数次ファン・デル・ポール（van der Pol）振動子（リミットサイクルの変形を調べるため）。
   • 分数次ロッスラー（Rössler）システム（周期倍分岐のカスケードと分岐閾値を調べるため）。
   • 分数次チェン（Chen）システム（カオス的レジームと周期的な窓の安定性をテストするため）。
     数値計算には、堅牢性を確保するために、分数次アダムス・バッシュフォース・モルトン予測子・修正子法と、拡散表現法（指数モードの重ね合わせ）の2つのスキームを用いた。

主な貢献と結果

1. 時計共変性の崩壊: 本論文は、$0 < \alpha < 1$ において、分数次システムの軌道は状態依存の時計に対して不変ではないことを証明している。

   • 時間則の不一致 ($t^{2\alpha}$): 最初の非消滅次数において、時計を適用した解と再パラメータ化された参照解との差は、軌道に接する方向である。この不一致は、普遍関数によって制御される：
     $$c(\alpha) = \frac{1}{\Gamma(1+2\alpha)} - \frac{\alpha}{\Gamma(1+\alpha)\Gamma(2+\alpha)}$$
     この関数は $\alpha=1$ のときのみゼロになり（整数次の定理を回復）、$\alpha \to 0$ のとき 1 に近づく。
   • 幾何学的不一致 ($t^{3\alpha}$): 次の次数では、横方向の成分が現れ、軌道が幾何学的に異なる曲線であることを証明している。これは、第2の普遍関数 $d(\alpha)$ によって制御されており、これは $\alpha \in (0, 1)$ のすべての範囲で非ゼロである。
   • 数値的確認: 係数 $c(\alpha)$ および $d(\alpha)$ は、理論的予測に対して高い精度（$c$ について誤差 < 0.1%、$d$ について誤差 < 1%）で数値的に検証された。

2. 局所安定性の不変性: 著者は、局所的な性質は不変であることを示している。平衡点、その線形安定性、およびマティオン（Matignon）の臨界次数（安定性が失われる閾値 $\alpha^*$ であり、固有値の偏角によって決定される）は、時計によって変化しない。これは、マティオンの判定基準が固有値の偏角に依存しており、正のスカラー $g(x^*)$ を乗じても影響を受けないためである。

3. 大域的ダイナミクスの非不変性: 局所的な安定性は保持されるが、大域的な構造は変形を受ける。

   • リミットサイクル: 分数次ファン・デル・ポール振動子において、時計を変更すると、中間次数においてリミットサイクルの形状と振幅が最大16%まで変形する。これは整数次では起こり得ない現象である。
     メニズム。
   • 分岐: 分数次ロッスラーシステムにおいて、時計は局所的な極大値をシフトさせるが、周期倍分岐のシーケンスやカオス開始の臨界次数を変えることはない。
   • カオス対周期性: 分数次チェンシステムにおいて、状態依存のアレニウス型時計（温度のような変数 $z$ に依存）は、同じ次数 $\alpha$ かつ同じ初期条件からであっても、周期的なレジームをカオス的なものへ（あるいはその逆に）変容させることができる。具体的には、$\alpha=0.94$ において、自明な時計の下では周期的なシステムが、特定のアレニウス型時計の下ではカオス的になる。

4. 一般則: 著者は、分岐のシフトに関する規則を確立した：

   • 局所分岐（平衡点の安定性、ホップ分岐）は移動しない。
   • 大域的分岐（周期倍分岐、クライシス、周期的な窓の生成・消滅）は、時計の強さに応じて線形にシフトする。
   • システムの「時計の変化に対する脆弱性」は、そのシステムが大域的分岐にどれだけ近いかによって決定される。

意義と主張
本論文は、分数次システムにとって「時計」は単なる数学的な便宜ではなく、モデルの物理的な構成要素であると主張している。

• 物理的解釈: 整数次モデルにおいて、ベクトル場に状態依存の因子（例：温度依存の速度）を乗じることは、時間の単位を再定義することと同等である。しかし、分数次モデルにおけるこれは、明確に異なる物理的仮説となる：

  • $D^\alpha x = g(x)f(x)$ は、メモリが物理的な時間において蓄積することを意味する。
  • $D^\alpha X = f(X)$ は、固有の時間または物質的時間（粘弾性やガラス転移の概念で使用されるもの）においてメモリが蓄積することを意味する。
    本論文は、これらが2つの異なるモデルであり、異なる予測をもたらすこと、そしてその差が普遍関数 $c(\alpha)$ および $d(\alpha)$ によって定量化されることを主張している。

• カオス理論への影響: 著者は、「カオスは分数次システムの不変量ではない」と結論付けている。最小カオス次数、分岐図、およびストレンジアトラクターの形状は、システム自体の不変量ではなく、特定の場と時計のペアの特性である。したがって、マティオンの判定基準に基づくカオスの必要条件は堅牢であるが、実際のカオスの発生は、動力学的因子の具体的な定式化に依存する。

• 実用的な帰結: 分数次モデルを用いる研究者は、メモリが物理的な時間において蓄積されるのか、あるいは固有の物質的時間において蓄積されるのかを明示的に宣言しなければならない。状態空間上に描かれた分数次アトラクターは、ベクトル場と、それを走らせるために使用される時計の両方の特性である。異なる時計を用いれば、異なるアトラクター（周期的なものかカオス的なものかといった性質の違いを含む）が得られることになる。