[论文解读] Population dynamics in a random environment

Irene Giardina, Jean‐Philippe Bouchaud+1 more|arXiv (Cornell University)|May 11, 2000

Stochastic processes and statistical mechanics被引用 3

一句话总结

本文研究了在随机环境中竞争性亚扩散（由随机力场引起）与超扩散（通过随机增殖引起）共同作用下的种群动力学，揭示了一种新型中间扩散行为。通过在一维空间中的单圈重整化群分析与数值模拟，识别出一种由‘增殖辅助势垒穿越’驱动的新稳定固定点，其异常标度行为可用类Flory论证解释，并发现质心的扩散行为为 $ t^{3/4} $，介于亚扩散与超扩散之间。

ABSTRACT

We investigate the competition between barrier slowing down and proliferation induced superdiffusion in a model of population dynamics in a random force field. Numerical results in $d=1$ suggest that a new intermediate diffusion behaviour appears. We introduce the idea of proliferation assisted barrier crossing and give a Flory like argument to understand qualitatively this non trivial diffusive behaviour. A one loop RG analysis close to the critical dimension d_c=2 confirms that the random force fixed point is unstable and flows towards an uncontrolled strong coupling regime.

研究动机与目标

• 理解随机力场引起的亚扩散与随机增殖引起的超扩散在种群动力学中的相互作用。
• 确定这两种机制竞争时是否会出现新的稳定固定点。
• 利用类Flory论证，对观测到的中间扩散行为提供定性和定量解释。
• 通过数值模拟与临界维数 $ d_c = 2 $ 附近的重整化群分析，验证该新行为的存在。

提出的方法

• 通过具有高斯分布的随机力 $ \vec{F}(\vec{x}) $ 与随机增长率 $ \eta(\vec{x},t) $ 的福克-普朗克型方程建模种群密度 $ P(\vec{x},t) $。
• 将方程转换为自由能描述 $ h(\vec{x},t) = \log P(\vec{x},t) $，得到一个非线性朗之万方程，其中 $ \lambda = \nu_0 $。
• 在 $ d=2-\epsilon $ 下应用单圈重整化群（RG）分析，研究随机力与增殖耦合常数的固定点稳定性。
• 采用类Flory的变分论证，估算有效增殖增益 $ \mathcal{G}(x^*,t) \sim \sigma_\eta \left( \int_0^t dt' / w(t')^d \right)^{1/2} $，其中 $ w(t) $ 为路径宽度。
• 在 $ d=1 $ 中进行数值模拟，测量质心 $ x_{\text{cm}}^2(t) $ 与宽度 $ \Delta^2(t) $，并分析势垒穿越时间。
• 将模型扩展至分数布朗运动势场，其赫斯特指数为 $ H $，并将预测与数值结果进行对比。

实验结果

研究问题

• RQ1 亚扩散（由势垒引起）与超扩散（由增殖引起）的竞争是否会在种群动力学中导致新的稳定固定点？
• RQ2 在同时存在随机力与随机增长的情况下，质心运动 $ x_{\text{cm}}^2(t) $ 的标度行为如何？
• RQ3 类Flory论证能否定量解释观测到的中间扩散指数？
• RQ4 势垒穿越时间如何依赖于势垒高度 $ U_0 $？增殖是否改变了标准的活化（指数）依赖关系？
• RQ5 长程相关无序（通过赫斯特指数 $ H $ 的分数布朗运动实现）对有效扩散指数有何影响？

主要发现

• 在 $ d=1 $ 的数值模拟中发现，质心以 $ x_{\text{cm}}^2(t) \sim t^{3/4} $ 扩散，表明存在一种非平凡的中间扩散行为。
• 单圈重整化群分析表明，当存在增殖时，亚扩散固定点（Sinai）变得不稳定，系统趋向强耦合区域。
• 类Flory论证通过平衡增殖增益 $ \mathcal{G} \sim \sigma_\eta \sqrt{t / w(t)} $ 与势垒高度 $ \sim \sigma_F \sqrt{x^*} $ 解释了 $ t^{3/4} $ 的标度行为，得出 $ x^* \sim t^{3/4} $。
• 有效势垒穿越时间以幂律形式增加 $ \tau \sim U_0^b $，其中 $ b \approx 5/3 $，与纯亚扩散中的指数（活化）行为形成对比。
• 对于赫斯特指数 $ H $ 的分数布朗运动，模型预测 $ \nu^* = 1/(1+2H) $（当 $ H > 1/4 $ 时），数值结果验证了该标度关系，其中 $ \nu^*(H=3/4) \approx 0.37 $，$ \nu^*(H=1/4) \approx 0.65 $。
• 增殖辅助势垒穿越机制通过将标准扩散概率乘以增殖因子 $ \exp(\mathcal{G}) $，增强了穿越高势垒的概率，从而补偿了能量势垒带来的指数抑制效应。