「作用不对等就会动」，这句话要怎样才能被否掉
两个组分混在一起，各自的浓度是一个守恒的标量场。牛顿第三定律要求它们对彼此的作用大小相等方向相反；可是当「作用」是通过某种介质间接传过去的——甲分泌一种东西吸引乙，乙分泌另一种东西排斥甲——这条对等就不再成立。已发表的结论很干脆：把这种不对等写进两个守恒场的动力学，静态的分相图样就会获得一个有限的速度，变成行波。1234
这句话听上去已经很具体了，但它还不能被一次计算否掉。要让它可判，得把三件事分开：
- 互易的那一侧真的不动吗？ 如果只是「没看见它动」，那是观察；如果是「只要能写成某个自由能的梯度流就一定不动」，那是定理，可以直接用来反证。
- 动起来的那一支，怎样与「被外面推着走」区分开？ 一团图样在显微镜下平移，可能是非互易，也可能是有一点点流、一点点温度梯度、一点点漂移。
- 「不对等」这个参数是从哪里来的？ 它从来不是基本作用，总是消去了某个介质之后才出现的有效量。消去这一步做了什么假设、什么时候不成立、以及被消去的那部分在付什么账，都必须交代。
本文把这三件事分别做一遍。用的是已发表的非互易 Cahn-Hilliard 模型（NRCH）及其非线性推广12，两个守恒场写成一个复场 $\phi=\phi_1+\mathrm{i}\phi_2$：
其中非互易强度 $\alpha(\phi)=\alpha_0-\alpha_1|\phi|^2$ 本身随局部浓度变化，在 $|\phi|=\sqrt{\alpha_0/\alpha_1}$ 处回到互易。1 全部数值由本文自写的谱方法求解器给出。
互易的一侧不动：这是定理，不是观察
若两个组分的化学势都能写成同一个自由能的变分，$\mu_a=\delta F/\delta\phi_a$，而动力学是守恒的 $\partial_t\phi_a=\Gamma\nabla^2\mu_a$，那么
$F$ 有下界，于是它必须停下来；在这个类里，持续移动的稳态是被禁止的，不需要看图样长什么样。反过来，一旦看见持续移动，就等于看见了「这套动力学写不成任何自由能的梯度流」。这是一条可以直接反证的陈述，也是本文全篇的支点。
数值上把两边都算出来（一维，随机初值，跑到 $t=300$）：
| 情形 | $F$ 从 → 到 | 单步最大上升 | 上升时间占比 | 后半段注入功率 |
|---|---|---|---|---|
| 互易 $\alpha=0$ | 0.00190 → −0.17151 | $-3.246\times10^{-6}$ | 0.0% | $-1.436\times10^{-21}$ |
| 线性非互易 $\alpha_0=2.5$ | 0.00190 → −0.11628 | $3.809\times10^{-4}$ | 23.7% | 1.3494 |
| 非线性 $\alpha_0=2.5,\ \alpha_1=4$ | 0.00190 → −0.09510 | $6.301\times10^{-3}$ | 47.3% | 0.034411 |
互易那一行的「上升」是负数，也就是从来没有上升过；注入功率在机器精度内为零。另外两行有确定比例的时间在爬坡，而且爬坡是靠持续的功率注入买来的——下面第六节会把这笔账算清楚。
精确解：一条能对表的色散关系
这个模型有精确的行波解。取 $\phi=\rho_q\mathrm{e}^{\mathrm{i}(qx-\omega t)}$ 代入，实部与虚部各给一条：
1 本文的求解器在六个波数上直接积分，振幅与频率与解析式的最大相对差分别是 $5.909\times10^{-11}$% 与 $5.112\times10^{-11}$%。这一节的用处是锚点：后面那些没有解析解的数（混沌态的功率、二维图样的漂移速度）都出自同一套流程。
顺便注意这条色散关系的形状：$\omega\propto q^2$，所以相速 $\omega/q\propto q$。长波几乎不动，短波跑得快。 这一点马上就有用。
判别：非互易的行波与被外面推着走的图样
外加平流（一点点流、一点点梯度）也让图样平移。两者在两个地方分得开。
第一，速度随波数变。 平流给出的速度与波长无关；非互易给出的相速随 $q$ 增大——由 $v=\omega/q=\Gamma q[\alpha_0-\alpha_1(1-q^2/q_0^2)]$，在 $\alpha_1=0$ 的线性极限下正好正比于 $q$，$\alpha_1\neq0$ 时还多一段 $q^3$。本文测到的相速比
而外加平流的同一个比值恒为 1.000。这条判据只需要在同一个样品里比较两个波长，不需要标定任何绝对量。群速度也一样：$V(q)=\mathrm{d}\omega/\mathrm{d}q=2\Gamma q[(\alpha_0-\alpha_1)+2\alpha_1q^2/q_0^2]$，与已发表的平流速度表达式逐字相同1，本文用测出的 $\omega(q)$ 做数值微分核对，内部四个点上最大差 1.92%。
第二，整体反号会掉头。 把 $(\alpha_0,\alpha_1)$ 同时反号，方程给出的频率精确反号：本文测到 $\omega=+0.230713$ 与 $-0.230713$，两者之和 $0$。这是「运动的方向由不对等的符号定」的直接读数。外加平流做不到这一点——把驱动反号，平流方向跟着变的是驱动本身，不是材料的不对等。
被消去的介质：不对等是从哪里来的
$\alpha$ 不是基本作用。它总是这样出现的：甲产生一种化学物质，乙对它有响应；乙产生另一种，甲对它有响应。每一步都是互易的（每个组分与每种化学物质之间都满足作用—反作用），不对等来自「谁产生什么、谁响应什么」这张表不对称，而且只有把化学场消去之后才看得见。已发表的处理里，化学场满足 $-D_c\nabla^2c+\gamma c=\sum_a r_a\phi_a$ 并被绝热消去（假定它比组分快得多）。5
本文照这个结构写一个显式的四场模型（两个守恒组分 $\phi_1,\phi_2$ 加两种化学物质 $c_1,c_2$），但保留化学场的时间导数：
令 $\lambda_1=-\lambda,\ \lambda_2=+\lambda$，绝热消去（$\gamma\to\infty$）给出纯反对称的耦合 $\alpha_0=\lambda r/\gamma$，正是 NRCH。问题是这一步丢了什么。把四场系统在均匀态附近逐波数对角化，与两场 NRCH 比：
| 化学场衰减率 $\gamma$ | 绝热耦合 $\alpha_0$ | 常数 $\alpha_0$ 的最大误差 | 改用 $\alpha_{\rm eff}(q)$ 后 |
|---|---|---|---|
| 1 | 1 | 115.9% | 22.23% |
| 2 | 0.5 | 52.43% | 10.95% |
| 5 | 0.2 | 19.96% | 4.598% |
| 10 | 0.1 | 9.917% | 2.386% |
| 20 | 0.05 | 4.953% | 1.22% |
| 50 | 0.02 | 1.981% | 0.495% |
| 100 | 0.01 | 0.9905% | 0.2487% |
丢掉的东西有名字：消去之后的耦合不是一个常数，而是一条随波数下降的曲线
把常数换成这条曲线，误差立刻降到四分之一左右。物理很直白：化学物质在衰减掉之前只能扩散 $\kappa^{-1}=\sqrt{D/\gamma}$ 那么远，所以比这更短的波长上，甲的信号还没传到乙那里就已经衰减了，有效的不对等就变弱。本文的例子里，最短波处 $\alpha_{\rm eff}$ 只剩绝热值的 0.7117 倍。
这条结论有一个直接的用处：报告一个非互易强度时必须同时报告它在哪个波数上成立。两台装置测出的 $\alpha$ 不能直接比，除非它们的图样波长可比——否则比的是同一条曲线上的两个不同点。而绝热消去成立的条件也就写清楚了：$Dq^2\ll\gamma$。本文的扫描里，常数耦合的误差一直到 $\gamma=10$ 都还有 9.92%，要到 $\gamma=20$ 才降到 5% 以下。
资源账：持续移动要一直花钱
守恒动力学里，电流 $J_a=-\Gamma\nabla\mu_a$ 耗散的功率密度是 $\Gamma\sum_a|\nabla\mu_a|^2$。在稳定的行波上 $F$ 不随时间变，于是注入功率恰好等于耗散功率。对精确平面波这个数有闭式：行波上平衡部分的化学势恒等于零（$-1+\rho_q^2+Kq^2=0$），剩下的全是非互易那一项 $\mathrm{i}\alpha\phi$，于是
本文在线性情形（$\alpha_1=0$，此时闭式最干净：$P=\Gamma\alpha_0^2q^2\rho_q^2$）上逐波数核对，数值与它最大差 $2.580\times10^{-11}$%。这条曲线在 $q=q_0/\sqrt{2}=0.7071$ 处最大，值 1.5625；而互易极限下它恒为零——下面那张表里互易那一行测到的注入功率就是 0.0000。
$\alpha_1\neq0$ 时峰会移位，而且多出一件事：当 $\alpha_1>\alpha_0$ 时，$\alpha_0-\alpha_1\rho_q^2$ 会在 $\rho_q^2=\alpha_0/\alpha_1$（即 $q^2=q_0^2(1-\alpha_0/\alpha_1)$）处穿零，功率曲线因此在中间也回到零——那正是局部恢复互易的那个波长。下面那个「混沌反而省能」的结果，用的是同一个 $\alpha(\phi)=0$ 的条件。
三个稳态的对照更有意思：
| 稳态 | 注入功率 |
|---|---|
| 行波 $(\alpha_0,\alpha_1)=(4,1)$ | 1.1506 |
| 起泡 $(2.5,4)$ | 0.0305 |
| 互易 $(0,0)$ | 0.0000 |
看上去更「活跃」的时空混沌态，功率反而只有规则行波的约 38 分之一。原因在闭式里已经写着：功率正比于 $\alpha^2$，而起泡态之所以混沌，恰恰是因为 $\alpha(\phi)=\alpha_0-\alpha_1|\phi|^2$ 在大片区域里被拉到接近零——局部恢复互易的那些「互易颗粒」不耗能。「看起来更乱」与「耗散更多」是两件事，这一条在把这类模型往真实体系上套时值得先算一遍。
什么时候它不再是一支行波
再往上加 $\alpha_1$，行波让位给时空混沌。已发表的判据是一维最大李雅普诺夫指数变号。1 本文重做这个扫描，先量分辨率：同一组参数换三个随机种子，指数的散布是 0.00265，小于这个量的读数只能算「零」。结果是
| $\alpha_1$ | 0 | 1 | 2 | 2.5 | 3 | 3.5 | 4 | 5 | 6 |
|---|---|---|---|---|---|---|---|---|---|
| $\lambda_{\max}$ | -0.0042 | -0.0009 | +0.0269 | +0.0423 | +0.0459 | +0.0539 | +0.0301 | +0.0063 | +0.0003 |
峰值 0.0539 落在 $\alpha_1=3.5$，两侧都回到零——形状与已发表的图一致。但有一处不一致，而且是定性的：本文读到的是「零转正」，不是「负转正」。
这不是精度问题，是对称性问题。方程对全局相位旋转 $\phi\to\mathrm{e}^{\mathrm{i}\gamma}\phi$ 严格等变——这条对称性是已发表工作自己写出来的1——所以只要稳态不是 $\phi\equiv0$，相位方向就是一条精确的中性方向，对应一个恰好为零的李雅普诺夫指数；因此 $\lambda_{\max}\ge0$ 恒成立，不可能为负。本文沿这个方向直接量了一次，增长率 $-1.608\times10^{-4}$，小于上面那条分辨率。
那已发表的负值是从哪来的？它写出的定义与本文算的其实是同一个量：
1 分歧只可能出在「$T\to\infty$ 这个极限在数值上怎么做出来」。本文用的是逐段把孪生轨道拉回初始距离的标准做法，每一段量的都是瞬时的分离率；而上式照字面只做一次——放一个小扰动，跑 $T$，量一次比值。在行波这一侧两者并不等价：收缩方向会把扰动压到中性流形上，$|\delta\phi(t)|$ 于是停在一个小于 $|\delta\phi(0)|$ 的常数上，读数就成了 $(1/T)\log(\text{常数})$——一个负数，而且随 $T$ 变大而消失。本文照字面把这条曲线量了出来：
| 测量时长 $T$ | 25 | 50 | 100 | 200 | 400 | 800 |
|---|---|---|---|---|---|---|
| $\alpha_1=1$（行波一侧） | -0.0382 | -0.0180 | -0.0084 | -0.0023 | -0.0008 | -0.0003 |
| $\alpha_1=3.5$（混沌一侧） | +0.0333 | +0.0519 | +0.0369 | +0.0481 | +0.0458 | +0.0231 |
行波一侧从头到尾都是负的，量级随 $T$ 按 $T^{-1.45\pm0.06}$ 衰减（纯 $1/T$ 给 $-1$，这里更陡，是因为被压掉的那些分量还在缓慢恢复）；混沌一侧则始终为正（$+0.0231$ 到 $+0.0519$）。负号因此不是一个负的指数，而是「$T$ 还不够大」的读数，它与 $\lambda_{\max}\ge0$ 并不矛盾。顺带一提，照字面做在混沌一侧同样会坏，只是坏法不同：到 $T=800$ 时两条轨道的距离已经涨到场本身的 1.69 倍，早出了线性区，于是读数也开始往下走——从 $T=200$ 的 $+0.0481$ 掉到 $T=800$ 的 $+0.0231$。这正是要逐段重整化的理由。
已发表工作没有写出它实际用的 $T$、初始扰动的大小、以及是否逐段重整化，所以本文只能说：这是那一支负值最可能的来源，而不是已证实的结论。但「行波一侧的真实指数应当是零而不是负」这一点，在本文的模型里是从对称性直接推出来的。
二维取三个代表性参数点——起泡 $(2.3,4.6)$ 与互易颗粒／液滴 $(5,4)$ 是已发表工作在状态图上标出的点1，行波那一组 $(4,1)$ 是本文自己在行波区里取的。本文不只给快照，还量了图样到底在不在动：跟踪主导傅里叶模的相位，得到
| 参数点 | 主导波数 | 该模占谱重 | 实测漂移 | 单波闭式 |
|---|---|---|---|---|
| 行波一侧 $(4,1)$ | 0.3125 | 5.0% | +1.0272 | +0.9680 |
| 起泡 $(2.3,4.6)$ | 0.2652 | 1.0% | +0.2345 | -0.5241 |
| 互易颗粒与液滴 $(5,4)$ | 0.2577 | 1.8% | +0.9453 | +0.3261 |
三块都在动，但只有第一块能被「一支行波」描述：它的主导模占谱重 5.0%，实测漂移 +1.0272 与单波闭式 +0.9680 差 6%。另外两块的主导模只占百分之一二，闭式给出的数就不再有意义——起泡那一行连符号都反了（实测 +0.2345，闭式 -0.5241）。「图样在动」这件事在这三块里都成立，但「它以多快的速度在动」只在第一块里是一个有定义的数。这正是把「持续移动的图样」当成单一现象来讨论时最容易滑过去的地方。
边界、反例与可检查的预期
没有覆盖的： 全部数值都在零温、确定性、无涨落的方程上（已发表工作里涨落会带来另一类效应，本文完全没有做）；没有流体力学耦合，没有非零平均组成，一维与二维各只做了一个盒子尺寸（$256$ 与 $128^2$，边长 $32\pi$），也没有把任何参数对真实体系标定。因此本文所有绝对数值只在模型内部有意义，能往外搬的是无量纲的结构性结论：梯度流禁止运动、相速随波数变、$\alpha_{\rm eff}(q)$ 的形状、以及功率闭式的形式。
四场模型的局限： 它是本文自己写的，不是任何一篇已发表工作的模型。它照 5 的结构（化学场由组分产生、满足带衰减的扩散方程、被绝热消去）搭起来，但把两种化学物质分开、并保留时间导数，是为了让「消去」这一步的代价可以被单独量出来。它没有对应任何具体实验体系，也没有包含化学反应的自由能收支——所以本文能说的是「消去这一步在什么条件下改变谱」，不能说「真实体系的非互易强度就是这条曲线」。
一处与已发表结论不同的读数： 一维最大李雅普诺夫指数，已发表工作报告在行波一侧为负1，本文读到的是零。两处用的是同一个定义，本文给出的理由是对称性（全局相位旋转是精确的中性方向，因此 $\lambda_{\max}\ge0$），并把「照字面不重整化、有限 $T$」这条路径直接量了出来——它在行波一侧确实给负数，且随 $T$ 按幂律消失。需要说清楚的是，本文并不知道已发表工作实际用了多长的 $T$、多大的初始扰动、有没有逐段重整化，这三项它都没有写出来；所以本文给的是一个可以直接去核的候选解释，不是对它实现方式的断言。
一处方法上的弱点： 二维那三块都是从随机初值跑到 $t=400$ 的单次实现，没有做系综平均，也没有检查是否已到稳态。主导模占谱重只有百分之几，说明第二、三块远没有被单个波描述——本文因此只把它们的漂移速度当作「有没有在动」的证据，不当作定量结果。
一处自查出来的假信号： 第七节的扫描第一版把暖机时间取成 $t=120$，行波一侧读到一个明显超出分辨率 0.00265 的正值。单独做收敛测试才发现要到 $t\approx600$ 才落回零：那个正值量的是初始暂态还没散掉时两条轨道的分离，不是稳态的指数。现行脚本取 500/500，行波一侧才落回零。（那次读数出自一个只做收敛测试的原型脚本，现行脚本不再产生它，所以这里不给具体数值。）这与第五节讲的是同一件事——一个没被写出来的约定（暖机多久、在哪个波数上），会让读数不是它看上去的那个量。程序本身也有过一处：积分因子四阶 Runge-Kutta 的最后一段少写了一个因子，症状很不显眼——振幅全对，只有频率错得离谱。这正是第三节那个解析锚点存在的理由：只有一半的量对得上，不足以说明程序对。
可检查的预期：
- 同一个样品里比两个波长。 这是本文最容易做的一条：若图样的移动来自非互易，相速应当随波数增大（线性极限下正比于波数；本文用的参数还带一段 $q^3$，6 倍的波数比给出 8.165 倍的速度比）；若来自外加平流，两个波长应当给同一个速度。这条判据不需要标定任何绝对量，只要能在同一次实验里同时看到两个波长的结构。
- 把介质的弛豫速率调慢。 化学介质的衰减率 $\gamma$ 可以通过改变酶浓度、缓冲液或温度来调。本文预言：$\gamma$ 变小时，短波处的有效非互易强度会先掉下来，也就是图样的特征波长会向长波移动，而且移动的方式由 $\alpha_{\rm eff}(q)=\lambda r/(\gamma+Dq^2)$ 定；若观察到波长不随 $\gamma$ 变，说明介质并没有被这样消去，本文第五节的图像在那里不成立。
- 功率与波长的非单调关系。 维持行波的功率随波长非单调：线性情形（$\alpha_1=0$）在 $q=q_0/\sqrt{2}$ 处最大，两端都趋于零。这意味着最费能的不是最细的图样——把波长压到接近界面宽度反而省能。若某个体系能同时测出耗散（比如量热）与波长，这条非单调性是可以直接对表的。
- 混沌态更省能。 本文测到起泡态的功率只有行波态的约三十八分之一，机制是局部把非互易拉回零。这条给出一个反直觉的预期：在 $\alpha_1$ 的扫描上，功率应当在进入混沌区时下降而不是上升。这一条本文只在一维两个参数点上量过，值得系统扫一遍。
如果被推翻： 若有人在一个确实能写成守恒梯度流的体系里观察到持续移动的稳态，本文第二节的定理就错了（更可能是那个体系其实写不成梯度流，而这正是判据的用法）。反过来，本文没有、也不能证明任何真实体系的运动来自非互易——本文给的是判据与代价，不是身份认定。
数据与代码
复算脚本 research/calculations/nonreciprocal_patterns_and_motion.py，输出在 research/assets/2026-09-17-nonreciprocal-coupling-and-travelling-patterns/（results.json 与 11 个 CSV、4 张图）。随机种子全部固定，两次完整运行在本文引用的所有有效数字上一致。原始文献在 evidence/2026-09-17-nonreciprocal-coupling-and-travelling-patterns/。