[论文解读] A fully solvable equilibrium self-assembly process: fine tuning the clusters size and the connectivity in patchy particle systems

Emanuela Bianchi, Emilia La Nave+2 more|ArXiv.org|Aug 14, 2007

Pickering emulsions and particle stabilization参考文献 22被引用 4

一句话总结

本文提出了一套无参数的理论框架，可完全预测具有可调簇尺寸和连通性的钉状粒子系统中的平衡自组装行为。通过结合Wertheim的缔合流体理论与Flory-Stockmayer的凝胶化理论，该理论在二元混合物（两功能与三功能粒子）中实现了对簇分布、渗透阈值及相行为的精确预测，经大量蒙特卡罗模拟验证，理论与模拟在所有热力学状态下完全一致。

ABSTRACT

Self-assembly is the mechanism that controls the formation of well defined structures from disordered pre-existing parts. Despite the importance of self-assembly as a manufacturing method and the increasingly large number of experimental realizations of complex self-assembled nano aggregates, theoretical predictions are lagging behind. Here we show that for a non-trivial self-assembly phenomenon, originating branched loop-less clusters, it is possible to derive a fully predictive parameter-free theory of equilibrium self-assembly by combining the Wertheim theory for associating liquids with the Flory-Stockmayer approach for chemical gelation.

研究动机与目标

• 开发一种完全可预测、无参数的理论，用于描述具有非球形、方向性相互作用的钉状粒子系统中的平衡自组装行为。
• 理解在强而局域化的吸引力作用下，粒子功能度、温度与密度的相互作用如何导致簇尺寸与连通性的形成。
• 识别并表征三功能粒子作为分支点时所引发的热力学不稳定性与渗透转变。
• 在广泛的热力学状态点上，通过大规模蒙特卡罗模拟验证理论预测。
• 证明键概率 $p_b$ 是簇结构的唯一有效控制参数，从而实现对材料性能的连续调控。

提出的方法

• 系统被建模为5670个两功能与330个三功能硬球粒子的二元混合物，具有特定位点的正方形阱吸引势。
• 通过深度为 $u_0$、范围为 $\delta = 0.119\sigma$ 的正方形阱势定义相互作用，由于空间位阻，每个位点仅能形成单键。
• 应用Wertheim的缔合流体理论推导状态方程，并预测双结线与渗透线。
• 采用Flory-Stockmayer方法计算键概率 $p_b$ 的函数关系下的簇尺寸分布及有限簇的数量。
• 在70多个状态点上进行蒙特卡罗模拟，通过改变温度 $T$ 与密度 $\rho$，利用粒子平动与转动采样平衡构型。
• 键概率 $p_b$ 定义为总势能与基态能量之比，从而实现模拟与理论的直接对比。

实验结果

研究问题

• RQ1 无参数理论框架能否准确预测具有非平凡连通性与分支结构的钉状粒子的平衡自组装行为？
• RQ2 平均功能度 $\langle M \rangle$ 如何影响此类系统中渗透与相分离的出现？
• RQ3 键概率 $p_b$ 在多大程度上可作为簇尺寸与网络拓扑结构的通用控制参数？
• RQ4 在全部热力学条件下，有限簇数量 $N_c$ 与键概率 $p_b$ 之间存在何种关系？
• RQ5 Wertheim理论与Flory-Stockmayer理论在捕捉相行为与簇统计特性方面，与大规模蒙特卡罗模拟的预测结果相比如何？

主要发现

• Wertheim与Flory-Stockmayer理论的预测结果在所有状态点上与蒙特卡罗模拟完全一致，验证了无参数方法的有效性。
• 键概率 $p_b$ 是唯一有效的控制参数：具有相同 $p_b$ 的构型表现出相同的簇尺寸分布，无论 $T$ 与 $\rho$ 如何。
• 系统在 $p_b \approx p_b^c$ 处表现出明确的渗透转变，此时首次出现贯穿整个体系的簇，模拟与理论结果均证实了这一点。
• 有限簇数量 $N_c$ 随 $p_b$ 线性减少，遵循 $N_c = N - p_b \cdot (2N_2 + 3N_3)/2$，模拟与理论结果高度一致。
• 相图显示了由Wertheim理论得到的双结线、比热的极大值以及渗透位点，所有结果均与模拟数据定量吻合。
• 在高密度与低温条件下，系统表现出类似气-液相分离的行为，非平衡态中密度涨落呈双峰分布。