颗粒排列看起来杂乱，为什么在大尺度上却能异常均匀？

不排成整齐队列，也能少一些忽疏忽密

向桌面撒一把颗粒，近看会有空隙、挤堆；拿一个固定大小的圆框在桌面上移动，框里的颗粒数也会时多时少。把圆框放大后，相对波动通常会减小，这件事本身并不特殊——平均的对象多了，偶然偏差就比较容易相互抵消。

真正需要解释的是：**有些排列没有明显的晶格，却比独立随机撒点更加有效地压低大尺度密度波动。**它们局部仍有不规则间距，也没有要求每个小区域恰好一样；但不同地方的疏密偏差之间存在联系，大区域不会任由这些偏差一起越积越大。

这类空间统计称为“超均匀性”。本篇问的是其中一条形成路线：只有局部相互作用、每次动作还带着随机性，系统怎样长出这样的关联？答案要同时追踪两件事：重排怎样平整原有的密度起伏，以及重排本身又怎样制造新的起伏。12

2026年的一项研究把几种随机组织模型连接起来，指出局部噪声之间的关联，能决定活动状态下是否保留长波波动。本篇从一个更容易复算的粒子搬运模型讲通这条机制，再回看新论文与真实乳液的证据。3

“均匀”需要一把随尺度变化的尺子

在统计均匀的点分布中，记半径为$R$的球形观察窗内粒子数为$N_R$。独立泊松点的数目方差等于平均值，均随窗口体积增长：

$d$是空间维数。超均匀要求大窗极限下，方差相对体积的增长受到进一步抑制：

这里比较的是数目方差。如果改为框内平均密度的方差，还要再除以窗口体积的平方，不能把两种幂指数混着读。

另一把尺子是结构因子$S(\mathbf k)$。它把空间起伏分解成不同波长：波数$k$越小，观察的波长$2\pi/k$越大。对周期盒子里的粒子，在非零允许波数处可计算：

角括号表示适当的构型或系综平均。对于这里讨论的统计均匀体系，超均匀的频谱判据是$k\to0$时$S(k)\to0$。它说的是逐渐变长、但仍非零的空间模式，不是把一个有限盒子的总粒子数固定之后，宣布零模式没有波动。1

若$S(k)\sim k^\alpha$，球窗数目方差有三类常见渐近行为：$0<\alpha<1$时约为$R^{d-\alpha}$，$\alpha=1$时约为$R^{d-1}\log R$，$\alpha>1$时通常由$R^{d-1}$的边界项主导。因此，$S\sim k^2$不意味着数目方差是$R^{d-2}$。在一维，对应的强超均匀数目方差可以趋于常数；在二维则随圆窗周长增长。1

理想晶格也可以超均匀，但这项指标没有要求周期性。反过来，局部排斥、不许颗粒重叠或者接触网络能够承重，也没有单独保证这个长波极限成立。排列方式、力学刚性与密度波动，是不同的问题。

一排小格子：平均都在摊平，为什么最后不一样

先把真实颗粒简化成一圈格点。每个格点允许放多个粒子，格点间距取为1。如果第$i$格有$n_i>1$个粒子，它以$n_i-1$的速率发生一次搬运，每次拿出两个粒子。

2017年Hexner与Levine研究的COMCon模型，让这两个粒子一个去左邻居、一个去右邻居。粒子总数不变，局部质心也不变：两个粒子原来都在$i$，现在分别在$i-1$和$i+1$，位置和仍是$2i$。在环上，跨边界时相应守恒按周期位置理解。2

本篇另设一个对照：以概率$\delta$，让两个粒子各自独立选左或右；其余时候仍一左一右。

• $\delta=0$：每次严格一左一右。
• $\delta=0.01$：1%的事件改为独立选向；其中仍有一半碰巧一左一右。
• $\delta=1$：每次都独立选向，可能一起去左边，也可能一起去右边。

三者在同一瞬时构型下，对每个邻居的期望搬运量相同，使用的事件速率也相同。若暂时忽略很低占据数处的截止，平均密度都满足同一个离散扩散方程：

也就是平均而言，较密的地方向较疏的地方输送粒子。只看这条平均方程，会以为三个模型最后没有重要差别。但一起向左的事件还会制造局部位置偏移，严格一左一右的事件消除了这一层随机偏移。同样的平均流，不代表同样的涨落。

