Linearization Principle: The Geometric Origin of Nonlinear Fokker-Planck Equations

本論文は、$dy/dx = y^q$ という成長則から出発し、$q$-対数座標系を用いる幾何学的アプローチにより、非線形ドリフト項や人為的な制約に依存することなく、非線形 Fokker-Planck 方程式とその定常状態である$q$-ガウス分布を導出する理論的枠組みを提示しています。

要約

1. 従来の「困った問題」：無理やりねじ曲げた地図

これまで、物理学者たちは「複雑なシステム（乱気流や生物の動きなど）」を説明するために、**「非線形フォッカー・プランク方程式」**という数式を使ってきました。

しかし、この数式には大きな欠点がありました。

• 問題点： 現象を説明するために、**「粒子同士が互いの存在に影響し合う」という、物理的に不自然な力（非線形なドリフト力）**を無理やり数式に組み込んでいたのです。
• 例え話：
  想像してください。あなたが山を登っているとき、**「自分の背丈が高くなると、重力が自分自身で強くなる」**なんていう法則があったらどうでしょう？それは不自然ですよね。でも、従来の理論は「確率（粒子のいる場所の密度）が高いところでは、粒子が自分自身を押し返す力」のような、少し奇妙な仮定をしないと現象が説明できませんでした。

2. この論文の「画期的な発見」：正しい「地図」を描き直せ！

著者の須江里（スヤリ）さんは、**「力（ドリフト）を変えて無理やり合わせるのではなく、そもそも『地図の目盛り（座標系）』を間違えていたのではないか？」**と考えました。

これがこの論文の核心である**「線形化の原理（Linearization Principle）」**です。

• 新しい視点：
  粒子の動きは、実は**「$y$ の成長速度は $y$ の $q$ 乗に比例する」**という単純な法則（$dy/dx = y^q$）に従っています。

  • $q=1$ の場合： 普通の指数関数的な成長（標準的な物理）。
  • $q \neq 1$ の場合： 複雑なシステム特有の「べき乗則（パワー・ロー）」の成長。

• 解決策：
  この成長法則に従う世界では、**「$q$-対数（$q$-logarithm）」という特殊な目盛りで測る方が、世界は実は「直線的（シンプル）」**に見えるのです。

  • 例え話：
    地球儀（球体）を平らな地図に描こうとすると、極地が引き伸ばされて歪んで見えますよね。でも、もしあなたが「地球儀そのもの」を直接扱えるなら、そこには歪みはありません。
    この論文は、**「複雑な現象は、$q$-対数という『正しい地球儀（自然な座標系）』で見れば、実は単純な直線的な動きだった」**と説いています。

3. 驚きの結果：シンプルさの維持

この「正しい地図（$q$-対数）」を使うと、どんなすごいことが起きるでしょうか？

1. 力はシンプルのまま：
   粒子にかかる「外力（重力やバネの力など）」は、相変わらずシンプルで直線的なままです。粒子が自分自身を押し返すような不自然な力は不要になりました。
2. アインシュタインの関係式が守られる：
   温度と拡散の関係を表す有名な「アインシュタインの関係式」が、複雑な修正なしにそのまま成り立ちます。
3. 「$q$-ガウス分布」が自然に生まれる：
   平衡状態（落ち着き状態）になると、粒子の分布は**「$q$-ガウス分布」という、しっぽが長い（または短い）特徴的な形になります。これは、従来の理論が「無理やり作ろうとした」形が、実は「自然な地図を使えば自然に導き出される」**ことを意味します。

4. 二面性（デュアリティ）：動きと静けさの不思議な関係

論文のもう一つの面白い発見は、「動きの指数 $q$」と「熱力学の指数 $2-q$」という双子のような関係です。

• 動き（ダイナミクス）： 粒子がどう動くかは、指数 $q$ で決まります。

• 静けさ（熱力学）： 最終的にどう落ち着くか（エネルギーやエントロピー）は、$2-q$ という別の指数で決まります。

