钙离子调控的两性离子/阳离子表面活性剂-反离子复合物的自组装：蠕虫状胶束生长被抑制的分子动力学洞察

《Molecules》：Calcium-Regulated Self-Assembly of Zwitterionic/Cationic Surfactant–Counterion Complexes: Molecular Dynamics Insights into the Suppression of Wormlike Micelle Growth

　　

蠕虫状胶束(Wormlike Micelles, WLMs)为清洁压裂液的粘弹性提供了结构基础，然而含钙卤水调控其生长的分子机制尚不充分明确。在此研究中，研究人员采用全原子分子动力学(Molecular Dynamics, MD)模拟研究了EAHSB/EHAC/NaPts表面活性剂-反离子体系在钙离子(Ca2+)调控下的自组装过程。在不含Ca2+的情况下，分散的表面活性剂和反离子分子首先形成小的球形聚集体，随后融合成较大的胶束、棒状中间体，最终形成蠕虫状胶束(WLMs)。增加局部表面活性剂/反离子负载进一步促进了一维胶束生长和网络状关联。当引入Ca2+时，组装过程变慢，最终聚集体形态从细长的蠕虫状胶束(WLMs)转变为较小的球形聚集体。径向分布函数(Radial Distribution Functions, RDFs)和配位数(Coordination-number)分析表明，Ca2+重新分布了混合表面活性剂/反离子聚集体内的带电基团相关性。水合(Hydration)分析表明极性水合壳层增大，而疏水尾链水合和旁式缺陷(gauche-defect)概率揭示了增强的尾部暴露和更松散的链堆积。这些结果表明，Ca2+通过静电重分布、界面水合和疏水链无序化的耦合作用抑制蠕虫状胶束(WLMs)的生长，从而有利于高曲率球形聚集体而非延伸的蠕虫状网络。

该研究发表于《Molecules》。研究背景聚焦于蠕虫状胶束(Wormlike Micelles, WLMs)作为清洁水力压裂液中的“活聚合物”，其宏观流变性受微观自组装过程调控。混合离子表面活性剂体系虽能通过静电作用促进WLMs生长，但其生长对电解质条件高度敏感，特别是油田卤水中普遍存在的二价阳离子如钙离子(Ca²⁺)可能通过静电、水合及离子特异性效应改变聚集体形态与流变性。目前，Ca²⁺调控两性离子/阳离子表面活性剂-反离子复合物自组装的分子机制，尤其是其如何抑制长蠕虫状胶束生长并使组装终点转向小尺寸高曲率聚集体的机制仍不明确。为此，研究人员选取EAHSB/EHAC/NaPts体系，利用全原子分子动力学(Molecular Dynamics, MD)模拟探究无钙及含钙条件下自组装路径的差异，旨在阐明Ca²⁺抑制WLMs生长的分子机理，为耐盐粘弹性表面活性剂流体的分子设计提供理论指导。研究人员发现，无Ca²⁺时体系经历小球融合、棒状中间体至WLMs的网络化生长；引入Ca²⁺后组装速率减慢且终态转为球形聚集体，该抑制源于静电重分布、界面水合增强及疏水尾链无序化三者的耦合作用。

作者为开展研究用到的主要关键技术方法包括：采用CHARMM通用力场(CGenFF)参数化两性离子表面活性剂EAHSB、阳离子表面活性剂EHAC及有机反离子Pts^?，搭配TIP3P水模型及CHARMM兼容离子参数；利用Packmol构建不同局部负载及Ca²⁺比例(0:1:1:1、10:1:1:1、20:1:1:1)的周期性模拟盒子；使用GROMACS 2025.3进行全原子分子动力学(MD)模拟，经最速下降法能量最小化及NPT系综(Berendsen弱耦合控温控压)平衡后，以1 fs步长运行210 ns生产轨迹；通过VMD 1.9.3可视化聚集体形态；计算径向分布函数(RDFs)、第一配位层配位数、极性基团与疏水尾链水合数及烷基链旁式缺陷(gauche-defect)概率以分析微观结构变化；结合经典堆积参数(packing parameter)理论解释形态转变。

研究结果如下。