本库用连续时间事件模拟实际执行了这些规则。模型在每一步严格保证粒子数非负、总数不变；不是把线性公式直接画成“模拟”。环长取64与128，平均每格20粒子，每组4个随机种子，在充分等待后取样。这个高密度选择是为了比较可解析机制，不代表真实乳液的体积分数，也不研究临界点。模型预设了格点，研究的是其上的占据数涨落；它并没有模拟无格点材料从头形成无序结构。

局部配对，怎样削弱新制造的长波波动

看一个波数为$k$的密度模式。单个粒子搬一步，它对模式的改变量在小$k$时从一阶项开始。两个粒子从同一格分别左右搬一步，改变量则为：

左右两边的一阶项抵消了。于是，这种事件对长波模式的随机冲击更弱：振幅从$k^2$开始，其平方从$k^4$开始。与此同时，扩散对原有模式的衰减率仍从$k^2$开始。

这给出一条很具体的机制链：局部配对约束改变随机流的低波数结构；密度平整仍然进行；新噪声补回大尺度起伏的能力变弱；稳态的长波密度方差随之下降。它不需要某个粒子知道全局密度，也不要求所有随机性消失。

令$\lambda_k=2(1-\cos k)$。把所有可能的左右搬运逐项枚举，本篇得到单次事件的二阶矩：

其中，独立选向带回了一个在小$k$处更占优势的$k^2$噪声强度项。高密度近似下，以总粒子数归一化的稳态结构因子为：

$\rho$是每格平均粒子数。$\delta=0$时，$S\sim k^2$；只要$\delta>0$，这个模型的大尺度极限就留下$(\rho-1)\delta/\rho$的平台。平台可以非常低，在有限尺度上很难发现，但低平台与趋于零仍是不同情况。

这里的推导保留了COMCon的高密度条件；真实跳跃模拟在$n_i\le1$时停用事件，不会允许负粒子数。采样中这种低占据情况极少但并非所有组都完全没有，因此不能把近似式宣称为任意密度下的严格解。2

在$L=128$的最小非零波数处，三组$S$分别约为0.00111、0.00907和0.946；解析预测为0.00114、0.01064和0.951。低波数处独立信息较少，数值并非逐点精确吻合，但三种规则的差别远大于这组采样误差。

直接看32格窗口的粒子数方差，三组约为19.06、22.98和462.85。严格配对组的局部粒子仍不断活动，大窗口的数目波动却没有像普通随机搬运那样不断累积。这是空间关联，不能从“每个粒子都动得很随机”直接猜出来。

只锁住整个系统的质心，够不够

不够。局部配对与一个全局约束的作用不同。

假设先让所有粒子按某种规则移动，再从每个粒子的位移中减去同一个平均位移，使整体质心不动。这一步只是把整幅新构型平移了一下。任何非零周期波数的傅里叶振幅只会乘上一个相位，绝对值不变，所以结构因子完全不变。

本篇实际计算也核对了这点。因而，“每一步总质心不变”若只是通过整体平移实现，不能把原本普通的空间波动变成超均匀。真正改变前述噪声谱的，是近邻搬运在事件发生的位置就相互配对。

这同样不是说局部质心守恒是所有超均匀结构的唯一形成条件。长程相互作用、特定排列约束、某些临界组织过程，可以通过别的机制抑制长波涨落。本篇的结论限定在所分析的活动搬运路线。12

2026年的进展：把几种不同的随机更新放到同一把尺子下

Anand、Zhang与Martiniani比较了三类相互作用粒子模型：随机组织RO、沿颗粒连心线偏向排斥的BRO，以及对粒子相互作用能量做随机梯度更新的SGD。它们的随机性分别涉及位移方向、位移大小或本轮选中了哪些相互作用；不是同一个算法换名字。3

这项工作的重点是活动相：颗粒一直在重排，形成统计稳态。作者从成对相关的微观噪声出发建立连续模型，再推到有涨落的密度方程；比较了离散粒子、连续时间粒子及连续密度模型。

在足够长波、密度线性化有效的条件下，归一化结构因子具有：

$c$是相应微观噪声的成对相关参数，完全反相关时$c=-1$；$M$仍含模型细节。它与本篇的事件混合概率$\delta$不是同一个参数，不直接设成$c=-1+\delta$。34

