返回周柯个人主页
近期工作 · J. Phys. Chem. Lett.

限域离子的结晶和溶解——限域调节离子水合状态构建高性能力敏纳米通道

发布日期:2026 年 1 月 30 日 研究团队:苏州大学周柯

该研究揭示,单层水极端限域下通道高度会改变离子水合与离子-离子作用的竞争,使离子在自由载流态与二维团簇态之间切换,并将机械压力转换为离子电流的百倍开关。

外力调节限域离子通道并构建力敏离子晶体管和逻辑门
图形摘要:外力改变二维通道高度 \(h\),使离子从水合导电态转为团簇绝缘态,形成力敏离子晶体管和机械 NOT 门。
J. Phys. Chem. Lett. 论文首页,已裁去下半页正文
论文首页,保留期刊标识、标题、作者、日期和摘要区域;点击图片可前往论文页面。

从构象门控到载流子门控:力敏通道的新路径

生物体中的 Piezo、TREK 和 TRAAK 等机械敏感离子通道,通常通过膜蛋白构象变化调节离子通过。人工纳米孔道若只依赖几何形变,通道尺寸的微小变化很难带来大的电流调制。论文提出另一条力学路径:不直接改变离子通道的几何通量,而是用机械力控制“真正能够导电的自由离子”数量 \(n_c\)。

把自由离子数 \(n_f\) 看作反应坐标,对其自由能面进行配分加权,可得到有效载流子数:

\[ n_c=\frac{\displaystyle\int_0^N n_f\exp\!\left[-\frac{\Delta F_f(n_f)}{k_{\mathrm B}T}\right]\,\mathrm dn_f} {\displaystyle\int_0^N\exp\!\left[-\frac{\Delta F_f(n_f)}{k_{\mathrm B}T}\right]\,\mathrm dn_f}, \qquad I\propto n_c \]

其中,成对束缚的阳离子和阴离子对净电荷贡献近似抵消,或者因尺寸变大而迁移率大幅下降,因此 \(n_c\) 比总离子数更接近实际的导电自由度。

极端限域:自由离子为何变成二维岩盐团簇

研究在 \(h=0.8\,\mathrm{nm}\) 的单层水通道中放入 30 对 NaCl 和 600 个水分子,并使用经过第一性原理修正的力场与微秒尺度增强采样分子动力学。即使浓度只有 \(2.4\,\mathrm{mol\,kg^{-1}}\),离子也会在约 100 ns 内从分散状态逐渐组装为二维岩盐样团簇;LiCl、KCl 和 RbCl 也表现出相同的趋势。

论文用水合数 \(S_H\) 和反离子局部配位数 \(S_O\) 描述离子状态。由增强采样得到的二维自由能面满足:

\[ F(S_O,S_H)=-k_{\mathrm B}T\ln P(S_O,S_H)+C \]

在体相水中,自由能最低点靠近 \(S_O\to0\),对应充分水合的自由离子;在单层水中,最低点转移到 \(S_H\to0\),对应水合壳层塌缩和高度有序的离子团簇。这是“限域越强,离子反而越容易结晶”的核心反直觉现象。

单层限域水中离子团簇化演化和水合自由能面
图 1. 离子由分散态演化为二维岩盐团簇,以及单层水与体相水中 \(S_O-S_H\) 自由能面的反转。

水合与离子-离子作用的热力学竞争

离子结晶不是单一作用力造成的,而是水合、离子-离子吸引和熵损失共同竞争的结果。将聚集过程中的焓变化拆开,可写成:

\[ \Delta H=\Delta H_{i-i}+\Delta H_{\mathrm{hy}}, \qquad \Delta H_{\mathrm{hy}}=\Delta H_{i-w}+\Delta H_{w-w} \]