2.1. Concentration-Dependent Self-Assembly of EAHSB/EHAC/NaPts。研究人员固定EAHSB:EHA⁺:Pts^?摩尔比1:1:1并改变总分子负载，发现随局部负载增加，平衡聚集体从紧凑球形或椭球形胶束转变为长圆柱形胶束，高负载下形成延伸交织网络，符合混合离子表面活性剂体系的一维生长与缠结特征，该系列作为局部富集模型用于比较负载对聚集及细长胶束形成的影响。

2.2. Time-Resolved Pathway from Spherical Aggregates to Wormlike Micelles。研究人员选取中间负载的System 3追踪组装路径，发现分散分子2 ns内形成小球，10 ns融合为大球，30 ns出现棒状中间体，90 ns可见首条蠕虫状胶束，150 ns残余小球并入形成单一长直蠕虫状胶束(直径约5 nm)，90–150 ns与150–210 ns窗口最大聚集体尺寸分别为175与180分子，波动小于3%，表明210 ns内达稳定平台；提出融合控制的生长路径及EAHSB、EHA⁺、Pts^?在聚集体内的分子排列示意图，指出带电基团缔合、反离子组织、界面水合与疏水链堆积的平衡主导球向蠕虫状转变。

2.3. Ca²⁺ Regulates Both the Rate and Endpoint of Self-Assembly。研究人员固定表面活性剂与反离子数各180，设置Ca²⁺比例0:1:1:1、10:1:1:1、20:1:1:1，发现无Ca²⁺体系组装为蠕虫状胶束，含钙体系终态主要为球形且尺寸随Ca²⁺负载增加而减小；无Ca²⁺时2 ns已现较大球聚体，含钙体系分别约10 ns与30 ns才出现可比尺寸，表明Ca²⁺调控自组装动力学路径与终态形态，源于界面带电基团相关性重组、极性水合增强及疏水链堆积减弱。

2.4. Charged-Group Correlations and Coordination Numbers。研究人员以20:1:1:1体系为重点分析径向分布函数(RDFs)与配位数，发现磺酸根-磺酸根、季铵-季铵、季铵-磺酸根等多组峰强降低，EAHSB_S-EAHSB_S配位数从0.97降至0.72，EAHSB_N-EAHSB_N从2.20降至1.76，EAHSB_N-EHAC_N从2.20降至1.81，EHAC_N-EHAC_N从1.90降至1.55，EHAC_N-Pts_S从2.46降至2.16，Pts_S-Pts_S从0.41降至0.05，EAHSB_S-EHAC_N从1.69降至1.40，EAHSB_S-Pts_S从1.08降至0.92，而EAHSB_S-EAHSB_N从2.95升至3.62，EAHSB_N-Pts_S近不变；表明Ca²⁺非简单非特异性屏蔽，而是重分布局部带电基团相关性，削弱跨物种分子间缔合，增强EAHSB中心关联性，削弱持续一维生长所需的分子间关联。

2.5. Hydration of Polar Groups and the Shift in Packing Tendency。研究人员分析极性基团与Pts^?水合数，发现Ca²⁺使表面活性亲水基团与Pts^?水合数升高，极性界面保留更多水分子，水合壳层增大使表观亲水头基面积(a)增大；依据堆积参数p = v/(a·l)，a增大会减小p并倾向高界面曲率，使低曲率圆柱或蠕虫状胶束不利，小而高曲率球形聚集体更稳定，从堆积参数角度解释形态偏移。

2.6. Hydrophobic-Tail Exposure and Conformational Disorder。研究人员分析季铵氮与尾碳径向分布函数(RDFs)发现含钙体系峰向短距移动，更多尾段靠近极性界面；尾碳水合数呈W形分布且含钙体系整体更高，疏水核屏蔽减弱；旁式缺陷(gauche-defect)概率分析显示含钙体系概率更高，表明尾链弯曲增加、有序性降低、堆积松散；长蠕虫状胶束需紧密有序疏水链沿轴堆积，Ca²⁺关联的尾暴露、水合增强与构象无序削弱核连续性，难以维持低曲率延伸形态，构成转向球形聚集体的结构原因。

