一个能逐项算清的最小例子
描述一杯水里的一颗小球，本来要跟踪水里每一个分子。实际上人们只写一个方程：黏性摩擦乘速度，加一个随机力。这个替换在很多场合好用到让人忘了它是个替换。
要看清它什么时候不成立，需要一个既真实又能算到底的系统。Ginot 与 Bechinger 2025 年的实验提供了一个：直径 2.73 μm 的二氧化硅小球，泡在 8 mM 的十六烷基氯化吡啶／水杨酸钠溶液里，这种溶液会形成蠕虫状胶束的缠结网络，是有明确松弛时间的黏弹性流体；小球被一束聚焦激光抓住，激光焦点的位置 λ(t) 由压电台按四步来回移动。1
四步是：①以 0.186 μm/s 匀速把阱向右拖 3 s；②停 $t_{neq}$；③以同样速度拖回；④等 50 s 让系统完全回到平衡。整套协议重复约 100 次取平均，相机 100 帧/秒，定位精度约 5 nm。
只有两个量能直接测：小球位置 x 和阱位置 λ。阱势是 $V(x,\lambda)=\tfrac12\kappa(x-\lambda)^2$，刚度 κ 由静止时小球位置的概率分布独立标定。由此可以算外界对系统做的功和势能，两者之差就是与环境交换的能量：
如果环境是普通黏性流体，$\delta E$ 就是散掉的热，只能单调增加。实验在 50% 水-甘油溶液里验证了这一点。1
在黏弹性流体里不是。步骤③一开始，$\langle E\rangle$ 先降了下去，降幅约 6.7 $k_BT$，1.7 s 后才转回上升——本文从作者存档的原始曲线上量得这个数，协议开始前那段平衡期的曲线起伏只有 0.052 $k_BT$。2
这个下降就是本篇要解释的现象。它不是热涨落——单条轨迹的 E 本来就上下跳动，甚至可以为负；这里下降的是上百条轨迹的平均。它也不是熵在减少：阱刚度恒定，小球位置的分布宽度在整个协议里不变（下文会算），所以系统熵没有变化。
消去环境，到底消去了什么
作者用来解释的模型只加了一个变量：一个虚拟的"浴粒子"，位置 $x_b$，通过一根弹簧 $\kappa_b$ 挂在小球上，自己有摩擦 $\gamma_b$。1
这是一个完全线性的两变量系统，可以精确消元。把第二个方程解出来代回第一个（对含 x 的卷积做一次分部积分），得到
这就是广义 Langevin 方程。它是精确的——不是近似，因为原系统是线性的。本文把这三项逐点积回去，与用矩阵指数精确推进的两变量模型比较：在整段协议里最大偏差 0.019 nm，而小球的位移幅度是 404 nm；这点残差还是 GLE 那边的时间离散造成的。
值得停下来看的是这个式子右边多出来的三样东西，而不是通常被强调的一样：
- 记忆核 $K$：摩擦不再是一个数乘当前速度，而是对过去速度的加权求和。
- 有色噪声 $\eta$：它的自相关不再是 δ 函数，而与 K 由涨落耗散关系绑定。
- 初值项 $K(t-t_0)\,[x_b(t_0)-x(t_0)]$：一个随时间指数衰减的确定性力，大小由起步那一刻环境被拉开多远决定。
第三项在通常的写法里看不见，因为默认起点是完全平衡态，那时它平均为零。而实际用约化描述时，人们几乎总是从某个非平衡时刻起步——"系统现在处在状态 x，往后会怎样"。这一步就把它丢了。
第二项也不是可以单独换掉的。本篇后面算的都是系综平均，噪声不进入均值轨迹；它的作用在涨落上。把记忆核留着、噪声却按白噪声处理，平衡分布就会错：用零频摩擦定噪声幅度，位置方差变成应有值的 12.8 倍；用裸摩擦定幅度，又只剩 0.65 倍；只有与核配套的有色噪声才给回 $k_BT/\kappa$。记忆和噪声必须成对替换。
代进这个实验（暂停 1 s 的那一组）：步骤③开始时，弹簧被拉开了 279 nm，储在环境里的弹性能是 10.05 $k_BT$，初值项给出的力是 296 fN——同一时刻阱本身施加的力是 302 fN。被略去的那一项，和唯一被看见的那个力一样大。 暂停越久这一项越小，衰减的时间常数就是 $\tau_b$。
在频率域上这件事更直白。对上式做傅里叶变换，摩擦系数变成
ω→0 时它等于 $\gamma+\gamma_b$，也就是零剪切摩擦；ω→∞ 时它退回裸摩擦 γ，因为环境完全跟不上。无记忆闭合做的事，是把这条曲线整体换成它的零频值。作者拟合的参数是 γ = 0.34、$\gamma_b$ = 6.33 μN·s/m、$\kappa_b$ = 1.06 μN/m，于是 $\tau_b$ = 5.97 s，而 $\gamma+\gamma_b$ = 6.67 μN·s/m 对应黏度 259 mPa·s，与论文另处引用的零剪切黏度约 300 mPa·s 相符，这是一个独立的自洽检查。
本协议的主频约 π/3 ≈ 1.05 rad/s。在这个频率上，$|\hat\gamma|$ 只有零频值的 0.166，相位角 63°。也就是说，无记忆闭合在驱动自己的频率上，把摩擦的模当成了实际值的六倍，而且完全丢掉了那个不同相的分量。要让这个替换的误差小于 5%，频率得低到 0.0088 rad/s——周期约 710 s。
一条不靠调参数就能用的判据
上面是模型内部的比较——好不好，取决于那个模型对不对。更有说服力的是：能不能不挑某一个具体模型，直接判断无记忆闭合这一整类描述在这个实验里成立不成立？
可以，而且判据就是那条下降的曲线。设可见变量服从任意无记忆闭合
阱刚度 κ 恒定，内能只有阱势 V。把位置拆成平均值与涨落 $x=\bar x+y$，涨落满足 $\gamma_{\rm eff}\dot y=-\kappa y+\xi$，与 λ 无关，所以只要一开始 $\langle y^2\rangle$ 处在平衡值 $k_BT/\kappa$，它就一直待在那里。于是
结论不依赖 $\gamma_{\rm eff}$ 取什么值，也不要求它是常数：只要它是正的，含时的一样成立（含时只是把涨落方程做了时间重标度，平衡方差 $k_BT/\kappa$ 不含 $\gamma_{\rm eff}$，所以照样不动）。若 $\gamma_{\rm eff}$ 还依赖状态，上面对涨落的拆分不再严格，但均值方程那一步 $\langle\dot E\rangle=\gamma_{\rm eff}\dot{\bar x}^2$ 仍然成立。这一整类描述都给出单调不降的 E。
数值上核过：把零剪切摩擦代进去跑同一协议，逐点算出的 $\mathrm d\langle E\rangle/\mathrm dt$ 最小值是 +0.000187 $k_BT$/s，与 $\gamma_{\rm eff}\dot{\bar x}^2$ 的最大偏离是 0.000329，而这个量的峰值是 19.1；让摩擦在协议中按正弦在 1.33 到 12.0 μN·s/m 之间来回变 9 倍，E 仍然逐步不降。
所以这不是"无记忆模型拟合得差一点"，而是它连曲线的形状都做不出来。把 $\gamma_{\rm eff}$ 当自由参数去最小二乘拟合整段实验 E(t)，最优值是 0.636 μN·s/m，残差均方根 5.16 $k_BT$；保留一个环境变量的两变量模型残差是 0.76 $k_BT$，小七倍。而无论怎么调，无记忆那一侧的下降量恒为 0。
顺带一提，这个最优的 0.636 μN·s/m 离零剪切值 6.67 很远，倒是接近上一节在协议主频上算出的 $|\hat\gamma|$ ≈ 1.11 μN·s/m。这正说明拟合出来的"有效摩擦"其实是在替协议的主要频率成分做局部近似，而不是任何一个材料常数。
实验测到的下降是 6.7 $k_BT$。作者自己给出的同一暂停时长下 ER 的标准误是 0.36 $k_BT$，所以这个下降约是其自身误差棒的 19 倍。它同时也可以换个说法：步骤①注入环境的能量是 20.9 $k_BT$，回收了 32%，与论文正文说的"约 30%"一致；在模型里这 5.69 $k_BT$ 相当于此刻环境中储存的 10.05 $k_BT$ 弹性能的 57%。
两变量模型对实验曲线的整体符合程度：功的残差均方根 0.65 $k_BT$（曲线幅度 28.5），势能 0.18（幅度 7.5），交换能量 0.76（幅度 23.3）。从存档的阱位置数据独立拟出的驱动速度是 0.1873 μm/s，作者代码里写的是 0.1862，差 0.5%。
反过来不成立：E 单调并不证明没有记忆。 暂停足够长时环境已经松弛完，同一个系统也给出单调的 E（图 3 左里 $t_{neq}/\tau_b$ 大的那一端，可回收能量趋于零）。这个判据只在"环境来不及松弛"的窗口里有判别力，而这个窗口正是要靠改变暂停时长去扫出来的。
需要说明判据成立的条件，因为它们不是自动满足的：阱刚度必须恒定（刚度随时间变，系统熵会变，$\langle\dot E\rangle$ 就可以为负）；位置方差要处在平衡值（本文用 20000 条带噪声的两变量轨迹核过，整个协议里 $\kappa\langle y^2\rangle/k_BT$ 在 0.998 到 1.013 之间，平均轨迹与确定性解的差不超过 0.71 nm——线性系统里对浴粒子作边缘化正好给出方差 $k_BT/\kappa$，与环境无关，这也是熵不变这条论证成立的原因）；以及 E 必须按 $W-\Delta V$ 定义，用的是真实的阱势而不是某个重新拟合的有效势——允许有效势自由变形，等于偷偷加回了一个储能变量。
缺的到底是几个数
既然无记忆不行，另一头的问题是：补回多少才够。
在这个系统里答案很干净：一个数。步骤③起点上，如果约化描述知道环境被拉开了多少，单变量的 GLE 就与两变量模型逐点重合（最大偏差 0.017 nm）；不知道，只带位置这一个状态重新起步，步骤③末的位置差 71.6 nm，可回收能量从 5.69 掉到 0.17 $k_BT$。
把"以为的拉伸量"从 0 连续调到真值：末位置误差严格线性地从 71.6 nm 降到 0（对直线的最大偏离在 0.001 nm 以下，线性系统本该如此），可回收能量则单调上升且下凸，从 0.17 升到 5.69 $k_BT$——它是力乘位移，两个因子都随拉伸量长。缺的信息就是这一个标量，而不是环境里那些数不清的分子。
这给出两条等价的补法，也正是这个问题的两种常见解法：
- 保留历史：从一个确实平衡的起点出发，把整段速度史都带着，初值项自动为零，记忆核描述就是精确的。代价是这个描述不是局域的——你不能中途接管。
- 加状态变量：把 $x_b$ 显式带上，方程重新变成局域的、无记忆的（这就是马尔可夫嵌入）。代价是多了一个测不到的量。
第二条在实验上不可行——$x_b$ 是虚拟的，直接测不到；但它可以从可测的历史反推。两条路给出同一个答案，因为它们编码的是同一份信息。
补几个，取决于你要回答什么
上面说"一个数"，是在一个环境变量的模型里。真实流体有几个模式，是另一个问题，而且答案不唯一。
作者自己给了两个版本：一个浴粒子和两个浴粒子。用他们的参数跑同一批暂停时长，与实验对照：
| 量 | 一个环境变量的残差均方根 | 两个环境变量的残差均方根 |
|---|---|---|
| 步骤③的功变化 ΔW | 4.42 $k_BT$ | 0.49 $k_BT$ |
| 步骤③的势能变化 ΔV | 1.69 | 0.39 |
| 步骤③的交换能量变化 ΔE | 2.75 | 0.31 |
| 可回收能量 ER | 0.66 | 0.40 |
加第二个变量把能量收支的误差压下去了一个量级，这和论文的说法一致。但同一次改进对不同的量效果差得远：ER 这个更细的特征只从 0.66 改善到 0.40 $k_BT$，在最长的两个暂停时长上，一变量模型偏离实验 6.0 个标准误，二变量仍偏离 2.9 个。也就是说，"够不够用"必须连着"用来回答什么问题"一起问；同一个约化描述可以在能量总账上准确，在某个更细的特征上仍然系统性偏小。
这里也遇到一处只能如实记下的不一致：同一个暂停时长（1 s），从作者存档的图 1c 曲线量出的 ER 是 6.73 $k_BT$，而图 2a 表里给的是 5.68 ± 0.36；那一组的 ΔE 也比图 1c 曲线小约 19%。两批存档数据的处理不完全相同，本文没有原始轨迹可以判断哪一套更该采信，因此模型对照一律用作者表里的值，而"下降存在"这一定性结论不依赖于这个差别。
需要多快的环境才能真的用无记忆闭合？固定零频摩擦、只改变 $\tau_b$，量无记忆闭合相对完整模型的最大位置误差：要小于 5%，需要 $\tau_b\lesssim0.20$ s，即 $\tau_b/t_d\lesssim0.068$；要小于 1%，需要 $\tau_b/t_d\lesssim0.013$。本实验所处的位置是 $\tau_b/t_d\approx2$，误差 46%。"环境比驱动快一些"远远不够，得快一两个量级。
反过来说，这个条件在日常情形里往往轻松满足，所以标准 Langevin 方程才那么好用。同一颗 2.73 μm 的球换到水里：涡量扩散过球半径要 1.9 μs，颗粒自身的动量松弛要 0.83 μs，比 3 s 的驱动短约六个量级——按上面扫描在小比值处的斜率外推，误差在 $10^{-7}$ 量级。需要说明的是，水里的记忆是流体动力学型的（长时幂律尾），与本篇的单指数核不同类，这里只做时标量级的比较，不是把上面的阈值直接套过去。真正的分界不在"有没有记忆"，而在环境的最慢过程与你关心的时间尺度谁更长。
拟合得好，不等于参数有物理含义
最后一件事关系到这类约化描述能不能跨条件复用，答案是要相当小心。
第一，观测到的松弛时间是耦合系统的本征值，不是流体单独的属性。 同一组人的回弹实验里，两个浴粒子的模型在"探针被光镊夹住"和"探针自由"两种条件下给出两套不同的时标：把他们报的参数代回本征值公式，探针自由时是 0.335 s 与 3.23 s，探针被夹住时是 0.733 s 与 4.4 s，与论文发表的 0.34/3.2 和 0.7/4.4 相符。4 原因是浴粒子之间只经由探针耦合，探针被夹住就等于把它们解耦。后续实验进一步表明，靠设计驱动协议可以选择性地压掉其中一个模式，甚至让回弹几乎完全消失——同一种流体，看上去有两个时标、一个时标还是没有时标，取决于你怎么推它。3
第二，被协议弱约束的那部分参数，不能当材料常数。 同一标称流体（8 mM、同一粒径、25 °C）在两篇 2025 年的论文里被拟合成三套参数：
| 拟合 | 裸摩擦 γ（μN·s/m） | 零频摩擦（μN·s/m） | 蕴含黏度（mPa·s） | 夹住探针时的松弛时标（s） |
|---|---|---|---|---|
| 一个浴粒子 | 0.34 | 6.67 | 259 | 5.97 |
| 两个浴粒子 | 0.38 | 18.63 | 724 | 4.85；60.94 |
| 回弹实验的两浴模型 | 0.44 | 12.28 | 477 | 2.53；17.77 |
直接受小球快速运动约束的裸摩擦，三套只差 1.3 倍；而最慢的那个时标差 10 倍，由它主导的零频摩擦对应的黏度差 2.8 倍。两浴粒子拟合里那个 60.9 s 的模式独占了零频摩擦的 79%，而整个协议只有 7 s 长——无记忆闭合唯一的那个参数，恰恰由实验最没有约束住的那个模式决定。 胶束溶液不同批次本身也有差别，所以这里不说哪一套错，只说这些数不适合当材料常数搬到别处。
第三，模式数目本身就不可辨识。 用两浴粒子模型生成的 ER 随暂停时长的曲线（0.45 到 6.43 $k_BT$）作为目标，去拟一个单浴粒子模型：最佳单一模式 $\tau_b$ = 9.92 s，残差均方根 0.158 $k_BT$，只有目标幅度的 2.5%。这个"最佳单一时标"既不是 4.85 s 也不是 60.9 s，而它蕴含的黏度是 306 mPa·s，与被模仿对象的 724 mPa·s 差 2.4 倍。在一个观测量上拟合得很好，参数却指向完全不同的物理图像。
边界、反例与可检查的预期
这个例子成立的条件。 整套结论建立在线性上：弹簧线性、摩擦线性、噪声高斯。作者自己核过在本协议的驱动速度下流体处于线性区间；但更宽速度范围的流动曲线显示这种溶液是持续剪切变稀的3，也就是说 γ 本身随速度变。s4 也转述了更早的工作：在 30 μm/s 量级观测到回弹幅度饱和，那是溶剂能储存的弹性能有限造成的，线性弹簧描述不了。4 所以上面每一个数都属于"某个驱动窗口内的有效描述"。
"加变量总能补回来"不是普遍成立的。 一个胶体耦合到高斯场的模型给出的有效方程带非线性记忆项，而且记忆核依赖于粒子受到的光阱约束——同一个粒子换个阱刚度，"环境的记忆核"就变了。5 这说明"把环境写成一个与系统状态无关的线性记忆核"本身就是一个近似，只是在本篇的线性系统里恰好精确。
初值项在一般情况下更麻烦。 用时间依赖的投影算符处理非稳态过程，得到的方程虽然形似 GLE，但记忆核本身含时，涨落力隐含地依赖初始条件，涨落耗散关系也要改写。6 本篇之所以能把它压缩成"一个数"，是因为线性加单指数核这个特殊结构。
约化描述的状态不完整，会以协议依赖的形式暴露出来。 一项 2026 年的理论工作研究了黏弹性介质里只重置探针、不重置介质的随机重置：稳态分布出现非指数尾巴，而且涨落依赖于返回时的速度——在普通布朗运动里它不依赖。7 这是同一件事的另一个面：只要环境保留自己的状态，"把可见变量放回某处"就不等于"把系统放回某处"。
没有验证的部分。 本文没有复算作者的梯度下降拟合流程，用的是他们发表的参数；没有触及原始轨迹（作者说明只应要求提供），用的是存档的平均曲线；三篇支撑文献里，回弹实验读的是 arXiv 版而不是正式刊出版；非稳态 GLE 那篇读的是 arXiv v5，正式版换了标题，本文没有逐字比对差异。上文那处 6.73 与 5.68 的不一致也没有解决。
可检查的预期。
- 在同一装置上做刚度调制协议（κ 随时间变而阱心不动），本文的单调性判据不再适用，$\langle E\rangle$ 在纯黏性流体里也可以下降。因此用 E 的非单调性判断记忆，必须先确认刚度恒定；这一点可以直接测。
- 若把暂停时长固定、只改变驱动时长 $t_d$，作者定义的工作比 $\Sigma$（步骤③的做功变化除以步骤①的做功变化，全部回收对应 $-1$、全部散掉对应 $+1$）在黏弹性一侧有极小（本文复算给模型 $t_d\approx2$ s、实验 3 s，七个点的残差均方根 0.067），而同零频黏度的无记忆对照随 $t_d$ 单调上升、没有极小；两条线在 $t_d\approx1$ s 交叉。这个"极小是否存在"比"曲线拟合得好不好"更能区分两类描述。
- 若两浴粒子拟合里那个 60.9 s 的模式是真实的，用一个长达数百秒的低速往返协议应该能直接看到它，并把零频摩擦独立定出来；若那只是拟合的冗余自由度，长协议测得的零频摩擦会接近 259 而不是 724 mPa·s。这是可以做的实验。
- 本文的单调性判据可以搬到任何"读数 + 恒定二次势"的装置上：在磁镊、AFM 或原子阱里，只要势的形状在协议中不变、初始涨落已平衡，平均交换能量的下降就是环境保有状态的直接证据，不需要先建模型。