水凝胶，最新ACS Nano：动态电荷调节让水凝胶离子传输“加速”

弱聚电解质水凝胶由化学或物理交联的聚电解质链构成，因其溶胀能力及对盐、温度、光等刺激的响应性，被视为脱盐、电黏附、药物递送和纳米力学等领域的可调材料平台。许多常用水凝胶由弱聚电解质构成，其电离度依赖pH，并会随盐浓度、离子价态及静电相互作用动态调整。这种电荷调节会影响水凝胶的pH响应溶胀、与客体分子的相互作用，甚至导致一级溶胀转变；但现有理论研究多集中在静态平衡性质或大尺度构象运动上，对小（多价）离子和带电客体分子如何穿过电荷调节聚电解质网络仍缺乏机制性理解。困难在于，电荷调节与盐浓度、离子结合、聚合物密度和网络结构深度耦合，实验难以解耦，全原子模拟成本过高，解析方法又常需大幅近似。为此，作者采用粗粒化、隐式溶剂分子动力学模拟结合反应系综蒙特卡洛方法，并通过逐步增加体系复杂度的分层策略，分离并考察各个效应。

加州理工学院Pierre J. Walker博士和斯图加特大学David Beyer博士合作，研究人员利用大规模分子模拟剖析了弱聚电解质水凝胶中小离子传输的耦合机制，揭示电荷密度、竞争结合和动态电荷调节如何共同控制离子迁移率。通过系统隔离这些效应，研究发现多价离子的自扩散对水凝胶pKa和盐浓度呈现显著非单调依赖；而一旦允许水凝胶骨架动态调节其电荷，这种非单调性消失，三价离子的扩散系数可提高2倍以上。作者将这一增强归因于“质子化辅助跳跃”机制：骨架持续重塑局部静电景观，在平移逃逸之外提供额外的分子尺度解绑路径，加速离子对动力学，却不改变自由离子比例。结果表明，动态纳米尺度电荷调节能够将离子迁移率与平衡结合强度解耦，为固定电荷材料之外提供新的设计杠杆，并据此提出下一代离子导水凝胶的分子设计原则。相关论文以“Ion Transport Enhancement in Charge-Regulating Hydrogels”为题，发表在ACS Nano上。

为厘清上述耦合效应，作者构建了类金刚石晶胞水凝胶（图1a），并设计三类模拟逐层分离关键因素：无盐电荷调节模拟结合分子动力学与反应系综蒙特卡洛获得平衡电荷构型，随后冻结骨架电荷并加入多价盐，以隔离电荷密度效应；加盐电荷调节模拟同时保留电荷密度、竞争结合和动态骨架效应；加盐竞争结合模拟则从电荷调节轨迹中取平衡电荷构型并冻结骨架电荷，保留竞争结合但去除动态电荷调节（图1b）。由此，作者得以在超过700个体系中逐步增加复杂度，考察电荷密度、竞争结合和动态电荷调节三种相互竞争的效应（图1c）。

图1：方法。(a) 使用类金刚石晶胞的模拟水凝胶结构。(b) 本工作中进行的模拟类型、所用方法及相应标签的总结。(c) 本工作研究的不同竞争效应：1. 水凝胶电荷密度变化；2. 添加盐诱导的竞争结合；3. 水凝胶骨架的动态电荷调节。

在电荷密度方面，无盐条件下随pKa升高，水凝胶平衡电离度下降；相邻电离珠之间的强静电排斥会抑制骨架电离，造成表观pKa偏移，即弱聚电解质的聚电解质效应（图2a）。冻结平衡电荷构型并加入多价盐后，阳离子自扩散系数总体随阳离子价态升高和电离度升高而下降，较高盐浓度可部分恢复迁移率；在高价阳离子、低盐和低pKa条件下出现明显扩散最小值（图2b）。该最小值并不出现在未结合离子分数中，而可用离子对弛豫时间解释：二价和三价体系在低盐下于中间pKa出现弛豫时间最大值，与扩散最小值对应（图2c）；高pKa、低盐时阳离子强结合于骨架，运动受慢骨架耦合；降低pKa使更多带电位点可用，阳离子更易沿骨架平移，扩散增加（图2d）；若电荷沿骨架均匀分布或盐浓度足够高，最小值消失（图2e）。

图2：电荷密度。(a) 无盐条件下水凝胶电离度随pKa的变化；理想极限预测以黑线显示；还包含电荷密度模拟中沿线性链段典型电荷分布的示意图。(b) 使用电荷密度模拟获得的阳离子自扩散系数随pKa、阳离子价态和盐浓度的变化；颜色表示价态（蓝：一价，黄：二价，绿：三价），线型表示浓度（实线：0.05 M，虚线：0.01 M）；点线表示完全电离水凝胶的扩散系数。(c) 在与(b)相同条件下，添加阳离子与水凝胶骨架之间的离子对弛豫时间，使用相同颜色和线型。(d) 静电关联对阳离子与水凝胶骨架结合影响的示意图。(e) 二价阳离子情况下阳离子自扩散系数的盐浓度依赖性随pKa的变化。

进一步考虑竞争结合后，加入盐会通过阳离子与骨架之间的有利静电相互作用使去质子化平衡前移，提高电离度；一价盐影响很小，二价尤其是三价离子可显著移动平衡，三价离子在高浓度下甚至使电离度超过理想溶液极限，有效逆转聚电解质效应（图3a）。为隔离盐诱导电离偏移而不引入动态电荷调节，作者冻结平衡电荷分布后进行竞争结合模拟。结果显示，扩散随pKa的定性依赖与电荷密度模拟相同，但拐点发生移动，并在低盐下出现更深的扩散最小值；在相同电离度下，竞争结合模拟的扩散系数持续低于电荷密度模拟（图3b、图3c），说明骨架电荷分布发生了结构重组，增强了离子对寿命并抑制阳离子扩散（图3d）。电荷分布自相关分析表明，多价盐不仅移动整体电离度，还使电荷沿骨架重新组织，形成可被单个多价阳离子桥接的交替电离模式，增加阳离子配位数并形成阻碍沿骨架运动的结合环境（图3e）。

