界面水分子取向分析：矢量定义、概率分布、空间条件化与转动相关函数的综合方法

说明：本文华算科技主要介绍界面水分子取向的矢量定义、概率分布、空间条件化、转动相关函数，以及表面状态和光谱信号对取向归属的约束。

界面水分子取向的基本定义是什么？

对轨迹中每个完整 H₂O，可把氧原子指向两个氢原子中点的单位矢量记为 û_μ，把由固体指向液体的表面法向记为 n̂。二者满足 cos θ_μ = û_μ·n̂，θ_μ 描述分子几何平分线相对界面的倾斜。法向若反向，cos θ_μ 的正负会整体互换，H-up 与 H-down 的命名也随之互换。

水有两根 O—H 键，仅记录分子平分线会丢失单根键的朝向。可分别计算 cos θ_OH,1 与 cos θ_OH,2：两根键都靠近平面对应近似平躺构型，一根指向固体而另一根指向液体则属于单氢朝下构型。H₃O⁺、OH⁻ 和解离吸附物具有不同的键数与电荷中心，应按化学物种分别统计。

各向同性水在 cos θ_μ ∈ [−1,1] 上呈均匀分布；以 θ_μ 为横坐标时，积分测度带有 sin θ_μ 权重，90°附近即使没有择优取向也会出现较高概率。P(cos θ_μ) 与 P(θ_μ) 不能共用同一概率基线。比较两套模拟时，横坐标、分箱宽度与归一化方式必须一致。

P₁ = ⟨cos θ_μ⟩ 保留上下方向，等量的 H-up 与 H-down 会互相抵消；P₂ = ⟨(3cos²θ_μ−1)/2⟩ 对正反方向不敏感，可区分沿法向排列与近似平躺。P₁ = 0 既可能对应各向同性，也可能对应两个方向相反的窄峰，说明不同构型可能产生相同均值，原始直方图仍需保留。

分子平分线绕自身旋转后，θ_μ 不变，两根 O—H 键相对界面的姿态却已改变。第二个欧拉变量 ψ、H—H 矢量或水分子平面法向可记录这类双轴取向。空气—水界面的 ρ_ORI(θ,ψ) 分布把平躺、指向液相和指向气相的构型分开，单变量 P(cos θ_μ) 会把其中一部分投影到同一区域。

经典刚性水模型中，几何平分线通常与模型偶极方向重合。AIMD 里的电子密度会随氢键、表面配位和电荷转移变化；若研究对象是真实分子偶极而非核坐标取向，可用最大局域 Wannier 中心或分子电子分区构造偶极矢量，并明确电子分配方案。几何取向与电子极化属于两种输出，不应共用一个符号后混报。

怎样从分子动力学轨迹得到取向分布？

周期轨迹中的 H₂O 应保持分子完整；若 O 与 H 被写在模拟盒两侧，需按最小镜像关系重组同一水分子。每一帧还要用同一表面参考面对齐 z 坐标；双界面模型的左右表面法向相反，右侧界面的 cos θ_μ 需按该侧 n̂ 计算。若把全局 +z 轴机械用于两侧，两个物理相同的界面会出现符号相反的假差异。

联合分布 P(cos θ_μ,z) 同时保存水氧距表面的法向位置与取向。常见做法是在每个 z 分箱内累积 cos θ_μ 直方图，再按该分箱的水分子数归一化；另行计算 ρ_O(z) 可标出第一、第二水合层。Au(100)/NaCl(aq) 的有限场模拟同时给出水密度、ρ_H₂Ocos θ_μ 和界面角分布，能区分阴极、阳极及其电位响应。

⟨cos θ_μ⟩_z 只统计已经位于该 z 区间的水偏向哪一侧；ρ_O(z)⟨cos θ_μ⟩_z 还计入该处水分子数，更接近偶极极化密度。稀疏区的条件均值可能很大，但它对总界面极化的贡献很小，图表应把样本数或 ρ_O(z) 一并列出。

两根 O—H 键的联合直方图 P(cos θ_OH,1,cos θ_OH,2|z) 可区分双氢朝下、单氢朝下、双氢朝上与平躺族群。为消除 H₁/H₂ 标签任意交换造成的镜像，可按 cos θ_OH 的大小排序，或对两个氢做对称化。构型比例应由明确的分区范围计算，分区线与表面法向定义需随数据发布。

Pt(111) 的 AIMD 轨迹显示，低氢覆盖时靠近表面的水可呈较平的化学吸附构型；高氢覆盖把最近的水氧推离表面，⟨cos θ_μ⟩_z 随 z 改变。误差棒在水数较少的分箱中明显增宽，局部峰值必须附带采样误差，否则单个分箱的符号翻转容易被写成稳定取向转换。