这个式子的第一项留下长波平台，第二项给出二次波数贡献。$c=-1$时平台消失；对$c>-1$，谱中两项相等的位置定义一个归一化交叉长度：

例如RO的$M=1/2$。当$c=-0.99$，该长度约为$9.97\ell_0$，其中$\ell_0=2\pi/k_0$是这套归一化里的参考尺度。看比它短的一段尺度，可能像二次幂；再看更长的尺度，会进入低平台。

本篇还从附录保留完整空间核的线性公式重算了二维RO：用圆形作用区的傅里叶核$F(k)=2J_1(ak)/(ak)$，得到$\widetilde S=2(1+cF)/(1+F)$，再与上面的二阶展开对照。$a$是作用截止距离。图中的完整核只是在线性理论内部保留更多波数信息，不是重新运行作者的三类大规模粒子仿真。4

新论文中的共同形式很有解释力，但“归一化后有共同规律”不等于真实长度、噪声幅度、形成时间都只由$c$决定。其线性理论主要适用于高密度软重叠模型，作者明确说明不能用它解释吸收相或临界点的空间统计。多体相互作用和一般神经网络也没有由此一并解决。3

结构因子与数窗，不能只把横轴取个倒数

一个波数对应一个空间波长，但一个大小固定的观察窗，并不只选中这一个波数。窗口会混合一段频谱；窗的边界还会引入独立的贡献。因此，从频谱读出的交叉长度，与从数目方差读出的交叉长度，需要各自定义。

本篇格点模型给出一个可直接核对的例子。在长度$L$、总数固定的环上，对$0<\ell<L$的连续格点窗口，高密度理论给出：

第一项随窗口大小增长，并带固定总数的有限盒修正；第二项是边界贡献。模拟图里的数目方差是直接对构型计数得到的，可以与这个频谱推导出的式子交叉比较。

在$\ell\ll L$下，让两个数窗贡献相等，得到$\ell\approx(1-\delta/2)/\delta$；而让谱中的平台项与二次项相等，再换成波长，得到$2\pi/k_c\approx2\pi\sqrt{(1-\delta/2)/(2\delta)}$。

例如$\delta=0.01$时，两种定义分别约为99.5格与44.3格。前者需要远大于100格的系统才能不受全盒修正地观察；本篇$L=128$并没有提供这样的无限空间验证。我们能确认的是数学定义不同，以及模拟在实际可用窗口内符合各自的有限尺寸预测。

在严格配对的高密度近似里，还能把抵消关系写得更直观：一个格点的占据数方差为$\rho-1$，它与左右相邻格点的协方差各为$-(\rho-1)/2$，更远的两点协方差在这层近似中为零。把一段连续格点加起来，内部的正负贡献相互抵消，留下边界贡献。于是，大尺度统计很特殊，并不要求两点协方差一定有缓慢衰减的长距离尾部。

这也解释了为什么不能从$S\sim k^2$直接推出数目方差的任意尺度幂律。2026论文图中的数密度方差由测得的$S(k)$积分得到，属于同一数据的另一种表示；不能把两张图当成两次独立实验。3

大尺度的均匀，不会瞬间长出来

局部规则每次只能搬动很近的粒子。更大范围的密度起伏需要许多次传递才会松弛。在前面的高密度模型中，模式均值的衰减时间为$1/\lambda_k$，模式方差接近稳态的时间为$1/(2\lambda_k)$；小$k$时均随$k^{-2}$增大。

如果初始点分布有普通的白色长波涨落，线性理论给出：

即使最终理论是超均匀的，有限等待时间内，足够长的波仍然主要保留初态。也就是说，先把时间等得很长再看长波，与先在固定时间追问无限长波，结果可以不同。2

本篇把系统长度从64增到128时，同时把等待和取样间隔增为四倍。严格配对组的最小波数结构因子由约0.00545降到0.00111，方向与二次波数预测一致；独立选向组仍在约0.94的平台附近。两种尺寸只是一项有限验证，并不是已经检验无限系统。

若再加入独立热扰动或其他不配对运动，长期平台也可能恢复。2026研究在其具体模型中同时加入热扩散与热噪声，得到长波涨落随温度增加；不能只改变噪声而忽略相应输运，再把自己的式子说成真实热体系。34

真实乳液看到了什么，没看到什么

