本研究由 S. Yanchuk、S. Wieczorek、H. Jardón-Kojakhmetov 和 H. Alkhayuon 共同完成，其中 S. Yanchuk 与 H. Alkhayuon 为共同第一作者。研究机构包括爱尔兰科克大学数学科学学院、德国波茨坦气候影响研究所以及荷兰格罗宁根大学约翰·伯努利数学与计算机科学研究所。该论文于 2026 年 1 月 30 日投稿，2026 年 7 月 24 日被接收，并于 2026 年 9 月 28 日发表在 Physical Review Letters 第 137 卷第 147202 期，同时被编辑推荐并作为亮点报道。
本文属于类型 a，即单项原创研究的学术报告。
一、研究背景与目的
在真实世界的复杂系统中，多重稳定性（multistability）普遍存在，即同一系统从不同初始条件出发，可能最终稳定到完全不同的渐近状态。这种多稳态现象广泛出现在癫痫与神经元模型、激光器、储备池计算以及气候系统等众多领域。在多稳态系统中，系统在受到扰动后能否回到某一给定的稳定状态，取决于该稳定状态的吸引域（basin of attraction）的形状和范围。吸引域的几何结构直接决定了系统对扰动的响应能力、不确定性与韧性（resilience）。
与此同时，现实世界系统往往在多个时间尺度上演化。多尺度系统（multiscale systems），特别是慢-快系统（slow-fast systems），会涌现出多种复杂的非线性现象。在分析这类系统时，研究者通常利用大时间尺度分离极限下的降维方法，例如准静态近似（quasistatic approximation）、绝热消除（adiabatic elimination）和时间平均（time averaging）等。这些方法的基本思路是：要么把快变量视为瞬间达到平衡而将其消去，要么把慢变量冻结以便集中分析快变量层的动力学。
然而，研究多尺度系统中吸引域的几何结构这一领域在很大程度上仍未被探索。本研究的目标正是在于：揭示慢-快系统的吸引域可能具有的特殊几何性质，以及这些性质如何影响系统对扰动的鲁棒性，并检验常用降维方法在处理此类系统时的适用性。
二、研究步骤与方法
本研究通过三类具有代表性的慢-快系统，详细展示了奇异吸引域与奇异漏斗的形成机制、几何结构及其体积尺度关系。
首先，作者研究了带慢适应分岔参数的叉形分岔（pitchfork bifurcation）范式系统。该系统由一个快变量 x 和一个慢变量 μ 组成，快方程描述 x 随 μ 与 x 的三次非线性关系的变化，慢方程描述 μ 随 x 的线性适应过程。当时间尺度比参数 ε 趋于零时，层系统（layer system）在每一个固定 μ 下都具有唯一稳定的准静态平衡点，且该平衡点随 μ 连续变化并吸引所有初始条件。因此，按照传统观点，对该系统进行绝热消除是完全合理的。作者将快变量 x 替换为其准静态平衡值，从而得到仅含慢变量 μ 的约化系统，并构造了相应的双势阱势函数。约化系统具有与全系统相同的平衡点结构，且其两个稳定平衡点的吸引域由不稳定平衡点所对应的单一 μ 值清晰地分隔开。
其次，作者研究了自适应相位转子（adaptive phase rotator）系统。该系统由一个快相位变量 φ 和一个慢适应变量 μ 组成。在层系统层面，当 |ω+μ|<1 时，系统具有唯一吸引的准静态平衡点；当 |ω+μ|>1 时，系统呈现周期旋转。作者对该系统进行了混合降维处理：在 |ω+μ|<1 的区域进行绝热消除，在 |ω+μ|>1 的区域对快旋转进行平均，从而得到约化系统的函数 f(μ)，并构造了双势阱势函数。约化系统同样是双稳态的。
第三，作者构建了平均场耦合的自适应相位转子网络系统，作为高维系统的代表。该系统由 N 个相位转子组成，它们的频率受到全局慢变量 μ 的自适应调控，而 μ 又受到系统平均场 x 的驱动。作者分别研究了 N=2 和 N=10 的情形。对于 N=2 的三维系统，作者通过固定某一相变量或固定 μ 值来获得二维截面，以观察奇异漏斗的拓扑结构。对于 N=10 的十一维系统，由于无法直接可视化高维截面，作者采用蒙特卡洛模拟来估计奇异漏斗的体积随 ε 的变化。
在所有上述系统中，奇异漏斗的边界均由鞍点平衡点 e2 的稳定不变流形（stable invariant manifold）构成。为了揭示奇异漏斗的形成机制，作者追踪了这些边界轨线随时间的逆向演化，发现它们会经过临界流形（critical manifold）的不稳定分支。在逆向时间中，该不稳定分支成为吸引的，导致两条边界轨线在向临界流形靠近时以指数速率相互逼近，从而形成指数窄的漏斗结构。基于这一几何图像，作者从理论上推导出了奇异漏斗体积的普适尺度关系：V(ε)∼exp(−cε−1)，其中 c 为正常数。此外，作者还通过蒙特卡洛模拟数值验证了这一尺度关系。
三、主要研究发现与结果
研究揭示了一个核心发现：在慢-快系统中，吸引域可能包含被称为“奇异漏斗”（singular funnels）的窄隧道状结构。这些漏斗虽然极其狭窄，却能延伸到相空间的不同区域，从而深刻影响系统对扰动的响应能力。
在自适应叉形分岔范式系统中，约化系统（即绝热消除后的系统）的两个吸引域由 μ=μb 这一单一分界值完全分开。然而，在全系统中，无论 ε 多么小（只要非零），这种清晰的吸引域分隔都不成立。奇异漏斗使得初始条件可以从 μb 的任意一侧收敛到左稳定平衡点 e1，这在约化系统中是绝对不可能发生的。从约化系统的视角看，这种由快变量方向的扰动引发的从 e3 到 e1 的转变，如同量子隧穿（quantum tunneling）或“瞬移”（teleporting）一般。这意味着即使层系统在每一个冻结慢变量值下都具有唯一稳定且连续变化的准静态吸引子，绝热消除仍然可能完全失效。
在自适应相位转子系统中，奇异漏斗展现出更为复杂的几何结构。与范式系统中漏斗在 x 方向上呈指数窄小不同，相位转子系统的奇异漏斗在 φ 方向上覆盖了整个相位变量范围，并且其指数窄化仅发生在有限的 μ 区间内。这导致该系统的奇异漏斗体积显著大于范式系统。奇异漏斗延伸到 μ 为负的区域，并形成一条通道，使得具有任意 μ 值但特定 φ 值的初始条件可以被周期旋转吸引子所吸引。这一行为在约化系统中是不可能的，因为约化系统要求初始条件必须满足 μ>μb 才能收敛到周期旋转。此外，奇异漏斗在穿过不稳定临界流形后还会发生多次旋转，以复杂的方式反复与另一个稳定态的吸引域相交。
对于奇异漏斗的体积尺度关系，作者发现低维系统（如自适应叉形分岔系统和自适应相位转子系统）均很好地遵循 V(ε)∼exp(−cε−1) 的普适标度定律。蒙特卡洛模拟结果与理论预测在数值上高度吻合。然而，在两转子高维系统中，当参数 δω 取某些特定值时，奇异漏斗的体积会偏离这一简单标度律，呈现出类似共振的偏差。这表明在更高维系统中，奇异漏斗可能具有更丰富的动力学行为，其体积标度可能受到系统内部共振或其他复杂机制的影响，这被作者列为未来研究的开放问题。
对于十维转子网络系统，蒙特卡洛模拟显示奇异漏斗体积仍随 ε 的减小而减小，这与奇异漏斗的存在一致。同时，两条从极相近初始条件出发的轨线——一条位于奇异漏斗内，另一条位于其外——最终收敛到完全不同的吸引子，直观地证明了奇异漏斗对系统长期行为的决定性影响。
四、研究结论与价值
本研究的核心结论是：奇异漏斗在双稳态慢-快系统中是普遍存在且具有鲁棒性的。只要存在鞍点平衡点的稳定不变流形与多个共存吸引子，奇异漏斗就可能出现，且其对小参数扰动具有结构稳定性。作者通过自适应叉形分岔范式系统、自适应相位转子系统和自适应网络系统三类截然不同的模型，证明了奇异漏斗的普适性。
这一发现具有深远的科学价值。首先，它从根本上挑战了多尺度系统分析中长期被广泛使用的降维方法。即便是准静态近似、绝热消除和时间平均这些成熟技术，在奇异漏斗存在时也可能得出严重误导性的结论。准静态近似可能完全掩盖多个吸引域的存在，而绝热消除则可能使原本在全系统中可达的稳定态变得不可达。因此，当研究者利用约化系统来推断完整系统的韧性或吸引域性质时，必须格外谨慎。
其次，本研究提供了重要的应用价值。真实世界的数学模型不可避免地会忽略某些快变量。本研究的结果意味着：即使模型未能捕捉到某些快变量，由这些被忽略的快变量所诱导的奇异漏斗仍可能使真实系统发生模型无法预测的临界转变。这一见解对气候系统、生态系统和神经科学等领域的模型构建与预测具有警示意义。
五、研究亮点
本研究的亮点体现在三个方面。第一，首次识别并系统刻画了奇异吸引域与奇异漏斗这一全新的几何现象，为理解多尺度系统中的全局动力学提供了新视角。第二，作者构造了奇异漏斗的低维范式系统，并推导出了体积随 ε 指数缩小的普适标度律，为后续理论研究奠定了基础。第三，本研究将奇异漏斗从低维扩展到高维自适应网络系统，并通过精心设计的数值模拟揭示了高维情形下可能出现的标度律偏差，开辟了新的研究问题。该工作因其重要性与新颖性被 Physical Review Letters 选为编辑推荐并作为亮点报道。