超越脱水：局域配位塑造聚电解质膜中的离子传输势垒

《Journal of Membrane Science》：Beyond dehydration: Local coordination shapes ion-transport barriers in polyelectrolyte membranes

　　

纳滤膜中的离子传输势垒常常以位阻排斥、Donnan 排斥和脱水来解释，然而在非均相带电聚合物网络中，其分子起源仍然难以解析。本文将分子动力学（MD）自由能计算与机器学习（ML）相结合，以识别控制聚（苯乙烯磺酸盐）/聚（二烯丙基二甲基氯化铵）聚电解质膜中传输势垒的局域离子环境描述符。研究人员采用伞形采样模拟重建了钠离子和氯离子传输的自由能剖面，同时通过第一与第二配位壳层的配位特征以及水合壳层几何结构表征局域离子环境。两种离子在膜内均表现出显著势垒，并伴随部分脱水、水合壳层尺寸减小以及壳层畸变增加。然而，这些仅基于水合的描述符本身并不能唯一决定局域自由能。值得注意的是，当水分子损失可由膜官能团及邻近离子的配位加以补偿时，某些弱水合构型仍可表现出相对有利的自由能。结合 SHAP 分析的 XGBoost 模型表明，在所采样的渗透路径内，局域自由能惩罚由水合、位阻限域、离子配对及聚合物配位的耦合作用共同塑造。邻近反离子以及带相反电荷的膜基团在降低与脱水相关的自由能代价方面似乎尤为重要。这些发现为聚电解质膜中离子传输势垒提供了具有路径分辨率的分子层面解释，并进一步表明，MD–ML 联合分析能够揭示控制非均相带电膜环境中自由能变化的主要描述符。

该论文发表于《Journal of Membrane Science》，聚焦于聚电解质膜中离子传输势垒的分子本质这一膜科学核心问题。当前对反渗透（RO）与纳滤（NF）膜分离机理的理解虽已显著进展，但针对离子为何在尺寸和电荷相近条件下仍表现出不同迁移选择性，传统的位阻排斥、Donnan 排斥与脱水理论仍不能给出充分解释。尤其在聚电解质多层膜这类高度非均相、强带电的聚合物网络中，局域水合状态、聚合物官能团分布、离子配对以及纳米尺度限域效应彼此耦合，使得离子在膜内经历的自由能势垒难以仅由单一参数刻画。实验手段又难以在膜孔内部同时提供足够的时空分辨率，因此需要借助分子尺度计算框架解析离子跨膜过程中的过渡态特征及其局域环境起源。正是在这一背景下，研究人员开展了本研究，旨在超越“脱水决定势垒”的单因素认识，解析局域配位环境如何共同塑造聚电解质膜中的离子传输障碍。

研究人员以聚（苯乙烯磺酸盐）/聚（二烯丙基二甲基氯化铵）膜（PSS/PDADMAC PEM）为模型体系，系统考察 Na⁺ 与 Cl^? 在膜内的自由能剖面及其局域环境描述符之间的关系。研究结论表明：离子进入 PEM 内部时普遍遭遇显著自由能势垒，并伴随部分脱水、水合壳层横截面积减小和各向异性增强；但局域自由能并非由脱水程度唯一决定，而是由水合重构、位阻限域、离子配对与聚合物配位共同耦合决定。特别是当失去的水分子可被带相反电荷的膜官能团或邻近反离子所补偿时，即使离子处于弱水合状态，其局域自由能仍可相对较低。这一结论的重要意义在于，它从路径分辨的分子层面解释了聚电解质膜中的离子传输势垒，为理解化学非均相膜中的离子选择性提供了超越体相水合热力学的新框架，也显示出分子动力学（MD）与机器学习（ML）联合分析在解析复杂膜环境结构—能量关系方面的潜力。

在技术方法上，研究人员首先构建了全原子 PSS/PDADMAC 聚电解质多层膜模型，在 0.168 M NaCl 条件下，对 4 个独立制备的膜构型开展模拟。随后通过牵引分子动力学（SMD）生成跨膜离子渗透路径，并利用伞形采样（US）沿膜法向坐标重建平均力势（PMF）自由能剖面，采用多状态 Bennett 接受比（MBAR）估算自由能及不确定性。围绕离子第一、第二水合壳层，提取水分子、离子、聚合物官能团配位数、水合壳层横截面积及各向异性等描述符，并以分块平均策略获得与 MBAR 一致的统计量。最后分别针对 Na⁺ 与 Cl^? 建立 XGBoost 回归模型，并结合 SHAP（Shapley additive explanations，沙普利加性解释）识别对局域自由能最具预测性的环境特征。

在研究结果部分，论文首先在“3.1. Free-energy barriers for ion permeation through the PEM”中表明，研究人员通过对单条渗透路径的一维自由能剖面重建，量化了离子由水相进入 PEM 的能量代价。不同路径的自由能曲线表现出明显差异，说明在非均相带电膜中，聚合物密度、电荷分布、水合程度及可移动离子含量的局域变化会形成彼此不同的渗透环境。所有离子进入膜后自由能总体升高，反映出其从体相水转移至致密聚合物环境所需付出的显著代价。研究人员进一步对不同路径的最大势垒进行 Boltzmann 加权，得到有效自由能势垒，其量级与文献中常规 NF/RO 膜的单价离子传输势垒相近，说明该模拟所得能量景观具有物理合理性。