图3：竞争结合。(a) 在添加盐浓度为0.05 M时水凝胶电离度随pKa的变化；包括理想极限预测；还包含电荷密度和竞争结合模拟沿线性链段电荷分布差异的示意图。(b,c) 使用电荷密度（细线）和竞争结合（粗线）模拟获得的阳离子自扩散系数随pKa (b) 和电离度 (c)、阳离子价态及盐浓度的变化；颜色表示价态（黄：二价，绿：三价），线型表示浓度（实线：0.05 M，虚线：0.01 M）；点线表示完全电离水凝胶的扩散系数。(d) 与(b)相同条件下添加阳离子与水凝胶骨架之间的离子对弛豫时间，使用相同格式。(e) pKa为-1时，不同单体间隔下的电荷分布自相关，使用与(b)相同方案。(f) 二价阳离子情况下电荷密度与竞争结合模拟中阳离子自扩散系数的盐浓度依赖性随pKa的变化。

当允许骨架位点在模拟过程中动态调节电荷状态时，所有考察条件下电荷调节模拟给出的自扩散系数均高于平均电离度相同的竞争结合体系，低pKa时阳离子迁移率甚至可超过电荷密度更低的固定电荷模拟；增强幅度强烈依赖价态：二价离子仅小幅增加，三价离子可加速2倍以上；更引人注意的是，电荷调节模拟中未结合离子比例与相应固定电荷体系相当，但此前的非单调扩散最小值消失（图4a）。其来源是更短的离子对弛豫时间（图4b）：离子对既可通过平移逃逸解离，也可通过配位骨架位点在反应系综蒙特卡洛移动中发生质子化而解离；后者一旦发生，预期比平移逃逸更快。作者将此解释为“质子化辅助跳跃”过程：某带电位点要具有可观质子化概率，附近须有带电且未被占据的单体，一旦质子化发生，多价阳离子可迅速迁移至相邻位点并重新结合，从而绕过直接质子化配位多价离子位点的高能量惩罚，加速结合与解绑动力学（图4d）。改变蒙特卡洛移动与分子动力学步数比例两个数量级后，离子对弛豫时间和自扩散系数的变化远小于速率变化本身，说明该跳跃机制并不强烈受本征质子化速率限制（图4c）。

图4：电荷调节。(a) 使用电荷调节（粗线）、竞争结合和电荷密度（细线）模拟获得的阳离子自扩散系数随pKa、阳离子价态和盐浓度的变化；颜色表示价态（蓝：一价，黄：二价，绿：三价）；点线表示完全电离水凝胶的极限扩散系数。(b) 与上述相同条件下添加阳离子与水凝胶骨架之间的离子对弛豫时间，使用相同颜色。(c) 改变蒙特卡洛移动与分子动力学步数比例对电荷调节模拟所得自扩散系数的影响。(d) 质子化辅助跳跃机制增强离子传输的直观示意。(e) 二价阳离子在pKa为-1时自扩散系数的盐浓度依赖性在三种模拟类型之间的变化。

基于上述认识，作者提出了若干假设性应用启示。对于水凝胶离子膜，若平均电荷分数相同，电荷调节水凝胶的离子分配预计与固定电荷等效体系相当，但增强的扩散系数应直接转化为更高渗透性，这一优势对高价离子尤其明显，可超过2倍，因而与重金属分离等目标离子常具高价态的场景相关（图5a）。对于共聚物水凝胶网络，作者提出在弱聚电解质水凝胶中引入中性间隔单体，可降低竞争结合影响，增加链间交联比例，并改善离子传输和弹性模量等力学性能（图5b）。对于生物凝聚体客体包封，作者假设电荷调节凝聚体的自适应环境可改善客体分子摄取并提高其扩散系数（图5c）。

图5：电荷调节对离子传输可能影响的假设。(a) 水凝胶离子膜：我们预期对于相同电荷分布，电荷调节水凝胶应具有与固定电荷等效体系相同的离子分配，但扩散增强，因而渗透性更高。(b) 共聚物水凝胶网络：我们提出在弱聚电解质水凝胶中引入中性间隔单体应减少竞争结合影响，增加链间交联比例，并改善离子传输和凝胶力学性能如弹性模量。(c) 生物凝聚体客体包封：我们假设电荷调节生物凝聚体的自适应环境可改善客体分子摄取并提高其扩散系数。

小结，作者用粗粒化模拟解构了电荷调节水凝胶中离子传输的物理机制：小离子自扩散由三个相互耦合因素竞争决定，即静电结合阻碍离子运动、静电屏蔽促进传输、电荷调节本身则视条件抑制或增强迁移率。动态电荷调节可通过质子化辅助跳跃消除非单调扩散最小值，并将离子迁移率与平衡结合强度解耦，为固定电荷材料之外提供新的设计杠杆。作者据此提出面向离子选择性膜、共聚物网络、温度依赖和生物分子凝聚体等方向的未来实验与理论研究思路，并强调将电荷调节与聚合物网络结构控制相结合，可为膜科学、能源材料、软物质力学和生物物理中的下一代离子导材料设计提供强大框架。