氢键的本质与量化分析：给体/受体定义、电子结构来源、量子化学证据与能量分解

说明：本文华算科技主要介绍氢键的给体与受体、电子结构来源、量子化学证据、能量分量，以及模型条件对计算结果的影响。

氢键的基本概念是什么？

氢键常写作 X—H···Y，其中 X—H 是氢键给体，Y 是受体原子或受体电子区。IUPAC 推荐定义把氢键描述为带有成键证据的吸引相互作用：氢原子属于较高电负性的 X—H 基团，并与同一分子或另一分子中的电子富集区域相互作用。虚线表示给体与受体之间的接触，并非新增一根普通共价单键。

O—H···O、N—H···O 和 O—H···N 是常见组合，C—H、S—H 等基团在合适的键极化与几何条件下也可供氢；受体既可能是孤对电子，也可能是 π 电子区。元素名称本身不会自动生成氢键标签，给体极化、受体电子密度、相对取向和吸引能量需来自同一个分子构型。

水二聚体给出最小模型：一分子以 O—H 基团供氢，另一分子的氧孤对电子区接受氢键。平衡构型中的 O—H···O 接近线形，体现了明显的取向选择。受体偏离给体键轴时，静电匹配和轨道重叠同步下降；两个分子靠得过近时，电子云重叠引起的交换排斥迅速上升。

给体一端沿 X—H 键延伸方向常呈正静电势，受体孤对电子附近常呈负静电势。静电势互补确定优选取向，受体孤对轨道与 σ^*_X—H 反键轨道之间的离域则改变 X—H 键级、电子占据和振动响应。色散与诱导贡献仍在同一平衡构型中参与总能量。

给体与受体分属两个分子时形成分子间氢键，位于同一分子的两个基团也可形成分子内氢键。后一种作用常固定分子构象，并改变后续聚集时可用的供受体数目。氢键没有统一的能量上下限，中性分子间的弱氢键、离子辅助氢键和短强氢键占据不同距离与能量区间；分类依据来自给受体身份、质子位置和势能面形状。

计算中怎样确认氢键已经形成？

结构优化后可记录 H···Y 距离、X···Y 距离、X—H 键长和 X—H···Y 夹角。距离小于相应范德华半径和、夹角靠近 180°，可把某次接触列为候选氢键。几何阈值属于筛选规则，不同元素组合、温度和轨迹采样会改变合适的截断值，单个距离并未给出吸引能来源。

孤立二聚体常用 E_int=E_AB−E_A−E_B 描述片段靠近后的电子能变化。单体能量若取复合物中的冻结几何，E_int 聚焦片段间作用；单体若分别充分弛豫，数值会混入构型形变。局域基组还会产生基组重叠误差 BSSE，counterpoise 校正可估计该项。

分子轨道或 NBO 分析可追踪受体 n_Y→σ^*_X—H 离域。受体孤对电子向反键轨道贡献后，σ^*_X—H 占据数增加，常伴随 X—H 伸长和伸缩频率降低。水二聚体的复合物轨道跨越两个单体，说明轨道混合确实进入氢键电子结构；轨道分解数值会随基组、局域化方案和分析方法变化。

QTAIM 在电子密度 ρ(r) 中寻找原子间键路径与键临界点，并读取 ρ_BCP、∇²ρ_BCP 和能量密度。NCI/RDG 则用低约化密度梯度区域定位非共价接触，sign(λ₂)ρ 的负值区常与吸引作用相连。键路径是电子密度拓扑特征，氢键归属仍须与给受体身份、几何和相互作用能保持一致。

差分电荷密度可显示 H···Y 附近的极化与电子重排，颜色含义由原图正负等值面约定控制。原子电荷分区给出 Bader、Hirshfeld 或 Mulliken 电荷变化，数值属于各自分区方案。空间重排、原子分区和轨道离域分别回答电子密度在哪里变化、分到哪个原子、哪些轨道参与，三类输出承担不同的计算含义。