在“3.2. Ion coordination behavior”部分，研究人员分析了离子由溶液进入膜中央区域后第一与第二配位壳层的环境变化。结果显示，Na⁺ 和 Cl^? 在第一壳层中均发生部分脱水，其中 Na⁺ 的失水更明显，而 Cl^? 因第一壳层半径较小，失水程度相对较低。与此同时，离子与膜官能团的配位增加，且呈现离子特异性：Na⁺ 更多与 PDADMAC 的胺邻位 CH₂/CH₃ 基团及部分 PSS 磺酸根氧发生配位；Cl^? 则主要与 PSS 的磺酸根氧相互作用。第二壳层中，两种离子的失水更为显著，尤其 Na⁺ 表现出更大程度的水配位下降及相应更强的膜官能团补偿配位。此外，研究人员还通过水合壳层面积和各向异性分析发现，离子进入 PEM 后壳层面积减小、形状由近球形转向更拉长的畸变构型，说明跨膜传输不仅伴随化学环境重排，还伴随显著的几何重构。

在“3.3. Free energy dependence on the local ion environment”部分，研究人员进一步将局域自由能与局域环境描述符联系起来。分箱分析表明，无论是第一还是第二配位壳层，水配位数降低、水合壳层面积缩小及各向异性增加，总体上都与更高的局域自由能惩罚相关；且第二壳层中的相关性比第一壳层更强，说明离子周围更外层的移动相重排对膜内能量学同样关键。然而，各类关系都伴随较大标准差，显示没有任何单一描述符能够充分解释全部自由能波动。研究人员随后选取一个反直觉子集：在 PEM 内既处于低水配位、又具有低自由能的 Na⁺ 构型。分析显示，这些弱水合但能量相对有利的构型具有更高的 PDADMAC 胺邻位 CH₂/CH₃ 配位，同时减少了与带负电 PSS 环境的不利关联。该结果直接证明，脱水所导致的能量代价可以被电荷互补的聚合物接触和局域离子配位所补偿，因此“弱水合”并不必然等价于“高势垒”。

在“3.4. Combining local descriptors with an exploratory ML model”部分，研究人员利用 XGBoost 模型整合第一、第二壳层的全部配位与几何描述符，对 Na⁺ 与 Cl^? 的局域自由能进行多变量诊断。模型结果显示，完整描述符集合比仅使用第一壳层水配位和壳层面积的简化模型预测效果更好，说明自由能变化确实不能简化为纯水合效应。SHAP 分析进一步揭示，对 Na⁺ 而言，第二壳层 PSS 苯基配位增多通常对应更高预测自由能，而第一壳层 Cl^? 配位增多则与更低预测自由能相关，体现了有利离子配对的稳定化作用；第二壳层 PDADMAC 相关配位在较多情况下也降低预测自由能。依赖关系图显示，某些有利的聚合物或反离子配位即使伴随局域失水，仍与较低自由能相关。由此可见，模型支持这样一种解释：在所采样路径中，真正支配离子局域自由能的不是单独的脱水程度，而是脱水、位阻限域、离子配对和聚合物配位之间的耦合平衡。

讨论部分表明，该研究将传统上强调的脱水屏障概念推进到了更精细的局域配位视角。研究结果并未否定脱水在离子跨膜能量学中的重要性，而是指出其热力学代价取决于失去的水配位是否能够被有利的局域相互作用补偿。在化学组成高度非均相的聚电解质膜中，膜官能团的电荷性质、空间接近性以及邻近反离子的存在，可显著改变局域自由能景观。因此，对于具备离子特异性补偿相互作用的膜体系，仅凭位阻排斥或体相水合热力学难以充分解释传输势垒。论文同时也强调，当前 ML 模型带有明显路径依赖性，说明更广泛的采样与更全面的描述符仍是建立一般化结构—能量关系所必需的。

研究结论部分可译为：研究人员采用伞形采样分子动力学模拟并结合探索性机器学习分析，考察了 PSS/PDADMAC 聚电解质膜中局域离子环境描述符与传输自由能势垒之间的关系。Na⁺ 和 Cl^? 在膜内均遭遇显著势垒，这些势垒平均上伴随着部分脱水、水合壳层面积减小及壳层各向异性增强。然而，局域自由能与基于水合的描述符之间的成对相关性存在较大标准差，表明仅凭水合不能唯一决定局域自由能惩罚。相反，失水的能量代价强烈依赖于更广泛的局域配位环境。当脱水可被膜官能团或邻近离子的有利配位所补偿时，弱水合的 Na⁺ 构型仍可能具有较低自由能。XGBoost 与 SHAP 分析从多变量角度支持了这一点：对于 Na⁺ 模型，第一壳层 Cl^? 配位和第二壳层胺邻位 CH₂ 配位增加，即使伴随较低水配位，也与较低预测自由能相关。因此，水分子损失并非在所有局域环境中都天然不利，其能量效应取决于被置换水分子周围是否存在电荷互补的聚合物接触或邻近反离子，抑或相反存在较不利、拥挤的聚合物环境。总体而言，传输势垒来源于脱水、位阻限域、离子配对和聚合物配位的耦合作用，而非任何单一描述符的独立控制。该研究说明，MD–ML 联合框架可用于解析非均相带电膜中的局域自由能环境，同时也强调需通过更广泛采样建立更具普适性的结构—能量关系。