Lie symmetry analysis and similarity reductions for the tempered-fractional Keller Segel system

技術要約：テンパード分数型ケラー・セーゲル系におけるリー対称解析と相似縮退

問題提起
本論文は、異常拡散を伴う化学走性を記述するテンパード分数型ケラー・セーゲル（TFKS）モデルの数学的解析に取り組んでいる。古典的なケラー・セーゲル（KS）モデルは通常の拡散を仮定し、しばしば有限時間での爆発（特異点形成）を示すが、分数型KS（FKS）モデルは長距離ジャンプを捉えるものの、依然として爆発の問題を抱える可能性がある。これに対し、TFKSモデルは指数関数的なテンパリング（緩和）機構を導入している。この機構は、極端な長距離ジャンプを遮断しつつ、メモリ効果を保持するものである。本研究が取り組む具体的な課題は、この系へのリー対称解析の適用である。標準的なリー対称手法は局所的な延長公式に依存しているが、これはテンパード分数型微分（時間および空間の両方）に固有の非局所性や、標準的な時間並進不変性を破る指数関数的なテンパリング因子と相容れない。

手法
非局所性とテンパリングの障壁を克服するために、著者らは二段階の縮退手法を採用している。

1. ゲージ変換: 変換 $v(x,t) = e^{\lambda t}u(x,t)$ および $w(x,t) = e^{\lambda t}c(x,t)$ を適用する。恒等式 $D^{\alpha,\lambda}_t u = e^{-\lambda t} {}^C D^\alpha_t (e^{\lambda t} u)$ を利用することで、テンパード・カプト（Caputo）時間微分を標準のカプト微分へと変換する。これにより、時間演算子に関する既存の分数型リー対称性の枠組みとの適合性が回復される。
2. 短距離ジャンプ拡散近似: 非局所的なテンパード分数型ラプラシアン $L^{(\lambda)}_\alpha$ を、局所的な二階拡散演算子 $c_2(\alpha, \lambda)\Delta$ によって近似する。この近似は、基礎となるレヴィ過程の特性ジャンプ長がマクロな空間スケールに対して小さいという仮定の下で有効である。このステップにより、元の積分微分系を局所的な分数偏微分方程式（PDE）へと変換し、古典的なリー点対称技術の直接的な適用を可能にする。

主要な貢献と結果
本論文は、一次元空間（$0 < \alpha < 1$）における縮退された局所TFKS系に対して完全なリー対称解析を行っている。

• 対称性の分類:

  • 一般的ケース ($\lambda > 0, \chi, r, \kappa > 0$): テンパリング因子 $e^{-\lambda t}$ によって導入された明示的な時間依存性と非線形結合項により、対称群は著しく制限される。許容される連続対称性は、空間並進（$X_1 = \partial_x$）のみである。
  • 条件付き拡張: 特定のパラメータ制約下でのみ、追加の対称性が現れる：
    • $\lambda = 0$（テンパリングなし）の場合、系は時間並進（$\partial_t$）と空間並進を認め、アーベル群を形成する。
    • 化学走性（$\chi$）およびロジスティック成長（$r$）が欠如している場合、修正されたスケーリング・並進対称性（$\partial_t - \frac{\lambda}{2}x\partial_x$）が現れ、局所近似におけるテンパリング効果を補償する。
  • 最適系: 著者らは、一般的ケース、無テンパードケース（$\lambda=0$）、および退化線形ケース（$\chi=r=0$）について、一次元部分代数の最適系を構築する。

• 相似縮退と厳密解:
  最適系を用いて、縮退されたPDEは常微分方程式（ODE）または分数常微分方程式へと変換される。

  • 空間的に均質な解: 一般的ケースにおいて、系は分数常微分方程式系へと縮退する。分数複素変換（FCT）を用いることで、求積法による厳密な解析解が導出される。
    • 無テンパード ($\lambda=0$): 解は $t \to \infty$ において環境収容力 $K_0$ に収束する。
    • テンパード ($\lambda>0$): 条件 $0 < v_0 < K_0$ の下で、解は一時的な上昇を示した後、ゼロへの単調減少（絶滅）を示す。テンパリング因子 $e^{-\lambda t}$ は、ロジスティック成長を最終的に圧倒する実効的な時間依存的損失項として機能する。
  • 線形退化ケース ($\chi=r=0$): 空間的に一様かつスケーリング不変なシナリオの両方について、厳密解が見出される。スケーリング解は、分数的な操作時間 $s(t) = t^\alpha/\Gamma(1+\alpha)$ と全体的な指数減衰因子を持つ、分数型テンパード熱核（空間的にはガウス分布）を与える。
  • 定常状態: 解析により、定常平衡 $(0,0)$ および $(K_0, K_0/\kappa)$、ならびに任意の細胞密度 $u$ が与えられた場合の化学誘引物質濃度 $c$ のグリーン関数表現が特定される。

• 生物学的ダイナミクス:
  導出された解は、テンパリングパラメータ $\lambda$ とロジスティック成長項の決定的な役割を浮き彫りにしている。ロジスティック成長項は、細胞増殖を制限することで個体群ダイナミクスを制御し、無制限な集積傾向を抑制し、適切な条件下では有限時間の爆発を防ぐ。対照的に、テンパリング機構は漸近的挙動を変化させる。テンパリングがない場合、モデルは環境収容力での生存を予測するが、テンパリングがある場合、初期密度が $0 < v_0 < K_0$ を満たすとき、指数減衰の累積効果により、モデルは最終的な絶滅を予測する。分数次数 $\alpha$ は、拡散の速度と集積の形状に影響を与え、より低い $\alpha$ 値は、劣拡散的でメモリ効果が支配的な運動に対応する。

意義と主張
本論文は、TFKS系にリー対称群を適用した最初の研究であると主張している。主な意義は、古典的な対称解析が可能な形式へとテンパード分数型演算子を変換することにより、厳密な数学的枠組みを提供したことにある。結果は、テンパード異常拡散を伴う化学走性系の長期挙動に関する解析的な洞察を提供しており、特に、テンパリング機構が（$0 < v_0 < K_0$ の条件下で）集積ではなく絶滅へと個体群を向かわせること、一方でロジスティック項が有限時間の爆発を防ぐ役割を果たすことを示している。この手法は他のテンパード分数型微分方程式にも適用可能であり、古典的なモデリング手法と分数型モデリング手法の架け橋となるものである。著者らは、短距離ジャンプ近似によってこれらの解析的結果が可能になったものの、変数とテンパリングパラメータ間の自然なスケーリング関係を不明瞭にする可能性があること、また現在の解析は劣拡散領域（$0 < \alpha < 1$）に限定されていることを注記している。