IUPAC 列出的成键证据还包括缔合后的能量降低、X—H 伸缩振动变化、核磁屏蔽变化与电子密度特征。晶体衍射给出供受体距离和取向，中子衍射对氢核位置更敏感，红外与拉曼记录键振动，量子化学则给出相互作用能、轨道占据和 ρ(r)。结构、能量、光谱和密度属于四类独立观测量，相符的变化把氢键命名限定在可复核的证据内。

氢键强弱由哪些能量分量共同决定？

SAPT 常把 E_int 分成静电、交换排斥、诱导和色散。静电项来自未扰动片段电荷分布的库仑作用；诱导项记录永久电场造成的极化及相关响应；色散来自片段电子涨落的相关运动；交换排斥随电子云重叠快速增大。平衡氢键能是吸引项与排斥项相加后的净值，任何单项都无法替代 E_int。

水二聚体距离扫描中，O···H 缩短会增强静电和诱导吸引，交换排斥的增幅在近距离区更陡。能量最低点出现在两类趋势相抵的位置。较短 H···Y 距离并未自动给出较大的净稳定化能，因为原子种类、键极化、受体轨道方向和片段形变代价都会改变势能曲线。

同一系列氮杂环—水或氮杂环—过氧化氢复合物中，N—H···O 距离与 BSSE 校正后的稳定化能呈近似线性相关。该相关式适用于结构和计算口径相近的样本。跨越 O—H···O、C—H···π、带电氢键或不同介质后，斜率与截距会随供受体化学环境改变，统一距离阈值会丢失电子结构差别。

多个氢键共享同一分子片段时，前一处极化会改变后一处给受体能力。给体极化增强后续受体作用时出现协同效应；电子竞争或几何挫折占优时出现反协同效应。氢键数量与总稳定化能常呈非线性，二聚体能量的简单相加难以复现环状簇、分子晶体或吸附层中的多体极化。

红移氢键中，n_Y→σ^*_X—H 离域和键极化常使 X—H 伸长，伸缩频率向低波数移动并伴随强度变化。蓝移氢键则可能出现 X—H 缩短与频率升高，局域再杂化、静电收缩和交换作用可超过反键轨道占据带来的伸长趋势。频移方向属于电子结构响应，它与氢键存在、净作用能大小分属不同物理量。

哪些模型条件会改变氢键判断？

交换相关泛函会改变极化、离域和短程排斥，色散修正会改变分子间距离及结合能。局域基组计算应报告基组大小与 BSSE 处理；周期平面波计算则受截断能、k 点、晶胞尺寸和有限浓度影响。泛函与模型误差可达到氢键能量本身的可观比例，数 kcal mol⁻¹ 的能量排序尤其敏感。

外部离子、表面电场和溶剂介电响应会重新极化给体与受体。Li⁺ 位于水二聚体不同位置时，O—H 键长、O···O 距离和 O—H···O 夹角出现不同变化，模拟红外谱中的伸缩峰随之移动。电荷位置决定局域电场方向，同一种离子可对应氢键增强、削弱或构型解离。

不同氢键体系的电荷转移 CT 与伸缩频移 Δν 分布在多组区域，说明二者并非一条普适标定曲线。O—H、N—H 和含芳环样本的轨道能级、极化率与受体类型不同，相同频移可来自不同电子重排幅度。红外峰位适合约束结构变化，定量相互作用能仍由同一计算口径下的总能或能量分解给出。

周期分子晶体和表面吸附模型还包含长程静电、晶格松弛与相邻周期像。截取二聚体适合分析局域轨道和成对作用，周期计算给出的晶格能、相对多晶型能量及晶胞参数包含多体环境。周期条件和晶格松弛会改变局域氢键势能面，团簇大小与晶胞大小应匹配被讨论的结构尺度。

有限温度下，氢键对应距离与夹角分布，而非静态优化得到的唯一构型。AIMD 可统计氢键占有率、寿命和自由能面；路径积分分子动力学把质子零点运动与隧穿涨落纳入核坐标采样。分子晶体、含水界面和质子导体中的同位素效应、宽化振动峰及质子位置分布，分别受温度、采样时长、核量子效应和周期模型尺寸约束。