凝胶化作为Smoluchowski流的动力学不稳定性

《The Journal of Physical Chemistry B》：Gelation as a Dynamical Instability of the Smoluchowski Flow

　　

研究人员将凝胶化定义为Smoluchowski凝聚流失去局部动力学稳定性的现象，并通过沿演化聚集动力学的雅可比矩阵（Jacobian）的时间依赖谱进行量化。针对形式为K(i,j) = (ij)α的齐次核及经典Smoluchowski核，研究表明凝胶化始终先于正实特征值的出现，表明局部动力学稳定性丧失；非凝胶核仅表现出瞬态有限尺寸效应，而凝胶核则伴随宏观凝胶形成的持续谱失稳。这些结果将凝胶化纳入聚集驱动相变的统一动力学框架，通过谱特征将凝胶化的动力学逃逸与Becker–D?ring动力学中成核的亚稳态联系起来。

研究背景方面，凝胶化在传统理论中主要通过Smoluchowski凝聚方程的矩发散或质量损失来表征，但这种静态描述难以揭示其集体机制。现代研究倾向于将其视为非平衡连续相变，由逾渗过程控制，并伴随临界指数和幂律分布。然而，数值模拟中有限尺寸效应与瞬态质量损失使得区分真实凝胶化和非凝胶核的行为变得困难。因此，研究人员从动力学系统角度提出了新的稳定性判据，旨在将凝胶化解释为Smoluchowski流本身的动力学不稳定性，而非仅依赖矩的奇异行为。这项工作发表在《The Journal of Physical Chemistry B》。

关键技术方法方面，研究人员采用截断离散Smoluchowski方程，引入最大簇尺寸截断p，构建有限维非线性常微分方程系统；选取两类核函数：齐次乘法核K(i,j) = (ij)^α（α > 1/2为超临界凝胶核，α ≤ 1/2为亚临界非凝胶核）及经典Smoluchowski扩散限制核（D_f= 1.8）；沿轨迹进行时间依赖线性化，计算雅可比矩阵及其特征值谱；通过截止尺寸缩放区分物理凝胶化与数值伪影；结合全局观测指标，如矩分析（M₁、M₂）与香农熵H(t)，追踪聚集状态与分布复杂性。

结果与讨论部分，首先在模型与稳定性框架中，研究人员建立了离散Smoluchowski方程并推导了雅可比矩阵的解析形式，指出由于无稳态解，稳定性分析需沿轨迹进行时变谱评估。其次，在核函数选择上，齐次乘法核的超临界与亚临界划分与经典理论一致，扩散限制核作为非凝胶对照。第三，有限维截断与收敛性分析表明，当截断尺寸足够大时，谱失稳是无限维流动的真实特征。第四，线性化方法揭示了正特征值对应局部扰动的指数增长，是动力学失稳的直接信号。第五，谱指标显示，超临界核在凝胶点前出现持续的正特征值，亚临界核则仅存在瞬态微扰，扩散限制核始终保持谱稳定。第六，数值设计与实现中，ODE求解器实时计算雅可比谱，并与矩分析和香农熵结合，形成完整的动力学诊断体系。

研究结果表明，在α = 0.2的非凝胶体系中，M₁守恒，M₂单调增长，雅可比最大实部Re(λ_max) ≈ 0，香农熵持续增长，簇尺寸分布呈指数衰减；而在α = 0.8的凝胶体系中，M₂在凝胶点发散，M₁下降，Re(λ_max)在凝胶点前跃升至峰值后转为负值，香农熵在凝胶后下降，簇尺寸分布形成幂律尾。扩散限制核始终无正特征值，验证了谱方法的特异性。

讨论与结论部分，研究人员指出凝胶化应理解为Smoluchowski流的局部动力学失稳，而非单纯的矩奇异性。这一发现连接了可逆Becker–D?ring方程的成核势垒跨越与不可逆Smoluchowski方程的动力学逃逸，表明这两类过程在谱拓扑上具有共性。谱失稳作为相变的普适标志，可推广至远离平衡的多种化学演化系统，为预测动力学逃逸提供了稳健工具。