2.7. Proposed Mechanism for Ca²⁺-Regulated Suppression of Wormlike Micelle Growth。研究人员综合轨迹提出耦合机制：无Ca²⁺时EAHSB、EHAC与Pts^?静电缔合与疏水尾链堆积驱动小球融合为蠕虫状胶束并形成网络；含钙体系中，径向分布函数(RDFs)与配位数显示跨物种带电基团相关性削弱，界面混合带电层重组；水合数显示极性基团与Pts^?水合增强、表观头基面积增大；尾水合与旁式缺陷(gauche-defect)显示疏水核更松散无序；三者耦合降低聚集体维持低曲率延伸形态能力，使组装变慢且终态由蠕虫状胶束转为较小球形聚集体，体现Ca²⁺特异性抑制一维聚集生长。

总结讨论部分，研究人员指出Ca²⁺对蠕虫状胶束(WLMs)生长的抑制源于耦合分子效应而非简单静电屏蔽：径向分布函数(RDFs)与配位数分析表明Ca²⁺重分布EAHSB、EHAC与Pts^?间局部带电基团相关性，削弱支撑混合表面活性剂/反离子组织的多种分子间接触；水合分析进一步表明Ca²⁺扩大极性头基与反离子水合壳，增加有效亲水面积并使堆积倾向转向高曲率聚集体；同时疏水尾链水合增加与旁式缺陷(gauche-defect)概率升高揭示链堆积更松、胶束核致密性下降。总体而言，含钙环境通过带电基团重分布、界面水合增强与疏水尾链无序化的耦合作用抑制蠕虫状胶束(WLMs)生长，削弱持续一维胶束生长并利于小尺寸高曲率球形聚集体。这些分子层面的洞察为设计盐响应或耐盐粘弹性表面活性剂流体提供指导。未来工作应结合更长时间更大尺度模拟与稳健定量形态描述符，并通过流变学、冷冻透射电镜(cryo-TEM)、小角中子/ X射线散射(SANS/SAXS)或核磁共振(NMR)实验验证将分子形态与宏观粘弹性关联。

研究结论部分翻译如下：在此研究中，研究人员采用全原子分子动力学(Molecular Dynamics, MD)模拟研究了EAHSB/EHAC/NaPts表面活性剂-反离子体系的钙调控自组装。在无钙体系中，分散的表面活性剂与反离子分子首先组装为小球形聚集体，随后融合为较大胶束、棒状中间体，最终形成蠕虫状胶束(Wormlike Micelles, WLMs)；增加局部表面活性剂/反离子负载进一步促进一维胶束生长与网络状关联。相比之下，引入钙离子(Ca²⁺)减慢了组装过程，并将最终聚集体形态从细长的蠕虫状胶束(WLMs)转变为较小的球形聚集体；这些结果表明Ca²⁺改变了表面活性剂/反离子自组装的动力学路径与形态终点。Ca²⁺对蠕虫状胶束(WLMs)生长的抑制源于耦合分子效应而非简单静电屏蔽；径向分布函数(Radial Distribution Functions, RDFs)与配位数分析表明Ca²⁺重分布EAHSB、EHAC与Pts^?间的局部带电基团相关性，削弱多种支持混合表面活性剂/反离子组织的分子间接触；水合分析进一步表明Ca²⁺扩大极性头基与反离子的水合壳层，增加有效亲水面积并使堆积倾向转向高曲率聚集体；同时疏水尾链水合增加与旁式缺陷(gauche-defect)概率升高揭示更松散的链堆积与较不致密的胶束核。总体而言，含钙环境通过带电基团重分布、界面水合增强与疏水尾链无序化的耦合作用抑制蠕虫状胶束(WLMs)生长，削弱持续一维胶束生长并利于较小高曲率球形聚集体；这些分子洞察为设计盐响应或耐盐粘弹性表面活性剂流体提供指导；未来工作应结合更长时间更大尺度模拟与稳健定量形态描述符，如回转半径、球形度、形状各向异性、形状因子与主转动惯量，尤其对跨越周期边界的聚集体，并通过流变学、冷冻透射电镜(cryo-TEM)、小角中子/ X射线散射(SANS/SAXS)或核磁共振(NMR)测量实验验证，进一步帮助连接分子形态与宏观粘弹性行为。

要不要我帮你把这篇解读里的关键分子机制（静电重分布、界面水合、疏水链无序化）整理成一张简洁的逻辑关系图来说明？