4张红外光谱图，发了一篇Nature Chemistry！

破解水合质子红外连续谱之谜

水合质子——简单来说，就是一个带正电的氢离子（质子）掉进了水里——是自然界中最基本的化学物种之一。往一杯水里滴入一滴盐酸，酸味来自于它；我们身体里每一个细胞代谢能量时，质子穿过生物膜的过程仰仗于它；甚至氢燃料电池和下一代储能电池的高效运转，也离不开质子在液态水中的快速迁移。可以说，它几乎渗透于化学、生物学与能源科学的每一个角落。

然而，这个看似简单的物种却让科学家们困惑了整整两个多世纪。早在1806年，德国化学家Grotthuss就提出了质子在水中通过“跳跃”方式在分子间传递的经典模型，但直到今天，关于水合质子究竟长什么样子、在水里如何运动，学界仍然没有统一的答案。最令人头疼的是它的红外光谱——如果把红外光打到酸性水溶液上，得到的吸收信号不是几条清晰的谱线，而是一片横跨约1,000至3,500 cm⁻¹的“大抹带”，仿佛一团迷雾笼罩了整个中红外区间。

为什么这么难？因为水分子本身就在不停地热运动、断裂和重组氢键，质子夹在中间，时而偏向一边，时而又被拉向另一边。科学家们早先提出了两个“定格”的模型来帮助理解（图1b）：一个是Zundel阳离子（H₂O₂⁺），想象质子像一颗三明治里的夹心，被左右两个水分子对称地夹在正中间；另一个是Eigen阳离子（H₉O₄⁺），想象质子成了一个小团伙的核心——中心是一个水合氢离子（H₃O⁺），周围再围上三个水分子作为“跟班”。这两个模型各自在气相实验中都表现良好，但放到液态水中，它们的简单叠加却完全解释不了那片“大抹带”式的连续光谱。原因很简单：液态水里的质子根本不“安分守己”，它在这两种结构之间来回切换，速度快到飞起。要理解这片连续谱，就必须超越“静态快照”的思维，在动力学层面上追踪质子如何随局域环境的涨落而实时演化。

如今，德国海德堡大学Oriol Vendrell教授、David Mendive-Tapia博士和英国剑桥大学Christoph Schran教授、法国勃艮第大学Dominik Marx教授合作，报告了一种全新的“频闪电影”策略——用计算模拟的手段，像放慢动作电影一样，一帧一帧地追踪质子在局域溶剂环境变化时，其红外响应如何连续演变（图1a）。他们精心挑选了一个“恰到好处”的模型体系：扩展Zundel配合物（eZ，H⁺(H₂O)₆），这个由六个水分子和一个质子组成的小团簇，恰好同时“内含”了Zundel和Eigen两种核心结构，是架起气态模型与液态真实环境之间桥梁的理想窗口（图1b）。他们发现，Zundel核心周围的第二层溶剂水分子——那些不直接接触质子、但像“旁观者”一样在第二圈排列的水分子——扮演着一位“隐形导演”的角色：当这层溶剂壳层排列得高度对称时，Zundel核心原本特征性的双峰（一种叫“费米共振”的量子力学借能量效应）被彻底压制，原本清晰的“双胞胎”峰变成了一团模糊的宽包；而一旦打破这种对称性——比如把第二壳层中的一个水分子稍稍拉远——光谱就会随之平滑地演变，从对称Zundel的模样逐步过渡到非对称的Zundel-like特征，最终呈现出类似Eigen的谱学面貌（图1a）。这项工作为理解液态水中质子的连续红外吸收提供了第一张“结构-光谱”之间的精准对应图谱。相关论文以“Deciphering the infrared spectrum of the hydrated proton using full-dimensional quantum dynamics”为题，发表在Nature Chemistry上。

