Theory of Nonequilibrium Crystallization and the Phase Diagram of Active Brownian Spheres

技术摘要：非平衡结晶理论与活性布朗球体的相图

问题陈述
硬球在平衡态下的结晶是熵驱动相变的典型范例。然而，最近的观测表明，活性（广泛存在于从运动细菌到自驱动胶体等系统中）显著改变了这种有序-无序转变。具体而言，有限的活性导致单分散硬球的固相密度迅速增加，趋向于最大面心立方（fcc）堆积密度（$\phi_{solid} \approx 0.74$），这与平衡态值（$\phi_{solid} \approx 0.545$）形成鲜明对比。尽管有这些观测结果，开发活性结晶的理论描述仍然难以实现，因为支撑被动冻结理解的传统热力学论据并不适用于活性非平衡系统。目前迫切需要一种理论来阐明活性结晶背后的驱动力，并预测固态-流体相图，而不依赖于平衡热力学。

方法论
作者开发了一种用于描述三维活性布朗粒子（ABPs）结晶的多尺度统计力学理论。该方法包含以下关键步骤：

1. 场论构建： 系统使用两个序参量场进行建模：一个守恒的数量密度场 $\rho(z)$（或体积分数 $\phi(z)$）和一个非守恒的、无量纲的“结晶度”场 $\psi_N(z)$（或结晶度密度 $\psi_V(z) = \psi_N \rho$）。这允许一种描述方式，即流体相和晶体相都具有均匀的密度场，而结晶度仅由序参量捕捉。
2. 推导共存准则： 鉴于不存在非平衡变分原理，作者利用基于稳态条件（$\dot{\rho}=0, \dot{\psi}_V=0$）的动力学框架开发了共存准则。
   • 他们定义了一个“结晶度伪势”（$u_\psi$）和一个“化学伪势”（$u_\rho$）。
   • 不同于平衡系统中吉布斯-杜恒关系允许严格的麦克斯韦构造，作者发现对于 ABPs，非平衡吉布斯-杜恒关系无法被精确满足。
   • 因此，他们提出了一个近似的广义麦克斯韦构造。这涉及求解一个变换张量以建立伪势的均匀性，并调用界面非平衡吉布斯-杜恒关系来近似最终的共存准则。
3. 状态方程 (EOS)： 在计算机模拟（布朗动力学）的指导下，作者开发了物理引导的经验状态方程，包括保守相互作用压力（$p^C_{bulk}$）、活性压力贡献（$U$）以及平衡结晶度路径（$\psi^*_V(\rho)$）。这些状态方程捕捉了系统对堆积率（$\phi$）和活性路程与粒子直径之比（$\ell_0/D$）的依赖关系。

主要贡献与结果

• 有序-无序转变的偏移： 推导出的状态方程表明，增加活性会将有序-无序转变向更高的堆积率方向移动。有序-无序体积分数（$\phi_{ODT}$）随活性单调增加，并在较低活性（$\ell_0/D \approx 5$）时饱和于随机紧密堆积极限（$\phi_{RCP} \approx 0.645$）。
• 重建相图： 通过将近似共存准则应用于开发的状态方程，作者构建了活性布朗球体的完整相图。该理论定量地重现了：
  • 固-流体共存曲线，正确预测了随着活性增加，固相密度向紧密堆积快速靠近的过程。
  • 液-气共存曲线（运动诱导相分离，即 MIPS）。
  • 固-液-气三相点以及随后在高活性下向固-气共存的转变。
• 对模拟的验证： 理论预测在各种中高活性范围内与模拟数据表现出强一致性。该理论成功捕捉了 MIPS 和活性结晶的 binodals（分界线），这一成就此前在不诉诸平衡热力学的情况下是无法实现的。
• 平衡准则的失效： 本文证明，天真地应用平衡共存准则（通过标准吉布斯-杜恒关系实现等化学势和等压力）会导致定性错误。平衡准则会极大地高估流体密度，并错误地预测在三相点之上存在固-液共存，而实际上系统表现为固-气共存。

意义
作者声称，这项工作提供了第一个能够定量重现活性布朗球体完整相图且不依赖于平衡热力学论据的活性结晶统计力学描述。通过识别特定的非平衡共存准则并开发相应状态方程，该理论阐明了导致有序-无序转变偏移的驱动力。这项工作为通向通用的非平衡结晶理论提供了一条具体的路径，强调了平衡吉布斯-杜恒关系对于活性系统在形式上是不正确的，且需要不同的非平衡准则才能准确预测相行为。