微米尺度胶体体系凝胶化中相互作用非均质性（Interaction Heterogeneity, IH）作用的建模研究

《Soft Matter》：Modelling the role of interaction heterogeneity in the gelation of micron-scale colloidal systems

　　

功能化超分子基团的微米尺度胶体是自组装的多功能平台。其中，单链DNA（single-stranded DNA）功能化已被证明是一种高度可调的诱导短程吸引粒子间相互作用的方法，可驱动自组装形成多种晶体结构。最近研究发现，胶体群体中表面功能化程度的差异会导致“相互作用非均质性”（interaction heterogeneity, IH），即一对胶体之间的结合强度随其表面DNA链密度而变化。研究人员在前期工作中已表明IH强烈影响胶体晶体成核与生长，但其对更远离平衡条件下的凝胶化的影响仍鲜有探索。本研究采用分子模拟系统研究了IH在热淬冷驱动的胶体凝胶化中的作用。研究人员考虑了四种IH分布：单分散（无IH）、高斯分布、双分散分布和均匀分布，并分析了它们对凝胶结构及凝胶化动力学的影响。结果表明，虽然IH对宏观凝胶结构影响较小，但通过配位数（coordination number, CN）分布可发现IH深刻改变了凝胶的局部环境。对CN矩的主成分分析（principal component analysis, PCA）突出了IH引起的独特结构趋势，强调了局部凝胶结构对IH的敏感性。研究人员还发现，IH导致强结合体与弱结合体的sequential aggregation（顺序聚集），其中强结合体首先形成凝胶“骨架”，随后弱结合体在其上修饰。这些发现凸显了IH可作为调节凝胶微结构而不显著扰动宏观组织的关键参数。

论文解读：微米尺度胶体体系凝胶化中相互作用非均质性作用的建模研究

一、研究背景与意义

微米尺度、短程吸引胶体是一类重要的自组装体系，广泛存在于食品、化妆品及蛋白质治疗药物中。为克服天然胶体批次间差异大、相互作用定义模糊的局限，研究人员开发了以DNA功能化胶体为代表的工程模型体系。该体系通过DNA互补链杂交产生可逆的、高度特异的粒子间吸引作用，类似于伦纳德-琼斯或莫斯（Morse）势。

此前研究表明，由于DNA接枝过程的不同，胶体群体中会存在DNA表面覆盖度的差异，导致“相互作用非均质性”（Interaction Heterogeneity, IH）——即不同胶体对之间的结合强度不同。研究人员此前的成果显示，IH能调节晶体成核与生长，通过按结合强度分层排序的粒子群体增强结晶鲁棒性并抑制凝胶化。然而，IH在远离平衡态的胶体凝胶化过程中的作用尚不清楚。凝胶化通常在高度非平衡条件下发生，涉及动力学逮捕，传统上难以精确控制。本研究旨在探究IH作为新的调控维度，如何影响热淬冷驱动的胶体凝胶化，论文发表于《Soft Matter》。

二、关键技术方法

研究人员采用Langevin动力学（LAMMPS软件）模拟4000个（部分16000个）粒子在NVT系综下的凝胶化过程。粒子间势采用Rogers-Crocker粗粒化模型描述DNA功能化胶体的吸引与排斥贡献。IH通过为每个粒子i分配结合调节因子b_i实现，两两作用势为U_ij(r)=b_ib_jU_ref(r)，其中U_ref(r)为参考势，分布归一化使< />²>=1。研究考察了四种IH分布：单分散（无IH）、高斯（σ_IH=0.15）、双分散（b_s/w=1.50，组分比可变）和均匀（半宽Δb=0.1–0.5）。通过线性热淬冷协议（dU/dt恒定）驱动凝胶化，引入尺寸多分散性（σ=0.05）抑制结晶。结构分析包括配位数（CN）分布、径向分布函数g(r)、静态结构因子S(k)、盒计数分形维数d_f、拓扑簇分类（TCC）及CN矩的主成分分析（PCA）。

三、研究结果

3.1 无IH时淬冷速率对凝胶结构的影响

