刘建川・GROMACS 实战｜氢键核心原理与氢键计算全解析

大家好，我是刘建川，14 年深耕 GROMACS 分子模拟领域，始终专注分子模拟实战。我将拆解分子模拟核心原理，帮你理清技术脉络、避开选择误区，真正做到学以致用。

一、氢键的基础定义

当氢原子与F、O、N形成共价键时，成键电子会偏向电负性大的F、O、N，使其带部分负电荷，氢原子则呈现近似氢离子的状态，可吸引邻近F、O、N原子的孤对电子。氢原子介于两个电负性大的原子之间形成的特殊相互作用，即为氢键。

判断氢键存在需满足两个条件：一是分子中含有F、O、N这类电负性大的原子；二是原子空间结构允许F、O、N与氢原子相互靠近，并非含这些元素就一定能形成氢键。

二、极性共价键与氢键的关联

原子间通过共用电子对形成的相互作用为共价键，分为极性与非极性。不同原子形成共价键时，因吸引电子能力不同，共用电子对会偏向电负性更强的原子，使其显负电性，另一方显正电性，这类共价键即为极性共价键。

氢键的形成，正是以极性共价键A-H为基础的。

三、氢键的标准形式与典型应用

氢原子与电负性大的原子X以共价键结合，再靠近电负性大、半径小的原子Y（O、F、N），以氢为媒介形成X-H…Y形式的分子间/分子内相互作用，就是氢键。X和Y可为同种原子，也可为不同原子。

典型场景：

• 蛋白质α-螺旋：N-H…O型氢键
• DNA双螺旋：N-H…O、N-H…N型氢键
• 水分子间：O-H…O型氢键

四、氢键的判定规则与本质

1. 存在与电负性极大的原子A形成强极性键的氢原子
2. 存在半径小、电负性大、含孤对电子且带部分负电荷的原子B（F、O、N）
3. 氢键通式：X－H…Y（X、Y为F、O、N）
4. 氢键本质：强极性键(A-H)上的氢核，与电负性大、含孤电子对且带部分负电荷的原子B之间的静电作用力

关于氢键的定义存在两种理解：一种将X－H…Y整体称为氢键，键长为X到Y的距离；另一种将H…Y称为氢键，键长为H到Y的距离，选用数据时需注意区分。氢键键能统一指将X－H…Y分解为HX和HY所需的能量。

五、氢键的饱和性与方向性

氢键区别于范德华力，具备饱和性与方向性：

• 饱和性：氢原子体积小，A、B原子体积大，一个A-H中的氢只能与一个B原子形成氢键，受负离子排斥无法结合第二个
• 方向性：A—H…B共线时相互作用最强，且氢键方向会与B原子孤对电子对称轴一致，保证氢键结构最稳定

六、氢键键能与强度分级

氢键结合能为2–8 kcal/mol，强度介于范德华力与共价键、离子键之间。

• 弱氢键：键能＜25 kJ/mol
• 中等氢键：键能25–40 kJ/mol
• 强氢键：键能＞40 kJ/mol

七、常见氢键平均键能数据

• F—H…:F：155 kJ/mol（40 kcal/mol）
• O—H…:N：29 kJ/mol（6.9 kcal/mol）
• O—H…:O：21 kJ/mol（5.0 kcal/mol）
• N—H…:N：13 kJ/mol（3.1 kcal/mol）
• N—H…:O：8 kJ/mol（1.9 kcal/mol）
• HO—H…:OH3⁺：18 kJ/mol（4.3 kcal/mol）

八、GROMACS氢键判定与分析命令

多数分子力场未包含氢键相关函数，GROMACS通过空间几何判定计算氢键，判定标准为：

• 供体与受体距离r ≤ 0.35 nm
• 氢键角度a ≤ 30°

氢键分析核心命令：

g_hbond -f md.xtc -s md.tpr -num num.xvg -ang ang.xvg -life life.xvg -ac ac.xvg

输出文件说明：

• num.xvg：氢键个数随时间变化曲线
• ang.xvg：氢键角度随时间变化曲线
• life.xvg：氢键寿命随时间变化曲线
• ac.xvg：氢键自相关函数随时间变化曲线

命令运行后，终端会直接输出氢键寿命的关键统计信息。