一组承重实验来自2015年的周期驱动乳液。研究者追踪受限微流道中的液滴，在不同往复幅度下比较运动是否能大致回到原处，以及空间疏密如何变化。约3000个液滴留在观察视野内；小幅度下结构更均匀，接近可逆性丧失的位置，数目波动受到更强抑制。5

但论文没有证明真实液滴进入一个永久静止的理想吸收态。所谓可逆，是按液滴是否回到原来的Voronoi区域定义；该区域内仍有不可逆的小位移与持续相互作用。实验中的均匀区也有有限范围，大致达到十个液滴直径的尺度，再往大处看就出现转折。

这件事很有意义：有限范围内异常均匀的结构确实能由驱动中的物质形成，但真实实验并不因此等同于理想吸收态模型，或已证明无限尺度的超均匀。

实验配套模型又做了一个有用对照：正式发表版把流体相互作用在超过两个液滴直径时截断，仍可形成这种均匀组织。它说明在该模型里，长程作用不是这段现象的必要成分。模型调过一个偶极强度参数来匹配转变幅度，不能称完全无拟合；截断后是否发生真正临界转变，作者也保留了需要有限尺寸分析的限制。5

因此，目前几段证据各回答不同问题：格点模型讲清局部守恒如何改变噪声谱；2026研究连接几种粒子动力学的活动相；乳液实验确认受驱动物质中的有限尺度组织。它们相互启发，但不构成对同一个微观机制的三次独立复现。

两个很容易把普通随机点看成特殊秩序的操作

第一个是拿接近整个样品大小的窗口去数。

如果盒中预先固定有$N$个相互独立、均匀分布的点，一个占盒体积分数$f$的窗口，其计数是二项分布：

于是方差除以平均值是$1-f$，当窗口接近整个盒子时趋于零。这里没有通过局部作用生成任何新秩序，只是整盒总数早已固定。检测真正的大尺度关联，需要在颗粒尺度与整个样品尺度之间留出空间。

第二个是对有限窗口减去样本均值，再追问比窗口还长的波。

本篇构造5000组独立随机点，每组128个，均匀放在长度$L$的一维窗口内。对任意连续波数，先减去同一窗口的样本平均密度，再计算谱。可直接求得该估计量的期望：

这里$\operatorname{sinc}(x)=\sin x/x$。它也出现漂亮的$k^2$！原因是减均值强制了这个有限窗口的零模式为零，低波数展开随之产生二次项；底下的点却始终是独立随机的。

在正规的非零周期傅里叶波数$k=2\pi m/L$处，上式恰好回到1。错误在于把窗外尺度的连续低波数曲线，当成已经测到无限系统的空间统计。这个反例针对一种明确估计方式，不意味着所有有限窗口分析都必然产生同一种假象。

2026年也有专门讨论观察窗和遮挡对频谱影响的新预印本。它提示这是活跃的方法问题，但其中的二次展开含有常数项，不能只看摘要里的“二次标度”就说所有观测谱都趋于零。本篇的反例由自己的明确取样定义推导并模拟，不依靠该预印本的普适说法。6

哪些变化最能区分解释

如果只看到颗粒变得均匀，可以有很多解释：局部配对噪声、长程排斥、接近某个临界状态，或者取样造成的表象。下面几种比较，比单看一条好看的拟合线更有区分力。

可以带走的机制图景

局部运动并不直接决定大尺度一定更乱。决定结果的，是原有密度起伏被搬运摊平的速度，与新噪声把起伏补回去的方式。局部质心守恒提供了一条具体路线：让长波随机冲击的一阶贡献抵消，同时保留密度松弛，于是形成强烈的空间反相关。

这条主线已有可解析模型、数值研究与相关驱动实验支持；其中，本篇自己核实了高密度格点模型的机制和对照，没有独立重跑2026论文的全部大规模模型，也没有重做乳液实验。离散模型的临界点、真实多体流体作用及实验无限尺度极限，仍在本篇的覆盖范围之外。

对实际材料，“在多大范围、形成多久、哪些扰动下仍然更均匀”往往比只给一个名称更有用。也正因为如此，超均匀性值得作为一种可测量、可改变的空间组织方式来研究：把局部随机动作与整体空间秩序之间的桥梁写出来，而不把一张杂乱但好看的颗粒图当作结论。