怎么分析界面水分子取向？

说明：本文华算科技主要介绍界面水取向的几何定义、倾角分布的归一化、不同水组分的取向统计，以及分子模拟与和频光谱之间的联系。

水分子靠近固体、薄膜或气相时，表面配位和氢键会改变不同朝向出现的概率。界面水取向描述分子相对界面的姿态分布，既包含氢端朝向，也包含分子平面的倾斜。单张 DFT 优化构型记录一个局部结构，有限温度下的取向分布通常由 AIMD、经典 MD 或机器学习势轨迹统计得到。

完整 H₂O 中，两根 O→H 键矢量分别描述两个氢的朝向，O 指向两个 H 中点的连线可作为几何偶极方向。选定由界面指向液体的单位法向 n，再用单位分子矢量 ê 定义 cosθ = ê·n。θ 小于 90° 时，该矢量沿法向的分量为正；大于 90° 时为负。

几何偶极方向与电子结构计算得到的分子偶极并非同一个计算对象。前者只使用原子坐标，后者涉及电子密度或局域电荷的分配；强极化、表面电荷转移和分子变形都可能使两者出现偏离。用几何矢量统计姿态时，结果适合描述结构朝向，不应自动转换成真实的电极极化强度。

同一个偶极倾角下，水分子仍可绕偶极轴转动，使两个 O—H 分别朝向不同的表面位点。单根 O—H、几何偶极和分子平面法向对应不同取向变量。粗糙或涨落的界面还存在局域法向；周期单元中两侧界面的法向相反，分别朝向各自液体后，才具有一致的取向正负约定。

三维空间中，各向同性的分子矢量按球面面积均匀分布。θ 以弧度计时，倾角概率密度满足 p_θ(θ) = sinθ/2，因此即使没有定向作用，按 θ 绘制的直方图仍在 90° 附近最高。90° 附近容纳的球面方向多于 0° 或 180° 附近，峰值中含有几何计数的贡献。

令 u = cosθ，归一化的倾角密度与余弦密度满足 p_θ(θ) = p_u(cosθ)sinθ，各向同性参照变为 p_u(u) = 1/2。按度数定义概率密度时，单位换算也须进入归一化。DOPC—水界面的原图以 θ 和距瞬时界面的距离作图，靠近膜与远离膜的亮区不同；解释这些亮区时，应保留原图使用的统计变量和权重。

取向随位置变化时，p(u|z) 描述某个水层内部的朝向概率，联合分布 p(z,u) 还包含水分子出现位置的概率。两个水层可能具有相似的归一化取向分布，却含有不同数量的水。把全部界面区域合并统计，分子数较多的水层会占较大权重；取向变化与水层占据变化须分别计算。

Co₃O₄(001) 表面的氧、氢密度峰可用于定位吸附区域，但原图的氧统计包含水和羟基。以完整水的氧位置选取样本，与把所有含氧物种纳入同一窗口，得到的取向分布会不同。若表面发生水解离，分子身份也要随成键状态更新，不能持续把一个羟基当作有两个氢的 H₂O。

ρ_N(z)⟨u⟩_z 是密度加权的几何取向量，曲线幅度受水分子数量和平均朝向影响；它与层内归一化的 p(u|z) 单位也不同。真实极化密度还包含偶极大小和电子贡献。比较不同电位或覆盖度时，氧密度、取向分布和加权量可分别检验水层占据、朝向重排及其空间贡献。

一阶取向参数 m₁ = ⟨u⟩ 记录朝向正负的不对称。等量分子分别取 u = +1 和 −1 时，m₁ = 0；所有分子取 u = 0 时，m₁ 也为 0。相反朝向相互抵消与全部偶极平躺给出相同均值，却具有完全不同的分布形状。各向同性分布同样满足 m₁ = 0，单个均值没有保留这些差异。

实际氧化物表面还包含化学组分的叠加。Co₃O₄ 上，完整水、源于水的羟基、表面氧质子化形成的羟基具有不同的 O—H 取向分布。将它们合并后，一个峰的增长可能来自该类物种数量增加，也可能来自同类物种重新定向；分组后的分子数和条件分布分别约束这两种变化。

二阶取向参数 S₂ = ⟨(3u²−1)/2⟩可区分部分排列：等量 u = ±1 时 S₂ = 1，全为 u = 0 时 S₂ = −1/2，各向同性时为 0。但等量 u = ±1/√3 也使 m₁ 和 S₂ 同时为 0，分布仍只有两组倾角。有限几个取向参数压缩了信息，完整 p(u) 以及两根 O—H 的联合分布保留了更多姿态差别。

NaOH 水溶液表面的模拟中，Na⁺ 邻近水与 OH⁻ 邻近水可以出现符号相反的净取向。原论文按离子邻近关系分组，并给出随深度变化的取向贡献；其曲线经过密度加权，部分还为便于绘图作了倍数缩放，曲线高度并非直接计数的分子占比。

离子配位水的平均取向改变，还可能来自配位分子与非配位分子之间的比例变化。同一配位类别内的 p(u|z)检验朝向是否重排，该类别的数密度检验组分占据是否变化。对于同一水分子的两根 O—H，合并成一维分布还会丢失两根键的关联；相近的偶极倾角可对应一个氢朝表面、另一个氢沿水层延伸等不同构型。

振动和频光谱 SFG 中，红外光与可见光在界面产生和频信号，其振幅涉及振动响应与分子的取向张量平均。常规强度检测主要得到 |χ⁽²⁾|²；外差检测 HD-VSFG 可提取复响应及其虚部。在统一的相位、法向和振动模式约定下，Imχ⁽²⁾ 的符号约束振动基团的净极性朝向，倾角反演还包含超极化率张量、偏振组合和分布宽度。

NaOH 溶液的模拟可分别计算 Na⁺ 邻近水和与 OH⁻ 配位的 O—H 振动谱，再与实验总谱比较。前向光谱把构型和振动响应转成可比较的谱线：某一水组分的数量、朝向及氢键环境都会改变其贡献。原图中的分项谱来自模拟归属，实验总谱没有逐个识别每个水分子的倾角。

O—H 伸缩频率受氢键和局域电场影响，光谱振幅还包含取向平均与振动耦合。因而，峰位变化与净取向符号是不同观测信息；同一频率范围内，朝向相反的响应还可能抵消。把某个频段的面积直接换算成固定倾角水分子的数量，会引入对单分子响应及分布形状的额外假设。

带电固液界面的有效和频响应还可包含静电场诱导的 χ⁽³⁾ 项，其幅度和相位受电场沿深度的分布与光学相位积累影响。盐浓度改变屏蔽长度后，更远处受场水也可参与有效响应。当屏蔽长度与光学相干长度可比时，不同深度的三阶贡献发生相位混合，测得的谱线可同时包含吸收与色散成分。

原论文的单振子理论算例固定 2,100 cm⁻¹ 的模型频率及振子相位参数，比较不同界面电势下的和频响应。吸收与色散项的混合改变峰形，计算得到的表观峰位随界面电势移动。