弹珠和陀螺的游戏—制造“超均匀”智能流体材料

无序超均匀态

物质一般可分为有序态和无序态。有序态在空间上分布均匀齐整，比如晶体以及准晶（见图1B）。而在无序态中，粒子的位置有很强的随机性，造成了分布的疏密不均（见图1A）。无序态涵盖了气态、液体、玻璃态等物态。2003年，美国普林斯顿大学化学系S. Torquato教授，提出了介于无序和有序之间的一种新的物态，即无序超均匀态。在小尺度上观察，这种结构显得毫无规律，但在大尺度上却具有晶体相似的均匀性（见图1C）。

在大自然以及数学世界中，科学家都发现了无序超均匀态的踪影，例如鸟类视网膜上的视锥细胞分布、质数的分布、随机矩阵的本征值分部、硬球的密堆积、宇宙大尺度结构等。无序超均匀态也具有很多奇特性质和用途，例如，它能像晶体一样产生光子带隙，并且对某些波长的光是完全通透的（超透明性）。不仅如此，这种分布也在数值抽样和计算机图形学领域运用广泛，比如，具有超均匀分布的伪随机数可以极大加速蒙特卡洛抽样的速度（随机打点法计算π值），同时用超均匀的无序点阵进行抽样绘图不仅可以克服图形失真（“反锯齿”），还可以用最少的像素点画出最佳的图像（见图2）[1,2]。

 图1：粒子分部的均匀性可由窗口内粒子数N的涨落（方差） 的标度率来反映（d为维度）。λ=0意味着体系的粒子数涨落正比于粒子数本身，即泊松过程或随机分布（图A）；而 λ=1则意味体系有着晶体般的超均匀性（图B）。0＜λ≤1都属于超均匀态。（图片编辑自[1]）

超均匀流体 

目前无序超均匀态的研究主要集中在无序固体。今年年初，新加坡南洋理工大学的一个理论研究小组首次提出了超均匀流体的概念，并且利用了自驱动胶体体系模拟验证了这种特殊流体的存在（见链接[3]）。但是这种特殊流体的形成机制，以及其与平衡态流体的本质区别仍然是一个未解之谜。与此同时，理论上也缺少描述这种奇异流体的流体力学方程。最近该小组的雷群利博士和倪冉教授对这个问题进行了深入的探索，他们通过构建一种简单的“弹珠”模型，巧妙地实现了从简单流体到超均匀流体的平滑过度，不仅如此，他们还发现超均匀流体在理论上具有非常简洁优美的流体力学方程形式。基于这套理论，超均匀流体的形成机制可以用简单的谐振子模型来理解。最后，作者还通过模拟展示了如何利用简单自驱动“陀螺”制造这种超均匀流体。相关工作于近期发表于《美国科学院院刊》上。

图2：超均匀点阵可用于消除图的水状波纹（上）和图形视觉优化（下）右下角是这种点阵的结构因子S（q）。（图片来自 [4]）

“弹珠”模型

作者研究的是一种非平衡态的硬球体系（图3）。这个模型类似于普通弹珠在有摩擦的平面上运动，即硬球会受到和速度大小成正比的阻尼（摩擦）力，并且两球碰撞满足动量守恒。该体系和弹珠碰撞不同的地方是，在碰撞过程中，会有额外的能量输入给相碰的两个硬球，使两个硬球在碰撞后的总动力增加？E，这被称作“激发碰撞”（activecollision）。在这个体系中激发碰撞提供能量输入，阻尼力又把能量耗散掉。这种“驱动-耗散”（drivendissipation）机制是活力物质体系的基本特征。体系的最终运动状态还是取决于两个特征长度。一个是耗散距离，即在？E的激发下，静止硬球可以运动的最远距离。如果我们固定？E 不变，那么这个距离完全由摩擦系数决定。所以，其实反映了硬球在阻尼系统中“惯性”的大小。另一个特征长度是“平均自由程”  （mean free path）即硬球在两次碰撞之间平均的运动距离。这个特征长度反映的是硬球体系疏密程度，越稀的体系，硬球之间越不容易碰撞，平均自由程就越大。

作者分析发现当在初始时刻赋予粒子随机速度后，如果   ，体系的碰撞事件数会随着时间指数减小，最终体系陷入所谓的“吸收态”（absorbingstate），即每个硬球的速度衰减为零。相反，如果，体系的激发碰撞会产生“链式反应”，使体系保持在一种持续碰撞但“温度”恒定的状态，即活力态（active state）。这里的“温度”恒定指的是体系的平均动能不变。作者发现，这种活力态就是一种超均匀态流体，其结构因子满足标度率S（q）~ q2（见图4A）。这里q 代表的是波矢量，S（q）反映的是体系在波长为2π/q 的尺度上密度涨落的大小。有意思的是，随着硬球惯性（）逐渐变大并趋向无穷（对应无摩擦），体系会升温，并且S（q）在小波矢区域会慢慢抬升（图4A），最后转变为普通平衡态流体的标度率，即S（q）~ q0（const）。硬球的速度分布也从非玻尔兹曼分布，转变为玻尔兹曼分布（图4F）。也就是说，通过调节摩擦力或者？E，我们可以使平衡态硬球流体平滑过度到非平衡态超均匀流体……

 图3：（A）研究模型。（B）二维体系的超均匀流体 （不同颜色代表不同速度方向）。

超均匀流体的流体力学方程和流体动力学 