在单层水中,限域增强的 \(\Delta H_{i-i}\) 下降幅度超过水合项的代价,驱动体系从自由水合态进入束缚团簇态;当 \(h\) 增大到双层水及以上,水合效应重新占优。对单个离子对的平均力势 \(W(d)\) 而言,接触离子对(CIP)和溶剂共享/分离离子对(SSIP)之间的势垒可表示为:

\[ W(d)=-k_{\mathrm B}T\ln g_{+-}(d)+C, \qquad \Delta W_{c-s}=W(d_{\mathrm{SSIP}})-W(d_{\mathrm{CIP}}) \]

当 \(h\lesssim8\,\text{Å}\) 时,\(\Delta W_{c-s}\) 可达到约 \(10k_{\mathrm B}T\),使已经形成的离子对难以重新解离;即便势垒达到 \(5k_{\mathrm B}T\),由 \(1-\exp[-\Delta W_{c-s}/(k_{\mathrm B}T)]\) 估算的阻滞率也已超过 99%。

埃级高度窗口:成核力学放大限域效应

自由能面显示,\(h\lt9.5\,\text{Å}\) 时 \(\Delta F_f(n_f)\) 近似单调下降,体系偏好 \(n_f\approx0\) 的完全团簇态;当 \(h\ge9.5\,\text{Å}\) 时,自由能曲线恢复体相水的抛物形状。\(h\approx8.5\,\text{Å}\) 附近出现最强的团簇稳定性,原因是单层水厚度、离子尺寸和二维盐晶格间距形成几何相容。

把有效可用宽度记为 \(h^*=h-3.4\,\text{Å}\),平面限域下的经典成核自由能可写成:

\[ \Delta G(r)\sim-4\pi r^2h^*\Delta\mu+2\pi rh^*\sigma, \qquad h^*=h-3.4\,\text{Å} \]

因此,离子团簇的自由能代价同时受离子对有效能量密度和可用厚度控制,论文得到现象学标度:

\[ \Delta F_f(N)\sim\Delta W_{c-s}h^*, \qquad \frac{\Delta F_f(N)}{Nh^*}\sim\Delta W_{c-s} \]
自由离子自由能曲线和有效载流子数随通道高度的开关行为
图 2. 自由离子自由能曲线随 \(h\) 的变化,以及有效载流子数在埃级通道高度上的跃迁。

从机械压力到离子 NOT 门

通道高度只需从 \(10.0\,\text{Å}\) 压缩到 \(9.0\,\text{Å}\),自由离子数就会从有限值降至近零,离子流随之关闭。由于 \(I\propto n_c\),开关比主要由载流子密度而非几何截面积决定:

\[ R=\frac{I_{\mathrm{on}}}{I_{\mathrm{off}}} \approx\frac{n_{c,\mathrm{on}}}{n_{c,\mathrm{off}}}\sim10^2, \qquad 10.0\,\text{Å}\xrightarrow{\,\text{pressure}\,}9.0\,\text{Å} \]

模拟中开放态电流约为 \(1.25\,\text{nA}\),受压后的团簇态接近 \(0\);若只考虑孔道几何缩小,面积比例带来的开关比约为 \(1.18\),远小于团簇相变提供的约 100 倍放大。机械输入“pressed/free”与离子输出“off/on”构成反向响应,因此可以直接实现离子机械 NOT 门,进一步通过串并联组合扩展到 NAND 和 NOR。

力敏离子通道的开放态、关闭态和机械 NOT 逻辑门示意
图 3. 机械压力改变限域高度,驱动离子结晶/溶解并控制电流,实现力敏离子晶体管和机械 NOT 门。

这项工作把“力改变孔道形状”的传统力敏图景推进为“力改变离子相态和载流子密度”:限域高度、离子水合、二维成核和电流门控由同一条自由能链条连接起来。

论文信息

Ke Zhou*, “Ion Clustering Regulated by Extreme Nanoconfinement Enables Mechanosensitive Nanochannels,” J. Phys. Chem. Lett. 17, 1798–1804 (2026). https://doi.org/10.1021/acs.jpclett.5c03717