图1 | 水合质子在扩展Zundel结构中的局域溶剂环境“频闪电影”及其红外响应随非对称性的变化。 a，从左到右随溶剂化非对称性增加，与eZ中水合质子相关的代表性结构示意图，以及中心氢键红外响应的频闪序列。计算所得中心质子的严格量子红外响应（见正文及图4说明）覆盖超过1,500 cm⁻¹的宽度，其随非对称性的变化解释了水合质子在真实涨落环境中极端宽泛的红外响应。b，各极限基元的分子结构：（从左至右）Zundel（Z）、扩展Zundel（eZ）和Eigen（E）配合物。在eZ配合物中，将Zundel核心的两个水分子（H₂O₂，橙色矩形）称为“内层”水分子，其余四个单体称为“外层”水分子。eZ配合物也可被视为被两个水分子溶剂化的E核心（H₉O₄⁺，绿色三角形），这两个水分子打破了裸E阳离子的三重对称性；这种非对称E-like构型源于过量质子沿eZ中中心氢键的完全转移（该质子也被称为“最活跃质子”），从而形成新生水合氢离子核心H₃O⁺。

研究人员以eZ配合物为起点，首先精确计算了它在气相中的红外吸收光谱，并与实验数据进行了严格的“面对面”比对（图3a）。一个值得注意的细节是，理论光谱没有做任何“调频”修正（很多计算工作为了贴合实验会额外平移频率），而它竟然在整个光谱范围内与实验信号高度吻合——这本身就说明计算所用的势能面和量子动力学方法是相当可靠的。

通过对光谱进行“拆解”——把每一个振动峰追溯到具体的原子运动模式——团队划分出了五个特征鲜明的光谱区域（图3b）。在低频段（区域I，低于800 cm⁻¹），信号主要来自外层水分子的“扭腰”和“摇摆”等慢速运动，以及外层氧-氧骨架的伸缩；而到了约1,000 cm⁻¹附近（区域II），出现的强吸收峰则明确指向质子沿着中心氢键“来回穿梭”的转移运动——你可以把它想象成质子像一颗弹珠在左右两个水分子之间来回振动。真正有趣的是1,600–1,800 cm⁻¹区间（区域III）：这里出现了质子转移运动与内层水分子“非对称弯曲”（类似于水分子以不等角度“歪头”）的强烈耦合——一个运动一发生，另一个就被“带节奏”。这种耦合在后续非对称溶剂环境中，恰恰成为塑造宽谱连续带的主要驱动力。而旁边1,618 cm⁻¹处一个相对尖锐的峰，则被证实与质子的转移运动几乎无关——它仅仅是第二壳层水分子的纯弯曲振动，“碰巧”出现在这个频率附近而已。

转向高频段（区域IV和V），内层Zundel核心的O-H键伸缩振动呈现出一个宽达500 cm⁻¹的“大鼓包”（图3b），其宽度和Eigen配合物中水合氢离子的伸缩带差不多。为什么这么宽？因为内层O-H振动与整个团簇骨架的“变形”运动（非谐角位移）广泛地“纠缠”在一起，相当于你拉了其中一根弹簧，整个架子都在晃。这种宽化意味着连接内外层水分子的那些氢键非常“灵活好动”——如果把这个团簇放进体相水里，它们随时可能因为热涨落一跃成为那个“最活跃的氢键”，也就是质子实际跳跃的通道。相比之下，外层水分子的O-H伸缩峰形要尖锐得多，因为这些“悬挂”在外面的氢原子与团簇内部其他部分的耦合较弱，相对“与世无争”。值得一提的是，a₀谱带的理论宽度比实验略宽一些，研究者认为这是因为外层水分子中那些对称的伸展模式与团簇骨架之间的关联比预期更强，导致理论谱线被额外展宽了——这并不影响主要结论，反而揭示了团簇内部更精细的相互作用。

图3 | 气态扩展Zundel配合物的红外光谱。 a，计算所得线性红外吸收光谱（黑色）和质子化水六聚体H⁺(H₂O)₆的IR³MS²实验振动光谱（蓝色）；灰色线条遮蔽了基于参考文献61中氩标记光谱所报道的2,000–3,000 cm⁻¹平坦基线所对应的分辨率不足的频率窗口。模拟光谱的标签用罗马数字（I、II、…）标示主要谱带组，用字母（a₁、a₂、…）标示相应的子带。b，使用局域振动激发对（灰色）全谱进行分解，并按振动运动对最相关峰进行归属（颜色，插图），详见正文。

