[论文解读] 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.

研究动机与目标

• 分析带有分数阶拉普拉斯扩散和非局部反应的非局部反应-扩散方程的渐近行为，以建模重尾突变下的表型进化。
• 将带约束的哈密顿-雅可比方法扩展至具有重尾突变核的模型，其中先前的粘性解收敛性失效。
• 在单调生长率条件下，建立表型密度随时间移动的狄拉克质量集中。
• 尽管突变核具有非局部性和重尾特性，仍严格证明小突变、大时间 regime 下的奇异极限。

提出的方法

• 使用小参数 ε 对模型进行标度，缩小突变步长同时保持代数尾部行为，得到以 ε∂tnε 为领先项的标度方程。
• 对标度方程应用霍普夫-科尔变换，将密度演化转化为对数变换后的哈密顿-雅可比型方程。
• 利用粘性解理论分析变换后解的极限，表明其收敛至粘性上解而非经典粘性解。
• 对极限施加最小性条件，以确保其捕捉正确的集中动力学，即使极限不是经典粘性解。
• 采用比较原理和障碍论证，证明极限满足所需哈密顿-雅可比不等式，并在上解类中具有最小性。
• 通过分析极限函数的零水平集及其在总种群大小连续点处的行为，推导出种群密度的狄拉克质量集中。

实验结果

研究问题

• RQ1 能否将带约束的哈密顿-雅可比方法扩展至重尾突变核模型，其中WKB变换不收敛至粘性解？
• RQ2 当突变服从重尾分布且突变效应被缩小后，表型密度的正确极限行为是什么？
• RQ3 在缺乏经典粘性解收敛性的情况下，如何严格证明种群密度在奇异极限下集中为狄拉克质量？
• RQ4 对生长率 R(x, I) 的何种条件可确保极限狄拉克质量按明确定义的动力学演化？
• RQ5 在非局部、非扩散设定下，最小粘性上解概念是否可替代经典粘性解以捕捉正确的渐近动力学？

主要发现

• WKB变换后的解收敛至哈密顿-雅可比方程的一个粘性上解，且在某一类上解中为最小解，而非经典粘性解。
• 在生长率 R(x, I) 的单调性条件下，表型密度集中为随时间移动的狄拉克质量，与极限哈密顿-雅可比方程的动力学一致。
• 极限函数 u 以弱粘性意义满足哈密顿-雅可比不等式，且在满足约束 (14) 的上解类中具有最小性，确保了集中轮廓的唯一性。
• 在极限函数 u(t, x) = 0 的点，总种群大小 I(t) 连续，且在这些点处的生长率 R(x, I(t)) 满足 R(x, I(t)) ≤ 0。
• 构造了一个反例，表明极限函数可能不满足 (14) 的第二条条件，从而确认其仅为上解而非经典解。
• 在测度的弱-* 拓扑下，建立了标度种群密度 nε 向狄拉克质量的收敛性，证实了奇异极限下的集中现象。