协同繁殖晶格Brownian bees模型的稳态、有限时间坍缩与扩散式铺展：流体动力学自由边界问题研究

《Physica D: Nonlinear Phenomena》：Lattice Brownian bees with cooperative reproduction: Steady states, collapse, and spreading

　　

研究人员将Berestycki等（2021，2022）的“Brownian bees（布朗蜂）”模型扩展至协同繁殖，kA→(k+1)A：由N个对称随机游走体组成的种群在每次出生事件中移除距原点最远的粒子。在N→∞极限下，研究人员为该模型建立流体动力学（hydrodynamic，HD）自由边界问题。利用该形式，研究人员确定所有k的稳态种群密度，并证明k≤2时线性稳定、k≥4时不稳定。在临界情形k=3，仅在繁殖率与扩散率的一个临界比值处存在整个连续稳态族。高于临界性时，种群经历向原点的渐近自相似有限时间坍缩。低于临界性时，种群呈扩散式铺展，但繁殖仍具有定量相关性。对于k≥4，不稳定稳态将有限时间坍缩与扩散式铺展两类区域分隔开。此处坍缩动力学为渐近自相似，种群密度呈现尺度分离，需要匹配渐近描述。研究人员的解析预测由HD自由边界问题的数值解和原始微观模型的Monte Carlo（MC）模拟所证实。

一、研究背景、问题与意义
N粒子分支Brownian运动选择模型是一类重要的随机系统：粒子扩散、分支，并按选择规则移除以保持种群规模固定。不同移除规则导致不同宏观行为，但在N→∞的流体动力学（HD）极限下通常归结为自由边界反应–扩散问题。既往研究多为简单分裂A→2A；Brownian bees模型在每次分支时移除距原点最远粒子，k=1时达到唯一球对称稳态。协同繁殖kA→(k+1)A在生物中见于群体感应介导繁殖、多配子受精和合作繁殖等。HD描述中，k体繁殖引入局部密度非线性项λu^k，k>1时从根本上改变空间结构并产生新的动力学区域。因此，研究人员将一维晶格Brownian bees推广至协同繁殖，系统研究所有k的稳态、线性稳定性、坍缩与铺展，以揭示高阶繁殖对空间选择种群宏观行为的影响。

二、研究工作、结论与意义
研究人员保留“最远粒子移除”规则，构建N?1一维晶格对称随机游走协同繁殖微观模型，取N→∞和连续极限得到HD自由边界反应–扩散方程及质量约束。研究得到所有k的精确稳态密度，发现k=3处发生线性稳定性转变：k≤2稳定，k≥4不稳定，k=3边缘且仅在临界α_c=λ/D处存在连续稳态族。k=3临界以上呈自相似有限时间坍缩，临界以下呈自相似扩散铺展但繁殖仍定量相关；k≥4不稳定稳态分隔坍缩与扩散，k=4坍缩具尺度分离并需匹配渐近描述。解析预测与HD数值解及原始微观模型Monte Carlo（MC）模拟一致。其意义在于：空间选择最小模型已强烈偏好低阶繁殖，k≤2可形成稳定宏观集体，k≥4无稳定配置，k=3位于边界，呼应生物中二分繁殖常见而高阶协作繁殖稀少。

三、关键技术方法（≤250字）
本研究为理论建模与数值模拟研究，无样本队列来源。方法包括：建立一维晶格N?1对称随机游走协同繁殖微观模型，取N→∞与连续极限得HD自由边界反应–扩散方程；用量纲分析确定k=3为质量临界指数及临界比α_c=λ/D；将稳态化为Lane-Emden方程，借助Euler beta函数、正则不完全beta函数、Jacobi椭圆余弦、Weierstrass椭圆函数和超几何函数获得精确剖面；将线性稳定性化为Schr?dinger型特征值问题，结合Sturm比较定理和射击法分析；对k=3、4采用自相似ansatz和匹配渐近展开；以HD自由边界数值解和原始微观模型MC模拟验证。