接下来，一个尤其令研究者着迷的现象出现了。将eZ配合物的光谱（图2c）与裸Zundel阳离子的光谱（图2a）放在一起对比，一个差异“跳”了出来：裸Zundel在约1,000 cm⁻¹处那对标志性的“双胞胎”峰（费米共振双峰），在eZ中居然完全消失了，取而代之的是一个宽而无结构的单峰（比较图2a与图2c）。这种费米共振到底是怎么回事？通俗地说，就是质子的伸缩振动（基频）和另一个低频振动（水分子的摇摆加O-O伸缩）的倍频恰好撞到了几乎相同的频率——就像两个人的脚步恰好合拍，于是发生了“能量借入”，产生了一对强度相当的双峰。那么，为什么在eZ中这对“合拍”的搭档被拆散了呢？

图2 | 三种具有递增复杂局域溶剂环境的质子化水团簇的实验红外光谱。 a–c，来自参考文献58的裸Zundel阳离子（Z: H₂O₂⁺或H⁺(H₂O)₂）（a）、来自参考文献59、60的裸Eigen阳离子（E: H₉O₄⁺或H⁺(H₂O)₄）（b）以及来自参考文献51的扩展Zundel配合物（eZ: H₆O₂⁺或H⁺(H₂O)₆）（c）的实验红外光谱。d，酸性溶液中对应水合质子H⁺(aq)贡献的实验差谱红外光谱，来自不同浓度下应用不同扣除方案的结果（见图例）；彩色频率窗口（黄色、蓝色、绿色、紫色）作为视觉引导，用于将此面板与图4中面板b和c进行对齐和比较。

研究者做了一个简单而巧妙的实验：将Zundel核心周围那四个“围观”的水分子，以完全对称的方式、像剥洋葱一样一层层慢慢拉远（图4a）。结果发现，从2.6 Å到7.5 Å的过程中，那对消失的费米双峰又像“复活”一样逐步重新出现，最后与裸Zundel阳离子的光谱完美重合。这意味着，溶剂壳层的存在并不是通过“掩盖”的方式让双峰看不见，而是从根本上破坏了费米共振的“合拍”条件。进一步分析揭示，溶剂壳层把内层水分子的摇摆频率从约100 cm⁻¹硬生生“抬”到了约800 cm⁻¹——就像把一个本来声音低沉的人强行变成了高音歌手，原本和谐的“和声”自然就散了。而有趣的是，内层水分子的另一种运动——非对称弯曲——几乎不受溶剂壳层的影响，其频率在裸核心和eZ配合物中基本保持稳定。这个细节后来证明至关重要：在非对称溶剂环境中，正是这个“倔强不变”的弯曲模式，与质子的转移运动持续耦合，持续贡献着光谱信号。

对称脱溶剂化的分析固然深刻，但它毕竟是一种理想化的“思想实验”，代表不了液态水中质子真实感受到的那种持续不断的“不对称”涨落——毕竟真实的溶液里，质子周围的水分子不会排得那么整整齐齐。为此，研究团队设计了一个更具现实意义的扰动：只把eZ配合物中第二壳层的一个水分子沿着连线方向向外拉，另外三个不动（图4e）。这一下，中心质子感受到的“拉力场”彻底失衡了：原本左右对称的势能面（一个对称的“双阱”结构——想象两个山谷中间隔着一座平缓的山脊，质子可以左右穿梭）迅速变成了一个倾斜的不对称双阱——质子被推着向一侧水分子偏移，逐渐长出Eigen-like的H₃O⁺核心的模样（图4d）。

伴随着这个结构变化，红外光谱也上演了一场“变形记”（图4b）。在1,500到2,300 cm⁻¹区间，一个双峰结构随着拉离距离的增加而不断“漂移”。尤其精彩的是，当拉离距离Rconst约等于3.9 Å时，双峰的两个分量强度变得相等，而它们的中心恰好落在约1,750 cm⁻¹——这与液态水合质子实验光谱中那个著名的蓝色窗口宽谱带的峰值位置几乎完美重合（比较图4b与图2d蓝色窗口）。当非对称性进一步加剧（Rconst > 5 Å）时，双峰逐渐融合并演化为类似裸Eigen阳离子的谱线形状——只不过因为还有一个“拖油瓶”水分子挂在另一侧，所以带有一点点红移（频率向低端偏移）。

