[论文解读] Effective interactions in active Brownian particles

Clare R. Rees-Zimmerman, C. Miguel Barriuso Gutierrez+2 more|arXiv (Cornell University)|Jan 27, 2026

Micro and Nano Robotics被引用 0

一句话总结

本文通过将不同方案的径向分布函数匹配，开发了一种反演方法来推导二维有活性布朗粒子（ABP）的有效成对势，并利用这些势来推断有效热力学量。

ABSTRACT

We report an approach to obtain effective pair potentials which describe the structure of two-dimensional systems of active Brownian particles. The pair potential is found by an inverse method, which matches the radial distribution function found from two different schemes. The inverse method, previously demonstrated via simulated equilibrium configurations of passive particles, has now been applied to a suspension of active particles. Interestingly, although active particles are inherently not in equilibrium, we still obtain effective interaction potentials which accurately describe the structure of the active system. Treating these effective potentials as if they were those of equilibrium systems, furthermore allows us to measure effective chemical potentials and pressures. Both the passive interactions and active motion of the active Brownian particles contribute to their effective interaction potentials.

研究动机与目标

• 通过基于有效平衡描述来促使对有活性物质结构的理解。
• 开发并应用反演方法从ABP模拟中提取 beta u_eff(r)。
• 证明一个类似于平衡的势能可以在 MC 模拟中再现 ABP 的 g(r)。
• 研究密度和活性对新 emergent 有效相互作用的影响。

提出的方法

• 使用两种方案计算 g(r)：距离直方图（DH）和测试粒子插入（TPI）。
• 以初始猜测 u_eff,0(r) = -kB T ln[g_DH(r)] 迭代更新 u_eff(r)。
• 应用 Schommers 校正：u_eff,j+1(r) = u_eff,j(r) - kB T ln[g_DH(r)/g_TPI,j(r)]。
• 使用 chi^2 = sum_i (g_DH(r_i) - g_TPI,j(r_i))^2 收敛。
• 在不同被动势（LJ、WCA、肩部势）和活性（Pe）下对 ABP 进行模拟以获得 g(r) 供反演使用。
• 通过用回收得到的 beta u_eff(r) 进行被动 MC 模拟并将 g(r) 与 ABP 结果对比来验证。

Fig. 1: Zoomed snapshots of the MD (top row) and MC (bottom row) simulations for the studied potentials for selected densities and Péclet numbers. From left to right: From the shoulder potential study at varying $\rho$ , we present $\rho=0.4$ (at $\rm{Pe}=120$ , see Fig. 5 ). Then, for the study var

实验结果

研究问题

• RQ1 是否可找到一个有效对势 u_eff(r) 通过平衡 MC 描述再现 ABP 的结构？
• RQ2 活性（Pe）和被动势类型如何影响 u_eff(r) 的形式？
• RQ3 平均力势 w(r) 是否为良好的初始猜测，以及它与 u_eff(r 的关系？
• RQ4 可以从推导的 u_eff(r) 推断出哪些热力学量（压力、化学势）以实现对 ABP 的热力学描述？
• RQ5 密度如何影响 ABP 的 emergent u_eff(r) 的密度依赖性并体现出非成对相互作用？

主要发现

• 通过比较 ABP 的 DH 和 TPI 方案得到的 g(r)，可获得有效对势 u_eff(r)。
• 反演方法在 LJ、WCA 和肩部被动势中产生的 u_eff(r) 具有定性相似性，活性通过在 u_eff(r) 中引入吸引特征实现。
• 随着 Pe 增大，u_eff(r) 的引力井初期加深，但在较高 Pe 下对于 WCA 势，井深度下降。
• 将有效对势在被动 MC 模拟中重新再现结构，表明它捕捉了主导的结构特征。
• 密度改变 u_eff(r)，并增强了两个局部极大，体现出新出现的密度相关的非成对效应。
• 由推导出的 u_eff(r) 可以推断有效化学势，以及一个概念性的有效压力，为 ABP 的热力学描述提供路径。

(a)