四、研究结果
3.1 稳态（Steady states）：通过稳态Lane-Emden方程、边界条件和归一化质量约束，研究人员以第一积分和Euler beta函数求出半宽L与峰值密度u_m。k<3时质量约束唯一确定u_m，稳态唯一；k=3时质量约束固定λ/D而非u_m，仅在α_c存在连续稳态族；k>3时稳态唯一但不稳定。k=1退化为已知余弦剖面。

3.2 线性稳定性（Linear stability）：通过小扰动u=u_s+εv e^σt线性化，得到Schr?dinger型特征值问题，边界条件要求v(L)=v(?L)，一阶质量守恒要求∫vdx=0。奇模由Sturm比较定理排除为不稳定来源；偶模用射击法求主特征值，得k<3稳定、k>3不稳定；k>3的不稳定性还可用质量斜率判据解析证明。

3.3 k=2：稳定稳态（Stable steady state）：通过无参数化方程?″+?²=0及边界条件，得到?_m≈0.224、ρ≈3.63，密度剖面亦可用Weierstrass椭圆函数表示。HD剖面与MC模拟吻合，质心位置?x?在HD中为零而有限N下涨落，验证HD极限及k=2稳态稳定性。

3.4 k=3：边缘稳定性（Marginal stability）：由质量约束得α_c=π²/2≈4.935；在临界α处存在以u_m参数化的连续稳态族，半宽满足Lu_m=K(1/2)/α_c，剖面为Jacobi椭圆余弦。α≠α_c时无符合质量约束稳态，动力学由自相似形式u=L^?1F(x/L)描述。3.4.1 自相似坍缩（α>α_c）：通过非线性特征值问题F″+αF³=γ(F+yF′)，得L～[2Dγ(T?t)]^1/2、u(0,t)～(T?t)^?1/2，HD与MC证实。3.4.2 自相似扩张（α<α_c）：通过F″+αF³+β(F+yF′)=0得L～(2Dβt)^1/2，铺展非高斯且繁殖保持定量相关。

3.5 k≥4：不稳定稳态、坍缩和铺展（Unstable steady state, collapse and spreading）：通过稳定性分析与HD/MC动力学，发现不稳定稳态在函数空间中起分界线作用：较窄初条件坍缩，较宽初条件扩散铺展；k=4示例中稳态半宽L_SS≈15.7。3.5.1 k=4坍缩的匹配渐近描述：体内繁殖主导、边缘窄扩散边界层主导，边界层厚度δx～L²；外部解含形状参数r，μ₀(r)=λ[2I(r)]^?3，得L～[3μ₀(T?t)]^1/3、u_max～τ^?1/3；边界层解H₀(z)=F_b(1?e^?μ₀z/D)，r由初条件决定。3.5.2 扩散铺展：宽初条件时λu⁴在体内渐次次要，体密度趋于高斯，方差σ²～2Dt，边界位置L(t)～(2Dt ln t)^1/2，MC与HD数值解证实。

五、讨论与结论
讨论部分指出，该工作将Brownian bees扩展至协同繁殖并建立HD模型，核心结果是k=3处线性稳定性转变：k<3稳定、k>3不稳定。k=3仅在α_c=π²/2处存在连续稳态族；低于临界值种群自相似扩散铺展且繁殖始终重要，高于临界值HD意义下发生自相似有限时间坍缩，坍缩率和扩张率由非线性特征值问题确定。k≥4时不稳定稳态分隔有限时间自相似坍缩与体内渐近高斯的扩散铺展；k=4坍缩表现为繁殖主导的体区与扩散主导的窄边界层共存，最终坍缩被连续HD理论未计入的晶格效应阻止；k=5经MC检查亦符合。解析预测由HD数值解和MC模拟定量证实。该研究与最左粒子移除规则模型比较：本模型无行波，但两者不稳定性阈值同为k=3，反映一维质量临界指数。作者Ohad Vilk与Baruch Meerson认为，空间选择最小模型偏好低阶繁殖：k≤2有稳定宏观集体，k≥4无稳定配置，k=3处边界；这与生物中分裂和二分繁殖常见、高阶合作繁殖稀少相呼应。