置信区间可由独立轨迹或块平均获得：把生产轨迹分成若干长度超过主要相关时间的区块，分别计算分布，再报告区块间标准误。若相邻帧间隔远短于转动弛豫时间，把每帧当独立样本会低估误差。轨迹长度、取样步长、分箱宽度与有效样本数应随 P(cos θ,z) 共同报告。

取向统计为何要绑定水层/氢键/时间？

第一水层与第二水层通常由 ρ_O(z) 相邻峰之间的谷值划分，固定取 3 Å 或 5 Å 宽度会把不同晶面上的配位区混合。O、H 数密度峰的相对位置还能辅助识别氢端偏向固体或液体，但峰位差是总体统计量，分子级分类仍由同一颗水的 O—H 矢量完成。

界面水既可向表面羟基供氢，也可从表面羟基受氢；电解质离子又会占据水的供受体位置。把每帧水分为表面氢键供体、表面氢键受体、离子配位水与未配位水，再计算条件分布 P(cos θ_μ|z,s)，可判断某个取向峰来自电场极化、表面氢键，还是离子水合构型。氢键距离与夹角截断值应固定，并对临界附近的结果做阈值敏感检验。

静态分布记录各构型出现的概率，却不记录一颗水维持该取向多久。驻留相关函数描述水是否停留在指定水层，转动相关函数 C_l(t)=⟨P_l[û(0)·û(t)]⟩ 描述分子方向记忆。C₁(t) 对 180°翻转敏感，C₂(t) 与许多转动弛豫测量的对称性更接近，两者数值不宜混用。

若只在 t=0 选取第一水层分子，后续离开界面的水仍进入 C_l(t)，所得衰减混有交换过程。另一种口径要求分子在 0 到 t 内持续留在该层，得到的相关函数偏向长驻留族群。文章应标明起始选择、末态选择或连续驻留权重中的哪一种，并把取向弛豫与水层交换分成两条曲线。

相关曲线常含快 libration 与慢氢键重排两个时间范围，单指数拟合会把二者压成一个时间常数。可报告积分相关时间 τ_l=∫C_l(t)dt，或使用具有明确时间窗的多指数/拉伸指数拟合。拟合上限不得超过轨迹中已有的稳定衰减区，长时间尾部应附区块误差，弛豫时间才有可比较的数值含义。

哪些模型条件会改变界面水取向？

晶面终止、表面羟基密度、粗糙位点和吸附覆盖度会改变水—表面配位。γ-Al₂O₃ [110] 与 [100] 的 P(cos θ_μ) 已有明显差异；加入 NaCl、乙酸铵、硝酸钡或乙酸钡后，第一、第二水合层的分布响应也依晶面而变。晶面与水层编号必须作为曲线标签，材料名称本身不足以定义界面。

电化学界面还受电极电位、表面电荷和特异吸附控制。固定电荷轨迹比较的是给定 σ 下的取向响应；恒电位或有限场轨迹允许电极电荷随液体构型涨落。报告结果时应写出电位参照、相对零电荷电位的位置、平均表面电荷密度，并列出离子浓度与吸附覆盖度，避免把电场效应和吸附物排水效应归到同一原因。

静态 DFT 优化可比较若干 H₂O 构型的局部极小点和能量排序，室温取向概率还包含构型熵与氢键交换。AIMD 提供电子结构精度下的短轨迹，经典 MD 或机器学习势可扩展到纳秒和更大界面。不同水模型、泛函、色散修正与核量子处理会改变 O—H 振动、氢键强度和转动速率，跨方法比较宜使用相同取向定义及相近热力学状态。

振动和频光谱的 Imχ⁽²⁾ 对界面净取向敏感，谱峰符号与偏振组合可约束 O—H 族群的朝向。谱—结构映射还包含分子超极化率、振动耦合、局域场和带电界面的 χ⁽³⁾Φ 项；θ 单变量模型还会平均掉 ψ 相关的双轴分布。模拟轨迹前向计算 χ⁽²⁾可检验取向集合是否产生相符的偏振与频率响应，单凭 ⟨cos θ_μ⟩ 的正负无法重建完整光谱。

取向数据发布时应写明 n̂ 的方向、û_μ/û_OH 的定义、z 分箱与水层范围、P(cos θ,z) 的归一化、P₁/P₂、氢键或离子条件分类、C₁(t)/C₂(t) 的分子选择规则；每个 z 区间的样本数与块平均误差也应一同列出。