昊梵体育网

界面水结构怎么证明?原位红外、SFG与分子动力学的综合解读

说明:本文主要介绍固液的空间分布、偶极取向、氢键配位与动态交换信号,以及、振动和频光谱和之间的证据转换。

界面水结构由哪些可观测量定义?

固体接触液态水后,靠近表面的 H2O 会在皮秒尺度持续交换位置和氢键伙伴。这里所说的“结构”是一组统计分布:氧、氢原子沿表面法向的数密度 ρO(z) 和 ρH(z)给出水层位置;水偶极与表面法向的夹角分布描述取向;O···O 距离、O–H···O 构型和供体/受体数描述氢键;驻留时间与重取向相关函数描述交换速率。优化后的单帧构型只保留一个局部极小点,液态界面的热涨落须由轨迹统计表达。

图1. 间隙可控 ATR-IR 装置中样品—水界面区、体相区和倏逝场的空间关系。DOI:10.1021/acs.analchem.5c02765

数密度和取向分布为“有序水”“无序水”给出可计算定义。数密度峰增强表示特定 z 位置的占据概率升高;偶极分布变窄表示取向涨落减小;平均氢键数上升只记录配位数量,氢键寿命、网络连通方式和水—表面键仍可能改变。结构至少包含位置、取向、配位和时间四类量。

ATR 红外倏逝场、SEIRAS 金属近场与 SFG 二阶选择定则具有不同空间权重。ATR 红外的倏逝场可延伸数百纳米,表面增强红外把权重压向金属近场;电偶极近似下,SFG 的二阶响应主要来自反演对称破缺区域。界面选择性与化学选择性是两件不同的事:O–H 共振锁定水的振动环境,光场分布和非线性选择定则决定空间权重。

同一组界面水还受电极晶面、表面终止、电位、离子浓度、pH、同位素组成和温度控制。互证的起点是同一界面状态。实验若在工作电位下采谱,模拟也应设置相近表面电荷或恒电位条件;实验使用浓电解液,模拟盒内的离子/水比例便要保持一致。否则,光谱和轨迹各自成立,却描述了不同的水层。

图2. RuN4、Ru(001) 与 RuSe 表面的界面水 AIMD 快照、氢键供受体分布和相关能量结果。DOI:10.1038/s41467-022-33984-5

原位红外怎样分辨界面水的振动环境?

ATR-SEIRAS 通过金属薄膜近场增强电极附近分子的红外吸收,并以某一参照电位生成差谱。水的 O–H 伸缩带通常位于 2800–3800 cm−1,H–O–H 弯曲带位于约 1600 cm−1。随电位可逆移动的峰位、面积与线形,记录氢键环境、电场 Stark 位移、表面吸附和水分子数量的共同响应。参照电位改变后,正峰和负峰的含义也会改变,谱线方向须和差谱定义一同记录。

图3. RuNC、RuSex–RuNC 和 Ru/C 在酸性、碱性与中性电解液中的 operando ATR-SEIRAS 电位序列。DOI:10.1038/s41467-022-33984-5

在同一材料和电解液内,较低波数常与较强氢键相关,3600 cm−1附近可出现弱氢键或悬空 O–H。宽 O–H 带来自连续构型分布与振动耦合,把 3200、3400 和 3600 cm−1固定命名为三种离散水会压掉电场、离子配位和费米共振的贡献。峰位归属应限定在同一实验口径和对照组内。

HOD/D2O 同位素稀释可削弱 O–H 振动之间的共振耦合,孤立的 O–D 或 O–H 振子更接近局域氢键环境。电位扫描若出现可逆频移,且 HOD 稀释后仍保留相同变化方向,氢键与局域电场的解释会获得更强约束;若峰面积随表面粗糙度同步变化,近场增强因子便是优先排查对象。

图4. 体相水 ATR-IR、PDMS/水界面的间隙可控 ATR-IR 与 SFG 的 O–H 伸缩谱对比。DOI:10.1021/acs.analchem.5c02765

温度跳变或电位阶跃可把静态谱转成时间响应。氢键断裂与重组会改变 O–H 频率涨落,水分子离开近场后又改变吸收权重。谱线恢复时间同时含有结构交换和探测体积更新,需用 HOD 稀释、不同入射偏振或独立扩散时间区分振动弛豫、取向转动与分子迁移。

原位红外能够锁定振动环境随工作条件的演化,但吸收强度还乘入跃迁偶极、近场衰减和取向选择。界面水数量、氢键强度与分子取向均可改变谱带。因而,红外最有分辨力的证据是同一电极上的电位序列、同位素对照、空白表面与重复扫描。一次峰分解得到的面积比例只保留拟合模型内的相对权重。

SFG怎样给出界面选择性/水分子取向

振动和频光谱把红外光与可见光在界面处时空重合,输出频率为两束光之和。电偶极近似下,具有反演对称的体相水二阶响应相互抵消,界面取向分布留下 χ(2)。ssp、sps、ppp 偏振组合选择不同张量分量,可约束 O–H 相对表面法向的取向分布。SFG 强度满足 ISFG ∝ |χ(2)|2IvisIIR,强度变化还包含不同振子之间的相干干涉。

