Quasi-two-dimensional dispersions of Brownian particles with competitive interactions: Dynamical clustering, non-Gaussianity and hydrodynamic correlations

技术摘要：具有竞争相互作用的布朗粒子准二维弥散体系

问题陈述
本文研究了受竞争性短程吸引（SA）与长程排斥（LR）相互作用（即 SALR 系统）影响的准二维（Q2D）布朗粒子弥散体系的动力学行为。虽然这些体系的结构特性——特别是有限尺寸簇（clusters）和微相（microphases）的形成——在三维（3D）和 Q2D 背景下已有充分文献记载，但对其动力学的系统理解仍不完整。关键的开放性问题包括：

1. 流体动力学相互作用（HIs）如何影响受限 Q2D 几何结构中的簇稳定性与集体输运？
2. 聚簇现象与 Q2D 系统特有的反常集体扩散（即当 $q \to 0$ 时，$D(q)$ 和 $H(q)$ 发散）之间的相互作用如何体现？
3. 控制聚簇相中非高斯动力学和动力学异质性的机制是什么？
4. 与硬球参考体系相比，短时惯性效应和流体动力学延迟如何影响聚簇相中集体相关性的建立？

研究方法
作者采用两种不同的模拟技术进行对比研究，以分离流体动力学的影响：

1. 朗之万动力学 (LD)： 用于模拟不存在流体动力学相互作用的系统。这作为过阻尼动力学的参考。
2. 多粒子碰撞动力学 (MPC)： 用于显式解析涨落流体动力学，捕捉惯性和过阻尼机制。溶剂由在碰撞单元内进行随机旋转动力学（SRD）的点粒子组成，确保了质量、动量和能量的守恒。

该系统由 $N_p$ 个相同的球形粒子组成，这些粒子被限制在三维牛顿溶剂内的平面单层中。相互作用势结合了陡峭的广义 Lennard-Jones 吸引力和屏蔽静电（Yukawa）排斥力。研究通过系统地改变吸引强度（$\epsilon^*$）和面积分率（$\phi_{2D}$）来绘制相行为图谱。

关键分析工具包括：

• 键相关函数 ($B(t)$) 和 簇寿命 ($\tau_b$) 分析，用于表征动力学停滞。
• 局部六角序参数 ($q_6$) 及其时间相关函数 ($g_6(t)$)，用于评估结构弛豫。
• 非高斯参数 ($\alpha_2$) 和自 Van Hove 函数，用于量化动力学异质性。
• 流体动力学函数： 包括静态 $H(q)$ 和随时间变化的 $H(q,t)$，并将其分解为自部分（self）和不同部分（distinct），以探测集体扩散和动量相关性。
• 纵向电流-电流相关函数 ($J(q,t)$)，用于研究声模式传播。

主要贡献与结果

• 相行为与热化： 研究确认了四种广义相：弥散流体（DF）、平衡簇（EC）、随机渗透（RP）和簇渗透（CP）。热化过程得到了验证，表明除了极低温度下出现的缓慢能量漂移外，系统能够达到独立于初始条件的稳态。

• 自扩散与次扩散： 随着吸引强度的增加，自扩散受到显著抑制。

  • 在流体相（DF, EC, RP）中，粒子表现出明显的短时和长时布朗扩散，这种减速是由瞬态聚簇而非笼效应（caging）引起的。
  • 在簇渗透（CP）相中，出现了两步弛豫，对应于两个不同的次扩散机制，这与凝胶形成相一致。

• 簇寿命与转变判据：

  • 研究发现，簇寿命 ($\tau_b$) 受吸引强度的影响比粒子浓度更为显著。
  • 确定了两种用于区分非聚簇相与聚簇相的动力学判据：
    1. 平均簇寿命阈值：$\tau_b / \tau_D \approx 10$ 区分了 DF 和 EC 相。
    2. 局部六角序的弛豫时间 ($\tau_6$)：交叉点 $\tau_6 \approx \tau_D$ 标志着 DF–EC 的转变。
  • 在 EC 相中，簇寿命符合 Kramers 势垒越过理论；而在 DF 机制中，则遵循 Smoluchowski 解离理论。

• 非高斯动力学： 聚簇 Q2D-SALR 系统表现出显著的非高斯动力学。平衡簇相中的自 Van Hove 函数呈现近似指数形式，这与“扩散性扩散”（diffusing-diffusivity）机制一致。这表明存在显著的动力学异质性，即簇内的粒子运动受限，而弥散相中的粒子保持移动性。

• 流体动力学相互作用 (HIs) 与集体扩散：

  • MPC 模拟显示，Q2D 系统中特征性的异常增强的大尺度集体扩散在 Q2D-SALR 弥散体系中得到了定性保留。
  • 然而，与 3D 对应体系相比，聚簇现象抑制了流体动力学函数中的中程有序峰。
  • 至关重要的是，对随时间变化的流体动力学相关性 ($H(q,t)$) 和声模式的分析表明，聚簇导致了流体动力学相互作用更早的出现，其时间尺度早于 Q2D 硬球参考系统。这意味着在聚簇 SALR 系统中，HIs 在惯性时间尺度上便已发挥作用。

意义与主张
本文旨在提供全面的动力学分析，以弥合受限 SALR 体系中结构微相分离与输运性质之间的鸿沟。通过结合 LD 和 MPC，作者证明了：

1. 吸引强度是动力学聚簇的主导控制参数，其影响力超过了粒子浓度的影响。
2. 流体动力学相互作用在 Q2D-SALR 弥散体系的集体动力学中起着至关重要的作用，既保留了反常的长程扩散特性，又修改了流体动力学函数的谱特征。
3. 聚簇加速了流体动力学的相关性，导致 HIs 在惯性时间尺度上显现，这一发现挑战了此类体系理论描述中关于瞬时 HIs 的假设。
4. 本研究提供了动力学判据（基于簇寿命和六角序弛豫），用以区分弥散相与聚簇相，从而补充了现有的结构判据。

该工作强调，若要全面理解受限胶体弥散体系中动力学聚簇和集体输运的微观机制，必须同时考虑竞争相互作用与流体动力学耦合。