• 例え話：
  川の流れ（動き）が激しくても（$q$）、川底の地形（熱力学）は、それとは異なるルールで安定しているようなものです。この論文は、この「動きと静けさのバランス」が、**「補正分布（エスコート分布）」**というごまかしを使わずに、自然に説明できることを証明しました。

5. 具体的な応用：バネと自由な粒子

この理論を実際に当てはめてみました。

• バネの振動（調和振動子）：
  バネに繋がれた粒子は、$q$-ガウス分布という形に落ち着きます。$q>1$ なら「しっぽが長く」、$q<1$ なら「範囲が限定された」形になります。
• 自由な粒子（拡散）：
  何もない空間を粒子が動くとき、この理論は**「異常拡散」**（通常の拡散より速かったり遅かったりする現象）を自然に説明します。
  • $q=1$：普通の拡散（コーヒーのミルクが広がるような感じ）。
  • $q<1$：遅い拡散（スポンジの中を水が染み込むような感じ）。
  • $q>1$：速い拡散（乱気流の中で粉が舞うような感じ）。

まとめ：なぜこれが重要なのか？

この論文は、**「複雑な現象を説明するために、物理法則を無理やりねじ曲げる必要はない」**と教えてくれます。

• 従来の考え方： 「現象が変だから、力も変なことにしよう」。
• この論文の考え方： 「現象が変に見えるのは、見る目（座標系）が合っていないから。正しい目盛り（$q$-対数）を使えば、世界はシンプルで美しい直線になる」。

これは、「複雑系」という謎を解く鍵が、実は「幾何学（図形や座標の性質）」にあったという、非常にエレガントで美しい発見です。これにより、量子流体や量子コンピューターなど、さらに高度な分野への応用も期待されています。

技術概要

論文要約：線形化の原理：非線形 Fokker-Planck 方程式の幾何学的起源

著者: Hiroki Suyari (千葉大学大学院情報学研究科)
日付: 2026 年 3 月 2 日

1. 背景と問題提起

複雑系において観測される異常拡散やべき乗則分布（Power-law distributions）の動的基盤は、統計物理学における中心的な課題です。これらの現象の静的性質は Tsallis エントロピーなどの一般化エントロピーでよく記述されますが、それらを生成する Fokker-Planck 方程式（FPE）の正確な形式については議論が続いていました。

従来の非線形 FPE の導出アプローチには、以下の物理的なジレンマが存在しました：

• 定常解の整合性: 観測される q-ガウス分布を定常解として得るために、状態依存の拡散係数や、現象論的な非線形ドリフト項（確率密度そのものに依存する力）を必要とすることが多い。
• 物理的意味の曖昧さ: 非線形ドリフト項の導入は、外部ポテンシャルからの粒子の独立性を損なう。また、「エスコート分布（escort distributions）」のような補助的な確率測度に依存することで、観測量の物理的意味を不明瞭にする。

本研究は、これらのジレンマを解決し、非線形 FPE に一貫した動的および幾何学的基盤を提供することを目的としています。

2. 方法論：線形化の原理（Linearization Principle）

本論文は、特定のエントロピー汎関数を仮定するのではなく、状態変数の「基本的な成長則（Fundamental Growth Law）」に根ざした線形化の原理を適用します。

• 基本成長則: 状態変数 $y(x)$ が以下の非線形成長則に従うと仮定します。
  $$\frac{dy}{dx} = y^q$$
  ここで、$q=1$ は標準的な指数関数的成長、$q \neq 1$ は複雑系に普遍的なべき乗則挙動を表します。
• 自然座標系: この成長則を線形化するために、**q-対数（q-logarithm）**を自然な座標系として導入します。
  $$z := \ln_q y = \frac{y^{1-q} - 1}{1-q}$$
  この座標系において、物理法則（拡散やドリフト）は線形な形式をとるべきだとする原理です。
