2025-11-25T17:10:18.317869

Random organization criticality with long-range hydrodynamic interactions

Jocteur, Nardini, Bertin et al.

Driven soft athermal systems may display a reversible-irreversible transition between an absorbing, arrested state and an active phase where a steady-state dynamics sets in. A paradigmatic example consists in cyclically sheared suspensions under stroboscopic observation, for which in absence of contacts during a shear cycle particle trajectories are reversible and the stroboscopic dynamics is frozen, while contacts lead to diffusive stroboscopic motion. The Random Organization Model (ROM), which is a minimal model of the transition, shows a transition which falls into the Conserved Directed Percolation (CDP) universality class. However, the ROM ignores hydrodynamic interactions between suspended particles, which make contacts a source of long-range mechanical noise that in turn can create new contacts. Here, we generalize the ROM to include long-range interactions decaying like inverse power laws of the distance. Critical properties continuously depend on the decay exponent when it is smaller than the space dimension. Upon increasing the interaction range, the transition turns convex (that is, with an order parameter exponent $Î²> 1$), fluctuations turn from diverging to vanishing, and hyperuniformity at the transition disappears. We rationalize this critical behavior using a local mean-field model describing how particle contacts are created via mechanical noise, showing that diffusive motion induced by long-range interactions becomes dominant for slowly-decaying interactions.

academic

Random organization criticality with long-range hydrodynamic interactions

基本信息

论文ID: 2506.14330
标题: Random organization criticality with long-range hydrodynamic interactions
作者: Tristan Jocteur, Cesare Nardini, Eric Bertin, Romain Mari
分类: cond-mat.soft cond-mat.stat-mech
发表时间: 2025年6月17日 (arXiv预印本)
论文链接: https://arxiv.org/abs/2506.14330

摘要

本文研究了软物质非热系统中的可逆-不可逆转变(RIT)，这是一种从吸收态到活跃相的相变。作者通过扩展随机组织模型(ROM)来包含长程流体力学相互作用，发现当相互作用衰减指数小于空间维度时，临界性质连续依赖于衰减指数。随着相互作用范围增加，转变从凹性变为凸性(序参量指数β > 1)，涨落从发散变为消失，转变处的超均匀性消失。作者通过局部平均场模型合理化了这种临界行为。

研究背景与动机

问题背景

软物质系统的动力学停滞: 软物质如泡沫、颗粒材料和乳液在外力驱动下会表现出动力学停滞，形成"吸收态"
周期剪切悬浮液的RIT: 在小振幅振荡剪切下，系统会达到极限环，频闪观察下动力学被冻结；超过临界驱动振幅后，频闪动力学变为扩散性
ROM模型的局限性: 现有的随机组织模型忽略了悬浮粒子间的流体力学相互作用

研究动机

理论挑战: 实验和数值模拟显示RIT为凸性转变(β > 1)，与ROM预测的保守定向渗流(CDP)普适类不符
流体力学效应: 粒子接触产生的长程流体力学相互作用会影响所有其他粒子，创造新的活性产生机制
普适性问题: 需要理解不同边界条件(体相vs受限系统)对RIT性质的影响

核心贡献

扩展ROM模型: 提出了包含长程流体力学相互作用的中介ROM(mediated ROM)，相互作用强度按距离的幂律衰减
临界指数的连续变化: 发现当衰减指数α < d时，临界指数连续依赖于α值
凸性转变的发现: 证明了长程相互作用导致转变从凹性(β < 1)变为凸性(β > 1)
超均匀性的消失: 揭示了长程相互作用会破坏CDP类转变的标志性超均匀性
理论框架: 建立了基于Hébraud-Lequeux类模型的平均场理论，解释了扩散诱导的活性产生机制

方法详解

任务定义

研究在长程流体力学相互作用影响下，周期剪切悬浮液中可逆-不可逆转变的临界行为，特别是临界指数如何依赖于相互作用的衰减指数α。

模型架构

原始ROM模型

系统设置: d维周期边界条件的立方体，包含N个直径为D的球形粒子
活性定义: 当粒子间距离r_ < D时，粒子i和j被定义为活跃的
动力学规则: 活跃粒子在每个时间步获得随机踢击δ_{a,i}，典型大小为Δ_a = D

中介ROM模型

核心创新在于引入被动粒子的长程诱导运动：

空间粗粒化: 将系统划分为大小为ℓ的立方网格，定义活性场： $A_b^t = \frac{1}{\ell^d}\sum_{i \in b} n_i$
长程相互作用核: $G(r) = \frac{c}{(1 + r^2)^α}$ 其中α控制相互作用的衰减速度
被动粒子动力学: 被动粒子i受到随机踢击δ_{p,i}，其振幅为： $\Delta_{p,i} = \sqrt{\sum_{b'} G(r_{b'b})A_{b'}}$
边界条件依赖:
- 体相系统: α = d-1
- 受限系统: α = d

技术创新点

长程噪声的空间化处理: 不同于之前的均匀化处理，本文考虑了空间分布的长程相互作用
系统性的α依赖研究: 首次系统性地研究了临界性质对相互作用范围的依赖
多尺度分析: 结合了粒子层面的动力学和连续场理论的分析

实验设置

数值模拟参数

二维系统: 系统尺寸最大L = 8192，粒子数最多N ≈ 1×10^7
三维系统: 系统尺寸最大L = 768，粒子数最多N ≈ 4×10^7
相互作用范围: α从0.5到3.5变化
初始条件: 粒子位置随机均匀分布

评价指标

序参量指数β: 平均活性⟨A⟩ ∼ δφ^β，其中δφ = φ - φ_c
涨落指数γ': 活性涨落⟨δA²⟩ ∼ δφ^{-γ'}
弛豫指数δ: 临界点处活性衰减A(t) ∼ t^{-δ}
结构因子: S(q) ∼ q^ζ (q → 0)，表征超均匀性

对比方法

标准ROM: α → ∞的短程极限
CDP普适类: 理论预期的临界指数
MROM∞模型: 活跃粒子进行无限大跳跃的平均场版本
Lévy-SHL模型: 基于异常扩散的理论模型

实验结果

主要结果

临界指数的演化

维度	α范围	β值	γ'值	转变性质
d=2	α=3	0.62(2)	0.39(3)	凹性(CDP)
d=2	α=1.5	1.05(2)	-0.10(2)	凸性
d=2	α=0.5	1.82(15)	-1.28(15)	强凸性
d=3	α=3.5	0.84(2)	0.16(6)	凹性(CDP)
d=3	α=2	1.30(3)	-0.44(3)	凸性