与此同时，外层那些氢键也没有闲着。对称环境下，所有外层O-H伸缩振动都挤在约3,200 cm⁻¹的一个单峰里；而随着非对称程度的升高，这些氢键开始“分化”成两个阵营（图4c）。朝向新生H₃O⁺核心的那几个氢键（图4e中绿色质子）逐步向更低频率移动（红移），说明它们变得更“结实”了；而连接核心与那个正在被拉走的水分子之间的氢键（图4e中紫色质子）则逐步向更高频率移动（蓝移），说明这根氢键正在被“扯松”。这种一红一蓝的“跷跷板”效应，恰好对应了实验差谱中2,600至3,500 cm⁻¹那一片巨大频率窗口内的强度变化（比较图4c与图2d绿色及紫色窗口）。换句话说，液态水合质子光谱中这个高频率段的复杂形状，直接反映了局域溶剂环境的非对称程度。

图4 | 扩展Zundel配合物的光谱响应随局域溶剂化环境结构的变化。 a，沿中心质子氢键激发的响应随四个溶剂化水分子对称脱附距离的变化。最上方光谱（橙色）为裸Zundel阳离子H⁺(H₂O)₂的参考谱。b，沿中心质子（e中大的黄色质子）激发的响应随单个溶剂化水分子（e中蓝色框右上方）非对称移除距离的变化。最上方光谱（绿色）为来自参考文献56的裸Eigen阳离子H₃O⁺(H₂O)₃的参考谱。c，外层质子沿与非对称移除后剩余三个外层水分子形成的氢键激发的响应（操作同b）。b和c中的彩色频率窗口标示了与图2d所示实验红外响应的相似性，在正文中讨论。d，基于神经网络CCSD(T)精度势能面（补充信息第2节），沿中心质子坐标（z，如e所定义）的势能剖面，对应于b和c中的非对称单水分子移除；注意对于共线的O₁…H…O₃结构基元，2z等于通常的δ = dO₁-H – dO₃-H坐标。e，蓝色虚线表示通过拉伸Rconst对蓝色框中外层水分子的移除示意。

综合所有这些信息，研究者终于拼出了液态水合质子连续谱的完整“画像”：约1,000–1,400 cm⁻¹的黄色窗口谱带（图2d黄色窗口），来自于中心质子在轻微非对称环境中的伸缩振动（图4b浅黄色谱线）；1,620 cm⁻¹那个相对尖锐的特征峰（图2d中位于蓝色窗口左侧的窄峰），来自第二壳层水分子的纯弯曲运动（图3a中a₀谱带）——它之所以窄而突出，是因为这个频率在Zundel-like和Eigen-like两种构型中几乎不变，相当于一个“稳定锚点”；而1,700–2,300 cm⁻¹的蓝色窗口宽谱带（图2d蓝色窗口），则是中心质子伸缩与内层水非对称弯曲强力“纠缠”的产物（图4b）——两者的耦合在非对称环境中尤其强烈，谱带峰值稳定地“钉”在1,750 cm⁻¹，与实验完美吻合。而约2,300–2,600 cm⁻¹那段相对“安静”的区间，在非对称eZ模型中完全是“无人区”——没有任何振动模式能在那里产生可辨识的信号。研究者推断，这部分光谱很可能来自体相水中那些数量较少但确实存在的Eigen-like构型：如果它们的最短氢键被邻近水分子稍稍拉伸，其伸缩频率就会从约2,700 cm⁻¹红移落入这个区间，从而填上这块“空白”。

这项研究的成功，不仅在于它为水合质子的连续红外光谱提供了迄今为止最为精确的量子动力学解释，更在于它展示了一种全新的研究策略——把高精度的气相团簇量子动力学计算，与液态涨落环境的模拟巧妙地衔接起来。通过系统地引入一个控制参数（拉远一个水分子的距离），研究者建造了一座从“静态极限模型”通往“动态真实图景”的梯度桥梁。展望未来，他们已经将目光投向了更具挑战性的“另一半”——水合氢氧根（OH⁻(aq)）。氢氧根周围的水分子排布比质子更加复杂：在静息态时它表现为超配位结构（配位数大于四），而在活性态时则切换为标准的四配位四面体结构，这个“变形金刚”式的行为将需要更高维度的计算方法和更精密的势能面描述。但一旦攻克，我们将首次获得一幅完整的、关于水相酸碱化学微观基础的统一图景——而这，正是两个多世纪以来无数化学家孜孜以求的目标。