🔬 condensed matter

Crystal to liquid cross-over for active particles with inverse-square power-law interaction

本文研究了处于谐振势阱中具有平方反比排斥力的一维奔跑-翻转粒子系统，揭示了由活性驱动的从具有尖锐密度峰值的类晶体态到类液体维格纳半圆轮廓的显著交叉转变，并得到了解析协方差计算与数值模拟的共同支持。

原作者： Saikat Santra, Leo Touzo, Chandan Dasgupta, Abhishek Dhar, Suman Dutta, Anupam Kundu, Pierre Le Doussal, Gregory Schehr, Prashant Singh

发布于 2026-08-18

📖 1 分钟阅读☕ 轻松阅读

CC BY 4.0

原作者： Saikat Santra, Leo Touzo, Chandan Dasgupta, Abhishek Dhar, Suman Dutta, Anupam Kundu, Pierre Le Doussal, Gregory Schehr, Prashant Singh

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一个这样的世界：微小的物体并非由无形的的热扰动推动，而是由它们自身的内部引擎驱动。这些是活性粒子（active particles），一类存在于从游泳的细菌到自驱动机器人等各种事物中的物质。不同于那些漫无目的地漂浮的尘埃，这些粒子消耗能量以在特定方向上移动一段时间，然后随机改变航向。这种持续的、自我驱动的运动打破了通常的平衡规则，创造出了一些科学家才刚刚开始理解的奇异且复杂的行为。当许多这类活性粒子被挤压在一起并受到限制时，它们并不会简单地混合成均匀的汤状，而是可以组织成截然不同的模式，在类固体结构和类流体流动之间切换，其方式是普通被动物质从未有过的。

一支研究团队现在已经精确地绘制出了这种转变是如何发生的，是在一个特定的、受控的环境中进行的。他们研究了一维线性的粒子，这些粒子被困在一个碗状的力场中，其中每个粒子都会推开邻居，且这种斥力随着距离缩短而增强。至关重要的是，这些粒子被建模为“奔跑与翻滚”（run-and-tumble）型代理，这意味着它们以恒定的速度向前冲刺，然后随机翻滚以选择新的方向。通过运行大量的计算机模拟并开发新的数学工具，该团队发现，随着这些粒子速度的增加，整个系统经历了一个剧烈且可预测的演变过程。它从一个刚性晶体开始，融化成一种平滑的液体，最后拉伸成一个宽阔的钟形云团。

在开始阶段，当粒子运动非常缓慢时，它们表现得像一个冻结的晶体。因为它们之间的排斥力如此之强，它们相对于彼此锁定在精确、固定的位置上。如果你观察一段时间内它们所在的位置，你会看到密度分布中有尖锐、清晰的峰值，每个峰值代表一个守住阵地的粒子。研究人员发现，这些位置对应于某种特定类型多项式的数学零点，这种模式在被动物理学中已被研究了几十年，但在这里却是由主动运动产生的。只要粒子的内部能量较低，它们只是在这些固定点附近轻微晃动，保持着有序的、类固体的结构。

随着研究人员增加粒子的速度，系统开始融化。粒子开始远离它们原始的位置，它们的运动也开始发生重叠。密度分布的尖锐峰值变得平滑，融合在一起，直到单个粒子的位置无法再被区分出来。在这个中间阶段，系统看起来像是一种液体。值得注意的是，这个液体云团的整体形状与一个著名的数学曲线——维格纳半圆（Wigner semi-circle）相匹配。这种形状通常与处于热平衡状态下的被动系统相关，但在这种由主动运动驱动、远离平衡态的系统中，它却出现了。研究人员通过测量粒子波动程度相对于它们之间距离的大小，确认了这一转变。当这些波动变得大到足以模糊邻里之间的间隙时，晶体便有效地转化成了液体。

然而，故事并未到此结束。当粒子被推向更快的速度时，平滑的半圆开始发生畸变。密度分布被拉伸，失去了其弯曲的顶部，并呈现出一种宽阔的、钟形的形态，其延伸范围远超液体阶段。在这种高度活跃的状态下，粒子不再仅仅是在原地抖动；它们正充满活力地探索整个陷阱，以至于云团的边缘是由它们的最大速度而非力的平衡来定义的。研究人员观察到，这个云团边缘处的密度并非骤降，而是遵循特定的幂律（power law）逐渐消散，这种模式在粒子接触时经历无限斥力的此类系统中似乎是一个普遍特征。

