Pathways to explosive transitions in interacting contagion dynamics

技术摘要：相互作用传染动力学中爆发性转变的路径

问题陈述
复杂系统通常表现出无法通过研究孤立组件来理解的集体行为。在传染动力学的背景下，理解伴随过程之间的动态相互依赖关系如何改变内在平衡和分叉，目前仍存在研究空白。具体而言，虽然孤立的流行病动力学通常经历二阶（连续）转变，而社会动力学由于复杂的复传染机制通常表现出一阶（爆发性）转变，但这些相互作用产生或抑制爆发性转变的路径尚未得到充分表征。这种爆发性转变的出现会损害系统的稳定性，因为微小的扰动可能导致截然不同的动力学结果。本研究旨在探讨两种传播过程相互作用时驱动此类转变的机制，其范围涵盖了从相互协作到负反馈循环的各种情形。

研究方法
作者提出了一个通用的平均场框架，将两个同时进行的传播过程（分别记为动力学 A 和动力学 B）耦合在一个由 $N$ 个个体组成的封闭种群中。该模型将标准的易感-感染-康复（SIS）分室模型扩展到了四个状态：易感（$S$）、仅感染 A（$A$）、仅感染 B（$B$）以及共同感染（$AB$）。

动力学由一组微分方程描述，用于描述种群中处于各感染状态的比例（$\rho_A, \rho_B, \rho_{AB}$）的演化。该模型包含以下要素：

1. 康复： 两种过程的统一康复率均为 $\mu = 1$。
2. 相互作用率： 感染率受耦合参数 $\epsilon_{AB}$ 和 $\epsilon_{BA}$ 的调节。
   • $\epsilon_{ij} = 1$：非相互作用过程。
   • $\epsilon_{ij} > 1$：协作/兴奋性相互作用（存在 $i$ 会增强 $j$）。
   • $\epsilon_{ij} < 1$：竞争/抑制性相互作用（存在 $i$ 会抑制 $j$）。
3. 传染性质：
   • 动力学 B 被建模为具有恒定传染率 $\lambda_B$ 的标准流行病，表现出二阶转变，在 $\lambda_B = 1$ 时发生。
   • 动力学 A 被建模为具有可调协同传染率 $\tilde{\lambda}_A = \lambda_A e^{\gamma(1-\rho_{AT})}$ 的过程，其中 $\rho_{AT}$ 是受 A 感染的总比例。参数 $\gamma$ 控制转变的性质：$\gamma = 0$ 对应标准流行病动力学，而 $\gamma < 0$ 则引入了具有社会传染特征的群体压力效应，能够诱发爆发性（一阶）转变。当 $\gamma < \gamma_c$ 时，这种效应尤为显著。

研究通过解析和数值方法分析了系统的稳态和分叉点。数值模拟采用控制参数的前向和后向连续化方法，以识别作为一阶转变特征的双稳态和滞后环。

核心贡献与结果

1. 流行病的相互协作 ($\gamma=0, \epsilon > 1$)：
   该模型重现了既有文献中的结果，即两个标准的流行病过程相互增强。这种相互协作将典型的二阶转变转化为一阶（爆发性）转变，从而创造了一个双稳态区域，使系统可以从无病状态突然跳跃到流行状态。相反，相互竞争 ($\epsilon < 1$) 会导致传染性较弱的病原体灭绝，但不会改变转变的连续性质。

2. 单向耦合（社会动力学驱动流行病）：
   当社会传染（动力学 A, $\gamma < 0$）单向驱动流行病（动力学 B, $\epsilon_{BA}=1, \epsilon_{AB} > 1$）时，B 的转变性质取决于社会压力 ($\gamma$) 与耦合强度 ($\epsilon_{AB}$) 之间的相互作用。

   • 如果动力学 A 本身不具备爆发性，则 B 的转变仍保持连续。
   • 如果动力学 A 本身是爆发性的，它可以在耦合强度超过临界阈值 $\epsilon^c_{AB}$ 的条件下，诱发 B 的爆发性转变。
   • 作者推导出了决定驱动流行病发生不连续转变的临界传染率 $\lambda^c_A$ 和临界耦合强度 $\epsilon^c_{AB}$ 的解析表达式。

3. 负反馈循环（抑制爆发性）：
   研究探讨了动力学 A 增强其抑制因子 B（负反馈，$\epsilon_{AB} > 1, \epsilon_{BA} < 1$）的情景。

   • 活跃的抑制因子 B 的存在可以抑制动力学 A 的爆发性。随着耦合强度的增加，与 A 的一阶转变相关的滞后环会缩小并最终消失，从而恢复平滑的（二阶）转变。
   • 这种抑制作用在抑制因子处于活跃状态且耦合强度足够高时最为有效。这揭示了一种提高系统鲁棒性的机制。

4. 带有复传染的相互竞争：
   在涉及两个过程间存在相互抑制相互作用（其中至少有一个过程表现出复传染/协同速率）的情景中，系统仍可能表现出不连续转变。具体而言，具有协同传输率的相互排斥竞争性病原体可以表现出爆发性转变，这一现象不同于标准竞争流行病中观察到的平滑灭绝现象。

意义
本文为相互作用传染过程中爆发性现象的路径提供了一个通用的视角。它表明，爆发性不仅仅是单个过程的内在属性，还可以通过同时进行的动力学之间相互依赖的性质而产生、传递或被抑制。

作者声称，其框架统一了以往关于相互协作和单向驱动的研究发现，同时揭示了一条新的稳定性路径：负反馈循环可以有效地“关闭”一阶转变，从而增强系统的鲁棒性。这一发现与基因调节网络中负反馈确保稳定性的观察结果相一致。该工作强调，在对同时传播的传染过程进行建模时，必须考虑相互依赖关系，因为这些相互作用会从根本上改变动力学状态和转变性质，使其不同于孤立系统。该研究局限于混合均匀的人口模型，作者指出，在复杂接触网络上探索这些动力学仍是未来的工作方向。