什么是“氢键网络”？从分子定义到图论及质子输运的说明

说明：本文华算科技主要介绍氢键网络的分子定义、图论表示、拓扑统计、动态重连，以及它与水结构、扩散和质子输运之间的联系。

氢键网络怎样由分子和氢键定义？

氢键网络是由氢键相连的分子集合。以水为例，每个 H₂O 分子可视为一个节点；某个 O—H 基团向邻近氧原子的孤对电子区供出氢键时，两节点之间形成一条有方向的边，方向由给体指向受体。多个节点沿边形成链、环和三维簇，研究对象由单次 O—H···O 接触扩展为连接方式与空间分布。

水分子可供出两个氢键，也可接受两个氢键。图论记号 A_mD_n 分别记录接受数和供出数，例如 A₂D₂ 表示四配位局域环境，A₂D₁ 则含一个未参与供氢的 O—H。供受体状态决定边的方向，配位数决定节点的度；链端、分叉点、闭环和贯穿周期盒的连通簇由这些局域状态组合产生。

轨迹后处理以原子坐标和分子身份生成氢键标签。常用几何规则把 O···O 距离小于 3.5 Å、O_A···O_D—H_D 偏离线形构型小于 30°的分子对标记为边；有些研究加入 O···H 距离、成对相互作用能或电子密度条件。距离截断和夹角截断会改变边数、环数及簇尺寸，论文中的网络统计应附带所用定义。

采用周期性复制的模拟盒时，位于盒子相对两侧的水分子仍可按最小镜像距离构成氢键。若一个连通分量沿 x、y 或 z 方向穿过周期盒，它属于贯穿网络；有限水团簇则只有若干彼此分离的连通分量。节点、边、方向和周期性共同构成氢键网络的计算定义，单张结构图中的虚线只是某一时刻的几何投影。

羧酸、酰胺、咪唑和含水多孔晶体也能形成氢键网络。节点可选整分子，也可选具体供受体基团；晶胞重复会把局域二聚体延伸为一维链、二维片或三维骨架。节点粒度必须在拓扑比较中保持一致，原子级节点和分子级节点会产生不同的度、环和最短路径。

网络结构为什么不能用氢键数代替？

平均氢键数描述每个节点连接多少条边，拓扑描述这些边怎样相接。两个体系都可具有平均 3.4 条氢键：一个体系可能富含五元、六元闭环，另一个体系可能含较多长链和分叉。平均配位相近时，环尺寸分布、最大连通簇、最短路径和空间孔隙仍可明显不同，密度涨落与分子迁移也会随之改变。

常用拓扑量包括节点度分布、供受体状态占比、连通分量尺寸、贯穿概率、最短路径长度和环尺寸分布。环统计还依赖算法：最短环、原始环和不可短路环对同一邻接图会给出不同计数。环的定义与最大搜索长度必须固定，温度、压力、组分或外场之间的比较才拥有同一计量口径。

分支结构可用割边、节点介数和最短路径识别。删除一条割边会把连通簇分成两部分；高介数节点频繁出现在多条最短路径中。这类量定位网络中控制远程连接的少数节点和边，平均配位数对它们没有位置分辨能力。

氢键具有多体协同性。一个水分子接受氢键后，局域电子密度和 O—H 极化发生调整，后续供氢能力可能增强或减弱；闭环中的相互作用能也不必等于各分子对能量的简单相加。DFT 的多体能分解、极化响应或高阶簇展开可刻画能量非加和，图拓扑负责标记协同作用发生在哪些连接单元。

液态水的氢键网络通常含五元至七元环以及更长闭环，晶态冰的环分布则受晶型约束。降温可提高六元环占比并压低长环占比，但环数增多本身不自动给出相态。平移扩散、转动关联、密度涨落和自由能盆地分别检验液体运动、取向松弛、长程结构波动及热力学状态。

动态重连怎样进入计算结果？

有限温轨迹把氢键网络写成随时间变化的邻接矩阵 A_ij(t)。热运动使 O···O 距离和分子取向越过截断值，旧边消失，新边生成；节点本身也会平移、转动或改变供受体状态。动态网络关注重连频率、持续时间和协同范围，它讨论的是一组构型的时间统计。

连续时间关联函数在氢键短暂越过截断时立即记为断裂，间歇关联函数允许边在稍后恢复，两者给出的衰减时间不同。邻居交换率记录节点在给定时间窗内更换了多少个氢键伙伴；环涨落 σ(n) 记录 n 元环数的时间波动。边已经重连而各尺寸环的计数仍保持不变时，网络具有动态连接和相对稳定的拓扑组成。

经典分子动力学可覆盖较大的体系和纳秒尺度，网络统计受水模型电荷、极化处理和势函数控制。AIMD 在运动过程中更新电子结构，适合电荷极化、断键和质子交换，但常见轨迹尺寸与时长较小。路径积分分子动力学把氢核离域纳入采样，O—H 伸缩、质子共享和氢键涨落对核量子效应敏感。

静态 DFT 优化给出势能面附近的一个局域极小构型，可比较网络候选结构的总能、晶格和电子密度。有限温网络则依赖轨迹长度、平衡区间和独立初态。环概率与邻居交换率应在多个时间窗内达到稳定统计，短轨迹中偶然保留的水序列不宜写成持续结构。

静电场作用下，液态水的局域配位与重连动力学会共同变化。外加静电场增强第一、第二配位壳层峰时，氧原子的均方位移斜率可明显下降。g_OO(r) 描述空间邻近概率，MSD 描述平移运动；二者来自同一轨迹，却测量不同物理对象。

网络拓扑怎样影响水和质子输运？

体相水中的三维网络容纳不同环尺寸、配位缺陷和局域孔隙，网络重连参与取向松弛、密度涨落与介电响应。表面吸附、离子水合或纳米限域会减少可用邻居并选择分子朝向，三维贯穿网络可能转为二维多边形或一维水序列。连通方向决定哪些分子可经连续氢键序列关联，孔径和表面官能团则控制序列出现的概率。

金属有机纳米管中的限域水可形成交替四聚体和八聚体，簇内 O···O 距离约 2.5—2.7 Å，簇间接触约 3.0—3.4 Å。晶体结构中的红色虚线标出簇内连接，浅蓝虚线标出沿孔轴连接各水簇的弱氢键。分子簇的重复排列把离散局域结构组织成一维网络。

质子结构扩散要求相邻给体和受体形成连续序列，质子转移后还会发生水分子取向调整与氢键重连。网络连通是质子迁移的结构条件，迁移速率还由质子缺陷浓度、自由能势垒、含水量、孔壁电荷和核量子效应控制。连续序列出现频率可由轨迹统计，转移自由能可由约束 AIMD、增强采样或反应型动力学计算获得。

有机分子晶体、氢键有机骨架和聚合物质子导体还包含固定供受体位点。水含量改变可移动节点数量，酸性基团改变质子缺陷浓度，晶界和孔口改变贯穿簇的连续程度。拓扑统计回答连接是否存在及延伸方向，自由能面和电导测量分别给出迁移势垒与宏观响应。

该纳米管在相对湿度由 40% 升至 95% 时，质子电导跨越三个数量级；95% 相对湿度下的 Nyquist 图和 Arrhenius 曲线给出电导响应与 0.22 eV 活化能。结构图确定一维氢键连接，阻抗谱测量电荷输运，NMR 约束水和质子的运动。这三类信号分别检验结构、宏观输运和分子运动。