为了理解为什么会发生这些变化，团队开发了一个理论框架，其核心思想是粒子正在对其理想的、冻结的位置进行微小的偏离。通过计算这些偏差如何相互关联，他们可以预测粒子位置的方差（即离散程度）。他们的计算表明，从晶体到液体的转变发生在粒子速度达到某个阈值时，且该阈值与系统中粒子的数量无关。第二个转变，即从液体到钟形云团的转变，仅在速度高到与粒子总数成比例时才会发生。这种分析方法与他们的计算机模拟完美契合，证实了从有序到无序、再到一种新型主动无序的转变是由简单、可预测的规则所支配的。

这些发现为活性物质在受限条件下如何自我组织提供了清晰的图景。研究人员表明，即使在一个远离平衡态的系统中，只要活性不是过于极端，粒子也可以稳定在一种模仿被动系统平滑、类液体行为的状态中。但一旦活性跨越了一个临界极限，系统就会脱离这些熟悉的形状，采用一种由粒子巨大动量决定的新形式。这项工作不仅阐明了陷阱中活性粒子的行为，也暗示了类似的转变可能在现实世界的实验中被观察到，例如在细菌或合成微型机器人实验中，前提是它们能够被限制且其相互作用可以得到控制。这项研究证明了将精确模拟与严谨理论相结合的力量，从而揭示复杂、自我驱动系统背后的隐藏秩序。

技术摘要：具有反平方幂律相互作用的活性粒子的晶体到液体跨越

问题陈述
本研究调查了由 $N$ 个在谐振阱中受限的运行-转向粒子（Run-and-Tumble Particles, RTPs）组成的单维系统的非平衡稳态性质。这些粒子通过排斥性的反平方幂律势（ $V(x) \propto 1/x^2$ ）进行相互作用，该模型在哈密顿动力学中被称为 Calogero-Moser 模型或反平方模型（ISM）。主要目标是理解引入活性（自驱动）如何改变系统的密度分布和结构有序性，特别是考察从晶体态到类液体态，并最终到一种独特的强活性机制的跨越过程。这项工作扩展了先前关于活性 Dyson 布朗运动（DBM，具有对数相互作用）的解析研究，旨在探索在非平衡驱动下，活性 ISM 与这两个可积系统之间的联系。

方法论
作者采用结合广泛数值模拟和弱噪声极限下解析推导的双重方法。

数值模拟： 系统根据具有反平方排斥力和谐约束势的 RTPs 的类 Langevin 方程进行演化。稳态密度分布 $\rho(x)$ 通过对粒子位置进行时间平均和多样本实现计算得出。关键观测指标包括：
- 全局密度分布及其与系统规模 $N$ 的标度关系。
- Lindemann 指数（ $\eta_i$ ），定义为粒子位置的标准差与平均粒子间距之比，用于量化晶体到液体的转变。
- 距离度量（ $\chi$ ），用于量化稳态密度与 Wigner 半圆（WSc）分布的偏差，以识别向强活性机制的转变。
- 粒子间隙的方差。
解析方法： 在小活性（ $v_0 \to 0$ ）和大持久时间（ $\gamma \to 0$ ）的极限下，作者围绕平衡构型（Hermite 多项式的零点）对运动方程进行线性化。他们利用了 ISM 势的 Hessian 矩阵 $H_{ISM}$ 。一个关键的理论洞察是利用关系式 $H_{ISM} = H_L^2$ ，其中 $H_L$ 是 DBM 的 Hessian 矩阵。通过使用 Hermite 多项式及其相关特征向量对 $H_L$ 进行对角化，作者推导出了粒子位移和粒子间隙协方差及方差的闭合形式表达式。随后，利用傅里叶空间技术，将这些结果扩展到有限的转向率 $\gamma$ 。

主要结果

