Prigogine & Lefever,1968 — 耗散结构(dissipative structure)的最小玩具模型
化学课上教的是:反应向着平衡走,浓度单调地趋近某个值,到了就不再变。可 1950 年代起,实验室里陆续出现了不肯这样收场的溶液——最有名的是别洛乌索夫-扎鲍廷斯基反应(Belousov–Zhabotinsky reaction,通称 BZ 反应),一锅溶液的颜色会周期性地来回变,像个钟一样走上几十分钟。这在当时被普遍认为不可能,别洛乌索夫的论文两度被拒。
于是问题变成:一个只由普通化学反应组成的体系,凭什么会自己振荡起来?需要多少条件才够?普里高津(Ilya Prigogine)和勒菲弗(René Lefever)1968 年在布鲁塞尔给出的回答,就是把条件削到不能再削的一个假想反应网——后来被 John Tyson 戏称为"布鲁塞尔器"(Brusselator,仿照当时的另一个模型 Oregonator 的构词)。
它不是要模拟 BZ 反应的真实化学。它要回答的是一个更干净的问题:振荡这件事,最少需要哪几样零件?
写成方程只有两行(X、Y 表示浓度):
dX/dt = A − (B+1)·X + X²·Y
dY/dt = B·X − X²·Y
第 2 条(原料浓度被外界按住)是全部关键所在,也最容易被读过去:正是它让这个体系永远到不了平衡态。有人在外面持续供料、持续清废,账单一直有人付——BZ 反应之所以最终会安静下来,就是因为烧杯里的原料终会用完,那时它才真正走进平衡。
把 A 固定,慢慢把 B 调大,会看到一次干净利落的切换:
B 小的时候,X 和 Y 无论从哪儿出发,都会晃几下然后停在一组固定的浓度上(X = A,Y = B/A)。这是个稳态:有物料在不停地流过,但读数不动。
B 越过 1 + A² 的时候,这个稳态失稳了。系统不再停下来,而是走上一条闭合的循环轨道,X 与 Y 一高一低地交替,周期和振幅都稳定。这种自我维持的循环在状态空间里叫极限环(limit cycle)——"极限"的意思是:无论你从哪个初始浓度出发,最后都会被吸到同一条环上去,振幅由方程决定,不由你怎么起步决定。
还有第二种玩法:给 X 和 Y 各加一个扩散项,让它们能在空间里散开。当两者扩散得快慢不同(抑制者跑得比激活者快),系统会长出静止的空间条纹——也就是图灵斑图(Turing pattern)。同一套反应,加不加扩散,给出的是时间上的秩序和空间上的秩序两种东西。
布鲁塞尔器的价值不在预测某个具体反应,而在于它把一件事证死了:持续的振荡不需要任何特殊的、神秘的成分,普通的化学反应加上"体系被持续供料"这一条就够了。它同时给出了一份最小零件清单:一个自催化的正反馈、一个带延迟的负反馈、以及一股不停的物料流。
照着这份清单去看,同一个骨架在很多地方出现:糖酵解过程中的振荡、细胞周期里蛋白质浓度的起伏、捕食者与猎物的数量循环、心肌与神经元的节律放电。它们的化学/生物细节毫不相干,但都装着同样的三个零件。这也是它的核心可证伪预期:只要拿到那个"离平衡多远"的参数并往上推,稳态就应当在某个具体的阈值上失稳,而不是逐渐变得摇晃。