• 化学ポテンシャルの一般化: 標準的な統計力学における化学ポテンシャル $\mu = k_B T \ln p + V(x)$ に対し、自然座標 $\ln_q p$ を用いて以下のように一般化します。
  $$\mu_q := k_B T \ln_q p + V(x)$$
  これにより、ドリフト項は確率密度 $p$ に対して線形に保たれつつ、拡散項に非線形性が自然に生じます。

3. 主要な成果と導出結果

A. 非線形 Fokker-Planck 方程式（NFPE）の導出

上記の原理に基づき、連続の方程式と結合させることで、以下の NFPE を導出しました。
$$\frac{\partial p}{\partial t} = \nabla \cdot \left( D p^{1-q} \nabla p + \frac{1}{\gamma} p \nabla V \right)$$

• 拡散項: $p^{1-q}$ の因子により非線形になります（多孔質媒体方程式 PME と同等）。
• ドリフト項: 外部ポテンシャル $V(x)$ に対する応答は $p \nabla V$ となり、線形のまま保たれます。これにより、アインシュタインの関係式 $D = k_B T / \gamma$ が標準的な形で維持されます。
• 既存モデルとの違い: 従来のアプローチ（Plastino & Plastino, Tsallis & Bukman など）が現象論的に非線形ドリフトを導入したのに対し、本アプローチは純粋に幾何学的原理から導出されるため、虚構の非線形力を必要としません。

B. 定常解と q-ガウス分布

定常状態（$J=0$）を解くと、以下の分布が得られます。
$$p_{st}(x) = \frac{1}{Z_q} \exp_q \left( -\beta_q V(x) \right)$$

• 調和振動子ポテンシャル $V(x) \propto x^2$ の場合、これはq-ガウス分布となります。
• $q > 1$ ではべき乗則の裾（heavy tails）を持ち、$q < 1$ では有限の支持域（compact support）を持ちます。

C. エントロピーと双対性（Duality）

本理論の最も重要な発見の一つは、動的指数 $q$ と熱力学的指数 $2-q$ の間の双対性です。

• 動的側面: 局所的な軌道は指数 $q$ によって支配されます。
• 熱力学的側面: 定常状態の安定性や時間発展の不可逆性（H 定理）は、指数 $2-q$ で定義される一般化エントロピー（Tsallis エントロピー $S_{2-q}$）に基づく Lyapunov 汎関数（自由エネルギー）によって支配されます。
• 意義: この双対性により、エスコート分布のような補助的な概念を必要とせず、物理観測量の標準的な期待値を用いた熱力学の定式化が可能になります。

D. 具体的な系への適用

1. 調和振動子: 調和ポテンシャル下で q-ガウス分布が得られ、これが線形復元力に対する唯一の定常状態であることを示しました。
2. 自由粒子: ポテンシャルがない場合、方程式は多孔質媒体方程式に帰着します。自己相似解を仮定することで、平均二乗変位が以下の異常拡散則に従うことを導きました。
   $$\langle x^2(t) \rangle \propto t^{\frac{2}{3-q}}$$
   • $q=1$: 正常拡散 ($\propto t$)
   • $q<1$: 遅い拡散（サブ拡散）
   • $1<q<3$: 速い拡散（スーパー拡散）

4. 結論と意義

本論文は、非線形 Fokker-Planck 方程式を「線形化の原理」に基づいて幾何学的に導出しました。

• 物理的整合性: 非線形拡散を自然に導き出しながらも、ドリフト項の線形性とアインシュタイン関係式を保持し、物理的な直観と整合性を取りました。
• 理論的統一: 動的指数 $q$ と熱力学的指数 $2-q$ の双対性を明らかにし、非拡張統計力学におけるエントロピー指数の選択基準を明確にしました。
• 応用可能性: この幾何学的枠組みは、1 次元量子流体の巨視的状態（共著論文 [7] で詳述）との橋渡しとしても機能し、情報幾何学（Information Geometry）の観点からは、双対平坦多様体上の線形勾配流として解釈できます。

本研究は、補正則的な制約や現象論的な力に頼ることなく、異常拡散とべき乗則分布を記述する一貫した理論的基盤を提供するものです。