分子自组装模拟：一个表面活性剂胶束体系的微观结构分析

自组装是自然界中最精巧的组织方式之一——表面活性剂分子在溶液中自发地排列成胶束、囊泡、层状相等有序结构，不需要任何外部指令。但在分子层面，这个”自发”过程到底是怎么发生的？哪些驱动力在主导分子的排列选择？分子动力学模拟能提供这些肉眼看不到的信息。

项目背景：胶束形成的动力学过程

团队处理过一个十二烷基硫酸钠（SDS）胶束的形成与结构分析项目。客户在做SDS基纳米乳液的配方优化，需要知道在给定的温度和盐浓度下，SDS形成的胶束是球形的还是棒状的，平均聚集数是多少，核心的疏水区域尺寸有多大——这些信息直接影响乳液的稳定性和药物包载效率。

实验手段（动态光散射、小角X射线散射）可以给出聚集数和尺寸的统计平均值，但无法直接观察到单个胶束的形成过程和内部分子排列。分子动力学模拟的切入点正是这个实验手段的盲区。

模型构建与力场参数

SDS是阴离子表面活性剂，分子由疏水尾链（十二烷基，C12H25）和亲水头基（硫酸基，SO4⁻）组成。团队使用GROMOS54a7力场——这个力场专门针对生物和表面活性剂分子优化过，对碳氢链的扭转能面描述比CHARMM力场更符合实验数据。

模拟体系的构建方式是：将100个SDS分子随机分散在一个10×10×10 nm³的立方盒子中，加Na⁺离子中和SO4⁻的电荷，再用TIP4P水模型填充。初始状态下SDS分子均匀分布，没有预组装的胶束结构——让体系自己”决定”怎么组装。

NaCl浓度设为0.1 M，温度298 K，NPT系综。体系总原子数约78000个。

从无序到有序：胶束自发形成

模拟的前20 ns是最精彩的阶段。随机分散的SDS分子在水中快速重排——疏水尾链为了逃离水环境而相互聚集，亲水头基朝向水相。大约在5-8 ns内，大部分SDS分子就已经参与到了某种聚集体中。但早期的聚集体形态不规则，大小不一，存在大量的”融合-分裂”事件。

到20 ns时，体系基本稳定下来——形成了约12-15个胶束，其中3-4个较大的球形胶束（聚集数30-45），其余为小胶束或预胶束聚集体（聚集数10-20）。这个分布与实验上用SAXS测得的平均聚集数62±8（在CMC附近）有合理的一致性——考虑到模拟体系的浓度（约100 mM）高于CMC（8.2 mM），小胶束的存在是合理的。

胶束内部结构的关键发现

径向分布函数（RDF）的分析揭示了胶束内部的精细结构。以胶束质心为参考点，疏水尾链的碳原子密度在距中心约1.2 nm处达到峰值，头基的硫原子密度在约1.8 nm处形成壳层。整个胶束的疏水核半径约1.5 nm，与实验中SAXS拟合的核半径1.48±0.05 nm吻合良好。

更有趣的是水分子在胶束中的渗透行为。RDF显示有约5%的水分子能够渗透到头基壳层内部，到达疏水核的边缘区域。这种”水的渗透”在热力学上是不利的，但在动力学上是不可避免的——它对胶束的动态稳定性有重要贡献，也是药物分子从胶束中释放的微观通道之一。

50 ns的模拟在128核CPU集群上运行了约3天。客户根据模拟结果调整了配方中的SDS浓度和盐浓度，使胶束尺寸更均匀，乳液稳定性提高了约30%。自组装模拟的价值在于：它不只是给出一组统计数字，而是展示了这些数字背后的微观过程——这种过程层面的理解对配方优化有直接的指导意义。

如需了解更多的科研计算案例与项目经验，可以参考 https://www.keyanxueshu.com/ 上的技术文章。