2.2 气/液界面微观行为的表征能力

2.2.1 界面结构与组分分布：数密度分布分析

图4为不同复配体系下z轴方向的数密度分布。由图4可知，沿z轴的数密度分布可清晰界定气相、界面区与液相的位置及界面厚度。通过该分析能直观地揭示氟碳链、氢碳链、水分子及各复配组分在界面区域的富集程度、分布轮廓及其与纯体系的差异，不同原子/基团（如极性头基与疏水尾链）在z轴上的位置偏移间接反映分子的取向偏好。

图4   不同复配体系下z轴方向的数密度分布

（a） PFHxA/SDS/Li+ （b） PFHxA/SDS/Na+ （c） PFHxA/SDS/K+

2.2.2 分子取向与有序度：动力学分析

为定量表征界面分子（特别是氟碳链与氢碳链）的排列有序性与取向，键序参数（P₂）和分子轴序参数（S）被广泛应用。这些参数能精确量化链的排列有序度、平均倾斜角度（直立、平躺或倾斜）以及头基（如羧基、磺酸基）插入水相的深度。通过直接对比发现，氟碳链通常具有比氢碳链更高的刚性与排列紧密性。

2.2.3 分子间相互作用空间关联：径向分布函数（RDF）分析

RDF可定量揭示界面区域分子间的特定相互作用。图5为复配体系的径向分布函数。其中，g（r）表示距离参考粒子r处找到目标粒子的相对概率密度。通过分析特定原子对（如F…F、F…H（水）、C（氟碳）…C（氢碳）等）的RDF，可探究氟碳链间潜在的氟⁃氟作用、氟碳/氢碳链与水分子间的特征水合结构（水合壳层）以及不同组分间的聚集倾向，为阐明复配体系中链间是倾向于分相聚集还是混合分布、是否存在特定分子间协同作用提供关键证据。

图5   复配体系的径向分布函数

2.2.4 微观构象与聚集状态：可视化与辅助分析

对界面微观结构与形貌的可视化分析（包括轨迹快照、二维密度俯视图）可提供分子构象、排列及聚集状态的直观图像，辅以界面涨落分析及团簇分析，可有效探测界面在横向（x⁃y平面）的均质性、微相分离现象（如氟碳/氢碳富集区）、微小聚集体（如半胶束）的形成以及界面粗糙度。图6为x⁃y平面内表面活性剂分子的二维数密度分布图。

图6   x⁃y平面内表面活性剂分子的二维数密度分布图

（a） PFHxA/SDS/Li+ （b） PFHxA/SDS/Na+ （c） PFHxA/SDS/K+

2.2.5 分子动态行为：动力学分析

通过计算均方位移（MSD）、转动相关函数、氢键动力学等指标，可动态评估界面分子的流动性（平动与转动扩散）、链的柔性（摆动、扭曲幅度与频率）、吸附稳定性（翻转、脱附难易度）以及界面水分子的结构动态特征。图7为不同复配体系中表面活性剂分子的关键动态性质随时间的演变曲线。

图7   不同复配体系中表面活性剂分子的关键动态性质随时间的演变曲线

2.2.6 宏观性质关联：界面张力计算

基于Kirkwood⁃Buff公式或压力张量计算模拟体系的界面张力，将微观模拟结果与关键宏观性质直接关联。通过对比纯组分模拟值、理论加和值与复配体系的模拟值，可定量验证复配体系降低界面张力的协同效应及其程度。

2.2.7 吸附热力学与协同机制：自由能计算

利用伞形采样（Umbrella Sampling）等自由能计算方法，可定量评估不同组分（氟碳/氢碳链分子或基团）在界面上的吸附自由能，揭示其自发吸附界面的热力学驱动力强弱。该方法能深入探究复配体系中一种组分的吸附如何影响另一种组分的吸附自由能，从而阐明协同作用（或拮抗作用）产生的热力学微观机制。

2.3 多组分体系协同机制的解析

MD模拟在原子尺度解析短链氟碳复配体系（如阴/阳离子型、氟碳/碳氢混合型）的多组分协同机制方面具有独特优势，能定量揭示界面分子排布、相互作用及微观⁃宏观关联。图8为不同n（PFB⁃MC）/n（1⁃OA）复配体系的最终MD模拟模型。

图8   不同n（PFB⁃MC）/n（1⁃OA）复配体系的最终MD模拟模型

（a） 1∶0 （b） 2∶1 （c） 1∶1 （d） 1∶2 （e） 0∶1

2.3.1 全原子模型揭示阴/阳离子静电互补驱动机制

X.L.QUAN等利用全原子模型（CHARMM力场）结合多尺度分析（密度分布、RDF、SASA），定量揭示了阴/阳离子表面活性剂在气⁃液界面的协同机制。结果表明，静电互补（-SO/-NH+）驱动n（PFB⁃MC）/n（1⁃OA）=1∶1的复配体系，可形成致密单分子膜，显著降低表面张力至28.0 mN/m（降幅>30%）；n（PFB⁃MC）/n（1⁃OA）=1∶1的复配体系中抗衡离子实现了均匀分散，且1⁃OA具有强氢键能力；通过动态可视化技术可揭示非等物质的量比下抗衡离子的竞争吸附行为，且可阐明比例调控机制。该研究可为环保短氟碳灭火剂设计提供原子尺度依据，指导复配比例优化。

A.SHARMA等采用全原子分子动力学模拟（OPLS⁃AA力场），系统研究了1～4头阳离子表面活性剂（吡啶盐类）在水溶液中的自组装行为，揭示了多头结构对聚集行为的调控机制。结果表明，随头基数目增加，静电排斥增强，导致聚集数显著降低，且聚集体形状偏离球形；4头表面活性剂因头基灵活性较高，可通过构象调整减弱静电排斥，形成较大聚集体（最大聚集数达15）。该研究从原子层次明确了阴/阳离子静电互补是驱动多头表面活性剂自组装及界面行为的主导因素，为理性设计离子型表面活性剂体系提供了关键依据。