为了研究这种奇特流体的流体力学机理，通过分析推导，作者发现了这种超均匀流体满足如下的流体力学方程：

这个流体力学方程的形式非常简洁：其中第一个方程代表粒子数守恒，第二个方程代表动量守恒。由于有阻尼力的存在（绿框），体系动量实际上是不守恒的。激发碰撞反映在红框的噪声项里。这套理论的预测和模拟结果高度吻合（图4的虚线），佐证了这套流体力学的正确性。动力学上，这套理论成功预测了超均匀流体的涨落形式不仅可以是扩散模（diffusive mode）也可以是声学模（acoustic mode）（图4G）。体系的动力学相图在图4E里给出。另外，作者还测量了超均匀流体的密度涨落的“色噪声”，发现色噪声的谱函数在低频满足特殊的f1/2的标度率（图 4B），这个发现有助于在频率空间探测和研究这种超均匀流体。

图4：（A）不同惯性（）大小的硬球体系的结构因子S（q），其中q 为波矢量，q2代表超均匀标度率；（B）密度涨落的谱函数 f 随频率的变化，其中f -1/2为超均匀标度率。；（C，D）平衡态流体和超均匀流体的谐振子模型。（F）超均匀流体中粒子速度分布随着的增加慢慢过度到玻尔兹曼分布。（E）体系的动力学相图。（G）动力学结构因子S（ω，q），其中三峰曲线的左右两侧峰是Brillouin峰，代表了声学模，蓝线的单峰代表了扩散模。注：虚线为理论预测。

超均匀流体的机理解释—谐振子模型

我们知道，根据傅里叶分解，任何形状的方程都可以看作为无数不同波长的正弦波的叠加。同样地，密度涨落也可以看成是众多的疏密相间的正弦波涨落（密度模）的叠加，结构因子 或者就是代表了在波长2π/q的尺度上的正弦波振幅。通过理论推导，作者发现满足一个非常熟悉的二阶动力学方程：

这个方程与热浴中谐振子的动力学方程完全一样，其中可以看成谐振子的“位移”q-2可以看成谐振子的“质量”，右边第一项括号里可以看成谐振子的阻尼系数，右边第二项是回复力（正比于位移），最后一项是热浴的噪声。也就是说，流体在某个波矢量q的密度涨落可以看成是一种抽象的谐振子振动（图4C,D）。我们知道，谐振子主要有两种运动模式：一种是在阻尼系数比较小时的“共振态”（图4C），另一种是在阻尼系数比较大时的“过阻尼态”（图4D）。共振态下，流体中疏密相间的密度模会表现出一种“弹性”并来回震荡，形成向前或者向后传播的“声波”；在“过阻尼态”下，密度模失去了弹性，产生之后就原地衰减，像冰淇淋融化一样慢慢铺平，形成一种“扩散波”。对于平衡态流体，在谐振子的阻尼系数（）和在热浴温度不变的情况下，根据能量均分定理，谐振子的振幅是不变的。也就是说，对于理想平衡态流体，不同波长下的密度涨落幅度是一样的，即S（q）~const. 这和图3A所示的结果完全一致的。然而，对于超均匀流体，我们发现随着波矢q慢慢减小到0，谐振子的阻尼系数会变为无穷大，导致谐振子的振幅为零。也就是说，密度涨落在无限长波长下为零，这正是超均匀态的定义。这个简单而具有启发性的模型同时也可以解释之前所述的超均匀流体的动力学性质。

“陀螺”模型

根据上面的理论分析，作者最后提出了一个实际模型来实现超均匀流体态，即自驱动的“陀螺”（spinner）集群。在外力矩（比如旋转磁场）的驱动下，陀螺会自发高速旋转。当两个陀螺碰撞时，陀螺的旋转动能会转变为平移动能，使两个陀螺相互弹开，从而实现“弹珠”模型中的“激发碰撞”（见图5A）。在有摩擦的条件下，作者发现这种“陀螺”集群会呈现出与之前“弹珠”模型一样的超均匀标度（图5B），并且体系中密度涨落的谱函数的标度率也和“弹珠”模型完全相同。这个实际体系证明了作者理论中所揭示超均匀流体机制的普适性。不仅如此，通过文献调研，作者也在“微型陀螺”的实验体系发现了超均匀流体存在的证据。

图5:（A）二维底板上，由外力矩驱动下高速旋转的陀螺满足“碰撞激发”的条件。（B）结构因子S（q）显示陀螺集群是一种超均匀流体。（C）陀螺集群的密度涨落的谱函数和简单硬球模型图4B基本相同。

结 语

总之，作者通过提出了一个简单的“弹珠”模型研究了超均匀流体的形成机理，发现只要满足相互的“碰撞激发”和“摩擦耗散”这两个条件就可以产生超均匀流体态。为了进一步地证明结论的普适性，作者还展示了如何利用自驱动的“陀螺”集群实现超均匀流体，这对实验具有极强的指导意义。在理论上作者也做出了突出贡献，不仅首次发现了超均匀流体的流体力学方程，还提出了启发性的谐振子模型来帮助理解超均匀流体态。从材料学角度上看，这种超均匀流体有望成为一种具有“自修复”和“自适应”的智能活性流体材料，比如活性“光子液体”。这种流体材料不仅具有像光子晶体一样的光学性质，而且活性“光子液体”即使受到机械损伤，也可以自我修复，同时也可随着外场驱动的变化改变自己的光学性质，是不是很科幻呢？