十秒钟的振动，三天后的裂纹
把 CaCO₃ 粉末和水按体积分数 44% 搅成浆，倒进直径 500 毫米的浅盘，放在恒温恒湿房间里让它干燥三天。通常得到的是各向同性的胞状裂纹——碎片大小相近、取向随机。但 Nakahara 和 Matsuo 2005 年报告了一个反直觉的做法：倒完盘后立刻把盘子水平振动十来秒，然后再也不碰它，最终裂纹会沿垂直于振动方向排列成层状图样；若改作角向振荡，裂纹则呈放射状。1
这组对照把"记忆"两个字坐实了。其一，层状方向与容器形状无关——圆盘也能产出单向层状纹，所以不是几何边界在引导。其二，记忆极其耐久：把振动过的浆体密闭存放（湿度饱和、干燥暂停），一个月后打开再干燥，裂纹仍垂直于一个月前的振动方向。1 其三，记忆存在浆体内部：振动后表面平滑无皱纹，且倒入 12 小时后切掉上三分之一表层，裂纹图样不受影响——表面痕迹假说被直接排除。同一种效应在高岭土浆体中复现，说明它不是碳酸钙的特例。1
什么时候写不进去？形态相图给出两个方向的边界。水太多（体积分数低于 25%）时浆体是牛顿流体，永远只得各向同性纹；振动强度也有一个窗口：太弱（低于由屈服应力画出的边界线）时浆体根本不动，无记忆；太强时表面波和湍流出现，宏观流动把刚写入的微观结构吹散，也无记忆。只有中间一段——振动足以让浆体发生塑性剪切、又不至于搅成宏观流动——记忆写入成功。1 与之对应，流变测量显示同一体系在 25% 体积分数以上出现非零屈服应力并随浓度急剧上升：记忆效应的存在区间恰好是"浆体是粘塑性流体"的区间。1
本文追问三件事：被写入的"东西"到底是什么（残余应力、颗粒网络各向异性，还是别的）；它为什么能挺过一个月、又为什么能被轻易改写；裂纹又怎样把如此微弱的内部差异读成一目了然的图样。框架性证据来自五篇全文精读的文献：Nakahara–Matsuo 的原始实验 1、Ooshida 的连续介质理论 2、Nakahara 等含记忆重写实验的 2018 年综述 3、Akiba–Shima 的独立复现 4 与 Nakahara–Matsuo 2006 年的方向转变快报 5；其余文献按摘要级引用、只取方向性结论（书目已逐条核验）。文中裂纹图样来自本库自建的弹簧网络示意计算，只核定性结构。
机制竞赛：写入的是什么
候选载体有两个，外加一个已被排除的干扰项。先说结论：两种载体目前都有直接观测，但它们由同一次振动同时诱发，因果贡献至今没有被实验分开——这是领域的真实现状，本节按证据分别称量。
塑性残余应力（连续介质语言）。这一支的理论链条最完整。Otsuki 在 2005 年与 Nakahara 实验同年独立提出黏塑性模型：振动写入塑性变形记忆，并可解析给出首条裂纹的位置。6 Ooshida 2008 年的连续介质版本给出一条干净的解析链条：浆体是各向同性弹塑性体，应力服从 Maxwell 型弛豫，弛豫时间在应变能低于屈服阈值时趋于无穷——低于屈服，形变被"冻结"。对坡面稳态流，模型证明流动必然产生流向张力参数 $\alpha>0$（由 $\varepsilon=2(e^{\alpha}-1)$ 与 $\varepsilon>0$ 立得），对应法向应力差 $\sigma_{xx}-\sigma_{zz}=2S\sinh\alpha$；停流后低于阈值的区域弛豫时间无穷大，一部分 $\alpha$ 被永久冻结。干燥时浆体各向收缩而盘底约束不动，流向已被预张紧：流向拉应力最大，裂纹面垂直于流向张开、裂纹线垂直于流向延伸——与实验观察一致（原文第 10–11 页用"键"表述这一环节时措辞有歧义，本文只保留应力与裂纹面的关系）。理论还顺手解释了相图内部结构：残余 $\alpha$ 按 $(\sigma_Y/S)^2$ 标度，浆越稠记忆越强，与实验一致。2 1 值得注意的是这个机制根本不需要振荡——坡面流的初步实验同样观察到记忆，写入的本质是塑性流动而非振动本身；Ooshida 次年把理论推广到了三维。2 7
这一支内部还有一次干净的自我裁决。准线性版本（Otsuki 式）要求变形在流向上不均匀、残留下 $\partial X/\partial x>0$ 的拉伸区；但实验里所有区域的裂纹都垂直于振动方向，几何上不可能处处都是拉伸区。非线性版本（Ooshida 式）里驱动 $\alpha$ 的项正比于 $\sigma_{xz}^2$、与剪切符号无关，振荡中单调累积，可以处处为正。Nakahara 等 2018 年的记忆重写实验把两者分开了：先沿一个方向振 60 次写入，再沿偏转 45° 的方向追加振动——追加不足两次时，只有"最后一次运动中被拉伸"的区域记忆被改写（准线性机制主导）；追加满两次，全场的记忆就被新方向改写完毕（非线性的 $\sigma^2$ 累积主导）。3
实验端的直接测量也跟上了：Kitsunezaki 等 2016 年用弹性底板的弯曲反演浆体应力，首次直接测到振动诱发的应力各向异性。8 但同一组人在 2018 年综述中自承一个时程上的麻烦：振动刚停时应力差 $\sigma_{xx}-\sigma_{yy}$ 并不大，是干燥过程中逐渐增大的。3 残余应力是记忆的"种子"还是"全部"，这条测量单独回答不了。
颗粒网络的统计各向异性（颗粒语言）。Nakahara 和 Matsuo 2005 年的微观猜想——振动诱发非弹性碰撞、产生沿振动方向的密度涨落、涨落改变局部键强——当时是假说级的，作者明说电镜看不到。1 但后续成像把颗粒各向异性从猜想变成了观测：2017 年的颗粒统计发现垂直于振动方向的近邻颗粒数增多9；2022 年 X 射线显微 CT 在三维上确认带振动记忆的浆体存在颗粒排列各向异性10。独立组的 Akiba 和 Shima 2019 年给出更强的因果性证据：同一 CaCO₃ 体系（体积分数 25%、加盐屏蔽颗粒电荷），振动的最大速度越过约 0.3 m/s 时，裂纹方向从垂直翻转为平行；相图按"最大速度 = 振幅×角频率（$V_{\max}=2\pi rf$）"组织得整整齐齐，换成按加速度组织就出现错点。他们把高速区归因于振荡流驱动颗粒重取向、连成沿流向的链，链间缝隙是力学弱面。4 颗粒相互作用也直接参与写入：带电胶体浆体记不住流动方向，加盐屏蔽库仑斥力后流动记忆恢复。12
边界干燥梯度主导——被两个事实排除：层状方向不随容器形状改变；切掉上三分之一表层不改变图样。特氟龙对照进一步说明边界只决定"裂不裂"（粘不住底就均匀收缩、一根裂纹都没有），不决定"往哪裂"。1
两个幸存载体并非互斥：残余应力是连续介质语言，颗粒排列是颗粒语言，干燥读出时两者终究都落在"哪根键先断"上。2018 年综述把"将颗粒各向异性纳入残余张力理论"列为待解决的挑战。3 对本文的问题而言，重要的不是给谁颁奖，而是记住：写入端是两个都被测到、尚未分离的载体；读出端是两者共同的终点——方向性的弱键分布。读出环节可以被单独检验，这正是下一节复算的位置。
读出端：断裂级联把方向性差异读进断键分布
写入之后剩下最后一个问题：即使内部各向异性微弱，裂纹为什么能把它读出来？本库用最标准的断裂最小模型——三角格子弹簧（熔断）网络——单独检验这个读出环节。节点间以弹簧相连、断裂阈值带 ±6% 的均匀噪声；写入步骤简化为把沿 x 方向的键弱化 w（两个候选载体的共同终点都是"某个方向的键变弱/先受力"）；读出步骤模拟干燥：边界固定（对应浆体粘住盘底），所有键的自然长度同步收缩，事件驱动地逐根断开超限的键，脱开边界的碎片释放应力退出。带结构无序的熔断/离散元网络用于收缩层裂纹已有先例15；这里只核一条定性结构：方向性差异如何被级联放大，不模拟任何具体材料。
数字本身记在这里（完整数据见 research/assets/2026-09-07-desiccation-crack-memory/desiccation_crack_results.json）：主对照（种子 7）早期断键 x 组占比从无写入的 0.49 升到写入 15% 的 0.59，全部断键占比从 0.31 到 0.40。w 扫描的早期占比增量：w=0.02 → −0.03±0.06，0.05 → +0.01±0.03，0.10 → +0.02±0.08（均淹没在种子间噪声里），0.20 → +0.07±0.06，0.40 → +0.11±0.06，0.60 → +0.22±0.06（三者 6 个种子全部为正）；同一运行全部断键口径的增量始终不超过 +0.04 且不随 w 稳定增长。在 ±6% 的阈值噪声与 6 个种子的样本量下，本模型要到约 20% 的键强差异才出现全种子一致的读出——这只描述本模型的行为，不是对真实浆体的定量预言。
这个模型只覆盖"读出"：写入的流变过程（塑性流动如何产生各向异性）不在模型内，那是 Ooshida 连续理论的地盘；实验上应力差在干燥过程中才被放大3 8，模型也不模拟这个增长，只把各向异性当作既有条件。模型的声明恰好是 Nakahara 猜想的可操作版本——即使微观差异"微弱到看不见"，断裂级联对弱键的忠实选择也会把它读进断键的方向分布；至于断键如何进一步连成宏观层状图样，本模型不检验。1
历史什么时候被洗掉
把实验相图、连续理论和读出模型合在一起，记忆的存亡由三道闸门决定，而"洗掉"本身已有直接实验：1 2 3
写入闸：必须发生塑性剪切，但又不能发展成宏观流动。低于屈服应力，浆体是弹性体，卸载即复原，什么也留不下；远超屈服应力，表面波与湍流出现，宏观流动把刚写入的微观结构搅散——在 2005 年报告的体系与参数区间内，浆体一旦整体流动就留不住方向记忆。1 窗口之外近年还有新发现：在另一套条件下（同种 CaCO₃ 但体积分数降到 25%、加盐屏蔽电荷），振动最大速度越过约 0.3 m/s 时写入的是翻转的平行（流动）记忆4——这不是 2005 年相图窗口的直接延伸，而是不同体系参数下的结果；流变仪上的振荡剪切实验则把两类记忆的转变判据定为振荡应变幅，并发现足够大的剪切不需要振荡就能写入流动记忆。13
保存闸：冻结后不再受扰——但受扰就改写。模型里塑性阈值让低于屈服的形变无限期冻结；实验上，密闭存放至少一个月后记忆仍在1——冻结机制与这一耐久观察相容，但实验只覆盖了该存放条件与时长。反过来，重新越过屈服的扰动预期会动这份记忆，在已检验的协议下改写快得惊人：偏转 45° 追加振动，两次就足以全场重写3；另有会议论文报告超声照射可擦除已写入的记忆（我们只核到书目，引为方向性证据，留作后续阅读线索）。11 添加多糖等化学物质也能增强或消除记忆效应——写入条件不只由力学决定，也由颗粒间作用决定。14 12
读出闸：浆体必须粘住边界、干燥收缩必须被约束。特氟龙对照显示粘附弱时浆体均匀收缩、一根裂纹都没有——此时没有裂纹可用于读出；仅凭这一对照无法判断内部记忆是否仍在。裂纹是读出器，读出失败不等于记忆被抹除。1 读出也不是纯粹的"显影"：振动刚停时应力差很小，干燥中才逐渐增大3 8——干燥过程自身参与了记忆的放大，这个时程目前只有测量事实、没有闭合的解释。
还有一类"换载体"的变体：把球形 CaCO₃ 换成盘状颗粒（直径 1.5 µm、厚 0.2 µm）的碱式碳酸镁，随体积分数降低，层状裂纹从垂直振动方向翻转为平行——同一振动写入了另一种记忆。5 Ooshida 称之为 Type-II，解释是盘状颗粒像液晶一样被层流取向。2 Akiba–Shima 的结果说明甚至不用换材料：同种 CaCO₃ 浆体在更低的体积分数（加盐屏蔽电荷）和更高的流速下，方向关系整个翻转。4 "裂纹垂直于振动方向"不是普适结论，而是"颗粒性质 + 塑性剪切条件"这对载体的输出；改变任一者，写入机制与读出方向都要重审。
缺口
- 两种载体的因果分离仍缺实验。残余应力各向异性（弹性底板反演8）与颗粒排列各向异性（近邻统计9、三维 μCT10）都被直接测到，但由同一次振动同时产生；2018 年综述明言把颗粒效应纳入残余张力理论是未完成的挑战。3 我们检索未见用可控手段单独改变其中一者（如固定颗粒排列只改残余应力）来裁决的实验。
- 时程问题没有闭合。振动刚停时应力差很小、干燥中才增大3 8：静置时记忆以什么形态存放、干燥中是什么机制把它放大，目前只有测量事实。这直接关系到"残余应力是种子还是全部"。
- 三维理论与三维微结构之间未对接。三维残余张力理论已有7，μCT 三维颗粒排列也有10，但两者的定量对照未见。
- 振荡写入与稳态流写入的定量等价性未建立。坡面流实验在原文中只是初步结果2；两类记忆（垂直/平行）的转变判据刚在振荡剪切下建立13，两种驱动产生的残余各向异性是否定量一致，没有系统的对照数据。
- 外场记忆与机械记忆的关系未决。电场、磁场也能定向胶体裂纹（如 Laponite 体系的电场记忆，2018 年综述 §1 有谱系3），与机械扰动写入的记忆是否同一机制，我们没有检索到定论。
小结
湿浆记得干燥前的扰动，靠的是三段接力。写入靠塑性：振动把浆体剪过屈服点，塑性流动在流向积累张力（连续介质语言）或颗粒排列各向异性（颗粒语言）——两种载体都已被直接测到，因果贡献尚未分离。保存靠冻结：模型里屈服阈值让残余形变在卸载后长期留存，实验上密闭一个月记忆仍在；但新的塑性扰动可以快速改写它，在已检验的协议下两次追加振动即全场重写，化学添加也能增强或消除。读出靠断裂：干燥收缩被盘底约束，断裂级联在早期忠实地先断弱键——本库复算显示，方向性的键强差异会被读进早期断键的方向分布，而实验上应力差本身也在干燥中才被放大。历史的存亡因此分三层：没写入（没过屈服）、被抹掉（宏观流动冲掉、新扰动改写）、以及读不出（无约束就无裂纹——读出失败不等于记忆消失）。"裂纹垂直于振动方向"从不是普适定律，而是颗粒性质与剪切条件共同选出的输出。1 2 3 4