以单分散体系为基线，固定体积分数?=0.20，改变淬冷速率dU/dt。快淬冷产生开放枝晶网络，慢淬冷形成更紧凑簇。CN分布从快淬的峰位CN=6偏移至慢淬的CN=7–8，在dU/dt≈?0.20 k_BT s^{?1附近发生转变。通过扩散时间τ_D=L2/D与凝胶化时间τ_g=ΔU/(dU/dt)平衡得到分选长度L≈(D·ΔU/|dU/dt|)1/2，取ΔU～2 k_BT。当L～3倍平均自由程L_MFP（约5粒子直径）时定义为“慢淬冷”阈值（dU/dt～?0.02 k_BT s?1），与CN峰移一致。}

3.2 IH对宏观凝胶结构的影响

在dU/dt=?0.020 k_BT s^?1下比较不同IH分布的最终凝胶。单色渲染下宏观形貌无明显差异；按b值着色可见强结合体（b>1）局部富集区。g(r)与S(k)对IH不敏感，无IH体系在中等k（≈2–4）略高，而双分散/均匀IH在低k(<1)略高，反映大尺度聚集微弱增强。分形维数d_f随淬冷减慢从~1.8（扩散限制簇聚结，DLCA）升至~2.3（反应限制簇聚结，RLCA）。中间淬冷速率下IH体系d_f稍高，但16000粒子大体系结果提示可能为有限尺寸噪声，宏观结构受IH影响甚微。

3.3 IH对局部凝胶结构的影响

CN分布对IH高度敏感。阶跃淬冷（无限快）下IH无影响；慢淬冷下，无IH为单峰，高斯IH使CN=4–8间平坦化，双分散IH（p(b_s)=0.50）产生明显双峰，均匀IH产生较宽双峰。平均CN随淬冷速率变化受IH影响小。拓扑簇分类（TCC）显示IH将中间尺寸（~5–8）簇推向更大更密簇（~10–13），均匀IH促一般致密化，高斯IH特异地稳定尺寸13的高配位局域堆砌。对CN分布前10阶矩做PCA，前两主成分（PC1, PC2）解释>99%方差。无IH数据沿PC1+PC2方向随淬冷速率移动；固定淬冷速率下，IH使数据沿垂直方向偏移，且不同IH类型与程度对应原无IH体系不可及的(PC1,PC2)空间位置，证实IH可定制新局部结构。体积分数?=0.20与0.30的趋势一致，略有偏移。

3.4 IH对凝胶化动力学的影响

以平均CN随参考势阱深U的演化表征动力学。IH使凝胶化过程变宽：平均CN在较低|U|开始上升，在较高|U|达终值。双分散IH清楚展示顺序聚集：强结合体在U～?3.0 k_BT开始成网（平均CN～9），弱结合体在U～?5.0 k_BT跟进装饰（平均CN～4.5），总平均CN呈两阶段平台（~6.5）。构型快照显示无IH体系在U=?3.5到?6.0 k_BT完成聚集；IH体系强结合体先成骨架，弱结合体后填充，延长有效凝胶化窗口。

四、讨论与结论翻译

凝胶化常是高度非平衡、难控的过程，传统控制参数为淬冷速率，可改大尺度形貌但粒子配位分布相对不变。

本研究将相互作用非均质性（IH）确立为热淬冷之外调控胶体凝胶化的新模态。研究人员系统研究了不同IH分布与淬冷速率的组合影响。

在体系整体层面，IH对凝胶结构影响较弱。g(r)与S(k)仅示微差，中间淬冷下分形维数差异在大体系中减弱。

反之，IH深刻影响局部凝胶结构。其效应源于将聚集过程在时间上展宽，类似IH介导的结晶行为。慢淬冷下特定IH类型显著改变CN分布，甚至产生单分散体系无法获得的双峰分布。CN矩的PCA可有效区分淬冷速率与IH的结构指纹，TCC显示IH促进连通拓扑基元形成，且IH分布类型决定网络细观特征。

动力学上，双分散IH将凝胶化分两步：强结合体先聚集成骨架，弱结合体后修饰。强/弱结合体配位环境迥异，为定制凝胶性能提供路径。

潜在应用包括各向同性非晶光子晶体（APCs）：利用IH顺序控制，强结合体先成短程有序畴，弱结合体填隙形成宏观无序凝胶，有望实现各向同性光子带隙。

综上，研究人员证明IH可调制凝胶化中粒子结合序列与动力学，在不显著改变宏观组织下控制局部结构特征，是设计胶体凝胶结构功能（尤指各向同性光子材料）的关键参数。未来需结合3D共聚焦显微等粒子分辨实验及IH定量度量，验证此机理并系统指导光学响应导向的凝胶化协议设计。