三个不同的活性机制： 研究根据活性参数 $v_0$ 确定了三个不同的机制：
- 弱活性机制 ( $v_0 \ll 1$ )： 系统表现得像晶体。密度分布在平衡位置（Hermite 零点）中心呈现出尖锐的峰值。涨落相对于粒子间距非常微小。
- 中等活性机制 ( $v_0 \sim O(1)$ )： 随着活性的增加，粒子涨落变得与粒子间距相当。尖锐的峰值消失，密度分布在体相部分平滑为 Wigner 半圆（WSc）形式。这种“液体”相在较宽的活性范围内持续存在。
- 强活性机制 ( $v_0 \sim O(N)$ )： 在极高活性下，密度分布显著偏离 WSc，转化为一个支撑集延伸至 $\sim \pm v_0$ 的钟形分布。值得注意的是，该分布的边缘表现出普遍的幂律衰减 $\sim 1/x^3$ 。
定量跨越准则：
- 从晶体（有峰值）到液体（平滑 WSc）相的转变发生在 Lindemann 指数 $\eta_i$ 超过阈值（约 0.5）时。解析结果表明，这种跨越发生在 $v_0 \sim O(1)$ 时，且与系统规模 $N$ 无关。
- 从液体（WSc）到强活性（钟形）相的转变由距离度量 $\chi$ 超过 0.1 来标记。这发生在 $v_0 \sim O(N)$ 时。
涨落的解析预测：
- 作者推导了在大 $N$ 极限下粒子位置方差 $\langle \delta x_i^2 \rangle$ 和间隙方差 $g_{i,n}^2$ 的精确表达式。
- 与活性 DBM 不同（在 DBM 中，边缘粒子表现出与体相粒子不同的标度行为），活性 ISM 显示对于体相和边缘粒子，方差均按 $1/N$ 衰减。这归因于涉及 $1/k^4$ 因子的级数展开（源于 $H_{ISM} = H_L^2$ 关系）比 DBM 中的 $1/k^2$ 因子收敛得更快。
- 即使对于有限的 $\gamma$ 以及高达 $v_0 \sim O(N)$ 的情况，陷阱中心附近的间隙方差解析预测也与数值模拟高度吻合。
边缘行为与翼部：
- 在中等机制中，活性 ISM 不会出现活性 DBM 中观察到的“翼部”（延伸尾部）。作者将其归因于在 ISM 中，边缘粒子的相互作用力由最近邻支配（标度为 $\sqrt{N}$ ），从而有效地钉扎了它们；而在 DBM 中，力由集体体相支配（标度同样为 $\sqrt{N}$ 但系数不同），从而允许更大的涨落。
- 然而，在小 $\gamma$ 的强活性机制下，由于驱动速度相对于谐振阱占主导地位，边缘附近可能会出现“翼部”结构。

意义与主张
本文声称对活性 ISM 进行了全面的表征，建立了与活性 DBM 的直接平行关系，同时强调了两者在涨落标度和边缘稳定性方面的关键差异。

Wigner 半圆的普遍性： 研究强化了这样一个发现：即 Wigner 半圆是谐振阱中相互作用活性粒子的鲁棒密度分布，即使在高活性水平下也能保持，直到发生第二次向钟形分布的转变。
解析可处理性： 通过利用 ISM 与 DBM Hessian 之间的联系，作者成功推导出了在广泛参数范围内（包括“液体”相）依然准确的粒子位置和间隙方差公式。
相图： 通过构建 $(a, 1/\gamma)$ 平面（其中 $v_0 = aN$ ）的相图，定量划定了晶体、液体和强活性相之间的边界。
实验相关性： 作者指出，如果能够设计出有效的幂律相互作用，所预测的晶体到液体熔融转变以及特定的密度分布，可以在当前的受限活性粒子系统（如细菌、Janus 粒子）中被实验观测到，尽管他们也指出在这些实验中相互作用效应尚未得到充分探索。

研究结论认为，虽然活性 ISM 和活性 DBM 具有定性的相似之处（三种机制、中间 WSc 相），但 ISM 特有的幂律相互作用导致了截然不同的定量行为（如涨落和边缘稳定性），这为理解可积系统中的非平衡统计力学提供了更深入的见解。

类似论文