[論文レビュー] On a kinetic description of Lotka-Volterra dynamics

Giuseppe Toscani, Mattia Zanella|arXiv (Cornell University)|Feb 28, 2023

Mathematical and Theoretical Epidemiology and Ecology ModelsMedicine被引用数 3

ひとこと要約

本稿では、出生プロセスを記述する再分配作用素を備えたボルツマン型方程式を用いて、ロトカ＝ヴォリエラの捕食者・被捕食者ダイナミクスの運動論的モデルを導入する。平均個体数は古典的ロトカ＝ヴォリエラ方程式に従うが、群のサイズの完全な確率分布は再分配機構に強く依存し、特定の政策下では被捕食者において肥尾分布が生じる。これは平均場的挙動とは顕著に異なる統計的挙動を示している。

ABSTRACT

Owing to the analogies between the problem of wealth redistribution with taxation in a multi-agent society, we introduce and discuss a kinetic model describing the statistical distributions in time of the sizes of groups of biological systems with prey-predator dynamic. While the evolution of the mean values is shown to be driven by a classical Lotka-Volterra system of differential equations, it is shown that the time evolution of the probability distributions of the size of groups of the two interacting species is heavily dependent both on a kinetic redistribution operator and the degree of randomness present in the system. Numerical experiments are given to clarify the time-behavior of the distributions of groups of the species.

研究の動機と目的

• 群のサイズの統計的分布（平均値のみでなく）を捉える、捕食者・被捕食者系の運動論的記述を構築すること。
• 出生イベントにおける異なる再分配政策が、個体数分布の形状に与える影響を調査すること。
• 租税に類似した再分配機構を通じて、富の交換の運動論的モデルと生物学的捕食者・被捕食者ダイナミクスを結びつけること。
• 特定の再分配ルール下で、被捕食者個体群に非ガウス的で肥尾な分布がどのように出現するかを分析すること。

提案手法

• 被捕食者および捕食者群のサイズの同時確率密度関数を記述する非マクスウェル型ボルツマン型方程式系を定式化する。
• 出生をモデル化するために、分布の一次および二次モーメントにのみ依存する運動論的再分配作用素を導入する。
• すり減らし極限（grazing limit）を適用し、小さな相互作用閾値に対してボルツマンダイナミクスを近似するフォッカー＝プランク系を導出する。
• 空間グリッド上でのフォッカー＝プランク系の数値解法として、4次スペクトル法と4次ルンゲ＝クッタ時間積分法を用いる。
• 再分配作用素のパラメータχを変化させ、比例的および反比例的出生政策を模擬する。
• 数値シミュレーションを実施し、平均ダイナミクスが古典的ロトカ＝ヴォリエラ系と一致するかを検証するとともに、分布全体の時間発展を分析する。

実験結果

リサーチクエスチョン

• RQ1 運動論的再分配作用素は、ロトカ＝ヴォリエラ系における被捕食者および捕食者群のサイズの統計的分布の形状にどのように影響を与えるか？
• RQ2 平均ダイナミクスが古典的ロトカ＝ヴォリエラ系に回帰する一方で、完全な分布が非自明な統計的特徴を示す条件は何か？
• RQ3 どのような条件下で、被捕食者個体群の分布に肥尾が出現し、再分配政策はこの現象にどのような役割を果たすか？
• RQ4 フォッカー＝プランク近似は、元のボルツマン型運動論的系の挙動をどの程度適切に捉えているか？
• RQ5 フォッカー＝プランク系の準定常解と、個体群分布の長期的安定性との関係は何か？

主な発見

• 被捕食者および捕食者群のサイズの平均値は、古典的ロトカ＝ヴォリエラ系の常微分方程式系に従い、既知のダイナミクスと整合していることが確認された。
• 群サイズの完全な確率分布は、再分配作用素の形式、特にパラメータχに極めて敏感であることが判明した。
• χ = -1（高サイズ群に新規出生を集約する）の場合、時間の経過とともに被捕食者分布に肥尾が出現し、t = 15における対数-対数プロットで明確に観察された。
• 被捕食者分布に肥尾が出現するのは再分配機構に起因するものであり、比例的または一様な出生政策ではこの現象は発生しない。
• N = 801グリッドポイントを用いた4次スぺクトル法による数値シミュレーションでは、運動論的モデルの平均ダイナミクスと古典的ロトカ＝ヴォリエラ系との間に良好な一致が得られた。
• フォッカー＝プランク近似は、分布の特徴（特に尾部の挙動）を詳細に分析可能にするよう、システムの時間発展を効果的に捉えていた。