图5. 空气/水界面的 ssp、ppp 复数 SFG 光谱以及均匀局域介电模型与实验谱的比较。DOI:10.1073/pnas.2204156119

常规同相检测得到 |χ(2)|2,相位和正负号在平方后消失。异相检测 HD-SFG 给出复数 χ(2),Imχ(2) 的符号在明确相位约定后可区分净取向方向。正负水分子取向共存时,振幅可因抵消而减小;较弱的 SFG 信号此时对应更小的净极化序,界面水数量仍可维持较高。

带电界面还会引入 χ(3)Φ 项,扩散层中的体相水在静电场作用下产生净取向。有效响应 χeff(2) = χs(2) + χ(3)Φ,离子强度、德拜长度和相位匹配均会改变第二项。将全部 SFG 振幅归给第一层水,会把表面化学响应与长程电场响应混在同一个谱带里。

图6. 空气/水界面的 MD 快照、氢键供受体类型、深度分布和模拟—实验 Imχ(2) 光谱。DOI:10.1073/pnas.2204156119

红外吸收对振动偶极敏感,SFG 另受界面非中心对称和 Raman 活性约束。红外峰位相近而 Imχ(2) 符号相反时,局域氢键强度可以接近,净取向却发生翻转;红外面积改变而 SFG 近乎不变时,水分子数量或近场权重可能变化,取向差保持较小。两类光谱的差异本身就是结构约束。

分子动力学怎样把谱峰还原为结构分布

MD 轨迹可逐帧统计 ρO(z)、ρH(z)、P(cosθ)、氢键供体/受体数、驻留时间和二阶重取向相关函数 C2(t)。ρO(z) 与 ρH(z) 的峰位差约束 H-up、H-down 或近似平行构型;P(cosθ) 的峰宽给出取向涨落;氢键寿命和 C2(t) 反映网络更新速度。固定电荷 MD、恒电位 MD 与 AIMD 对极化和成键的处理不同,模型类型必须进入结果说明。

图7. 本征与含孔石墨烯/水界面的 AIMD 快照、数密度、局域介电函数和电势剖面。DOI:10.1073/pnas.2204156119

轨迹还能前向计算光谱。IR 吸收可由总偶极矩时间相关函数的傅里叶变换得到;SFG 的振动响应常写成极化率与偶极矩时间导数的互相关,χ(2)(ω) ∝ ∫dt eiωt⟨α̇(0)μ̇(t)⟩。前向算谱把结构假设送回实验坐标:模拟须共同复现峰位、线形、Imχ(2) 符号及其电位变化趋势。只匹配一处峰位会留下多组可能构型。

频率映射、经验水模型、交换关联泛函和核量子效应都会改变 O–H 谱。经典刚性水可覆盖较长时间统计,却缺少 O–H 振动自由度;AIMD 保留成键和电子极化,轨迹长度常受限于几十皮秒;机器学习势可延长采样,但训练集须覆盖表面电荷、离子配位和水解离构型。采样长度决定统计误差,势能面决定物理偏差。延长轨迹只能缩小统计置信区间,无法补偿错误的水—表面相互作用。

固定总电荷与恒电位模型对界面电荷涨落采用不同处理。固定总电荷模型在水分子靠近电极时维持同一电子数,恒电位模型允许电极电荷随溶液构型涨落。表面电荷密度—电位标定应与实验的参比电极或零电荷电位相连;离子层化、局域电场和水偶极取向随后才具有相同电位标尺。非对称 slab 还涉及偶极修正、真空区或双电极周期条件。

图8. Au(111)/21 m LiTFSI 界面中 O、H、Li+ 的数密度及水分子法向、偶极和 O–H 取向分布。DOI:10.1038/s41467-022-33129-8

互证可写成逐项约束。原位红外的 O–H 红移若对应 MD 中 O···O 距离缩短、氢键供体数增加,氢键解释获得结构支撑;HD-SFG 的 Imχ(2) 翻号若对应 P(cosθ) 主峰跨过平行取向,偶极翻转得到空间支撑;线宽收窄若伴随频率涨落减小和 C2(t) 衰减变慢,动态网络趋于缓慢。三条证据必须使用相同晶面、电解液、电位区间和同位素组成。任一实验条件改变,都应重新生成匹配的轨迹和模拟谱。

可复查的界面水结论应列出表面原子模型、水层厚度、离子浓度、电位参照、温度和轨迹长度,并给出氢键的 O···O 截断、O–H···O 几何范围、取向矢量定义、差谱参照电位、SFG 相位约定和 χ(3)Φ 处理方式。实验峰位、Imχ(2) 符号、ρO/H(z) 与 P(cosθ) 在同一组条件下闭合,界面水的空间排布与动态交换才获得可追溯的原子尺度约束。