[論文レビュー] Singular limits for models of selection and mutations with heavy-tailed mutation distribution

Sepideh Mirrahimi|arXiv (Cornell University)|Jul 27, 2018

Evolution and Genetic Dynamics参考文献 38被引用数 9

ひとこと要約

本稿は、重い尾を持つ変異を伴う進化的ダイナミクスをモデル化する非局所反応拡散方程式を、分数ラプラシアンと非局所反応項を用いて分析する。特異的摂動スケーリングを適用することで、集団密度が移動するデルタ関数に集中することを示し、ハミルトニアン・ジャコビ法を、従来の粘性解収束が失敗する重い尾を持つ変異核に拡張する。WKB変換された解は最小粘性超解に収束し、単調な増殖率のもとで集中を厳密に導出可能となる。

ABSTRACT

In this article, we perform an asymptotic analysis of a nonlocal reaction-diffusion equation, with a fractional laplacian as the diffusion term and with a nonlocal reaction term. Such equation models the evolutionary dynamics of a phenotypically structured population. We perform a rescaling considering large time and small effect of mutations, but still with algebraic law. We prove that asymptotically the phenotypic distribution density concentrates as a Dirac mass which evolves in time. This work extends an approach based on Hamilton-Jacobi equations with constraint, that has been developed to study models from evolutionary biology, to the case of fat-tailed mutation kernels. However, unlike previous works within this approach, the WKB transformation of the solution does not converge to a viscosity solution of a Hamilton-Jacobi equation but to a viscosity supersolution of such equation which is minimal in a certain class of supersolutions.

研究の動機と目的

• 重い尾を持つ変異を伴う表現型進化をモデル化する非局所反応拡散方程式（分数ラプラシアン拡散と非局所反応項を有する）の漸近的挙動を分析すること。
• 従来の粘性解収束が失敗する、太い尾を持つ変異核を有するモデルへ、制約付きハミルトニアン・ジャコビ法を拡張すること。
• 単調な増殖率条件の下で、表現型密度が時間的に移動するデルタ関数に集中することを確立すること。
• 小変異・大時間領域における特異極限を、非局所性および変異核の重い尾の存在にもかかわらず、厳密に正当化すること。

提案手法

• 小さなパラメータ ε を用いたスケーリングにより、変異ステップを縮小しながら代数的尾の性質を保存する。これにより、ε∂tnε が主要項となるスケーリング済方程式が得られる。
• スケーリング済方程式にホプフ＝コール変換を適用し、密度の時間発展を対数変換された関数に関するハミルトニアン・ジャコビ型方程式に変換する。
• 粘性解理論を用いて変換された解の極限を分析し、古典的粘性解ではなく、粘性超解に収束することを示す。
• 極限が正しい集中ダイナミクスを捉えるために、極限に最小性条件を課す。これにより、極限が古典的粘性解でない場合でも、正しい動的挙動を保証できる。
• 比較原理とバリア法を用いて、極限が必要なハミルトニアン・ジャコビ不等式を満たし、超解のクラス内で最小性を満たすことを証明する。
• 極限関数のゼロレベル集合の分析と、集団総数の連続点における挙動を検討することで、集団密度がデルタ関数に集中することを導出する。

実験結果

リサーチクエスチョン

• RQ1 WKB変換が粘性解に収束しない重い尾を持つ変異核を有するモデルに対して、制約付きハミルトニアン・ジャコビ法を拡張できるか？
• RQ2 変異が太い尾分布に従い、変異効果がスケーリングされた場合、表現型密度の正しい極限挙動は何か？
• RQ3 古典的粘性解収束が成立しないにもかかわらず、集団密度が特異極限においてデルタ関数に集中することをどのように厳密に証明できるか？
• RQ4 増殖率 R(x, I) にどのような条件を課すと、極限におけるデルタ関数が定義された動的ダイナミクスに従うようになるか？
• RQ5 この非局所的・非拡散的設定において、最小粘性超解の概念が、正しい漸近的ダイナミクスを捉えるために古典的粘性解の代わりに使用可能か？

主な発見

• WKB変換された解は、あるクラスの超解の中で最小である粘性超解に収束するが、古典的粘性解には収束しない。
• 増殖率 R(x, I) に単調性条件が成り立つ場合、表現型密度は時間的に移動するデルタ関数に集中し、極限ハミルトニアン・ジャコビ方程式のダイナミクスと整合的である。
• 極限関数 u は弱い粘性解の意味でハミルトニアン・ジャコビ不等式を満たし、制約 (14) を満たす超解のクラスの中で最小であるため、集中分布の一意性が保証される。
• 集団総数 I(t) が極限関数 u(t, x) = 0 である点で連続であり、その点における増殖率 R(x, I(t)) は R(x, I(t)) ≤ 0 を満たす。
• 反例を構成することで、極限関数が (14) の2番目の条件を満たさないことが確認され、これは極限関数が超解ではあるが、古典的解ではないことを裏付ける。
• 弱-*測度位相において、スケーリング済集団密度 nε がデルタ関数に収束することが確立され、特異極限における集中現象が確認される。