1NH2-4C60 1 4 0 300 100

4NH2-8C60 4 8 0 300 100

4NH2-16C60 4 16 0 300 100

1NH3-2C60 1 2 10 300 100

1NH3-4C60 1 4 10 300 100

4NH3-8C60 4 8 40 300 100

4NH3-16C60 4 16 40 300 100

3。3氢键的力场模型

通常把氢原子成键的原子称为氢键给体（Donor），将与氢形成氢键的原子称为氢键受体（Acceptor）。一些力场通常用Lennard-Jones10-12势能函数来描述氢键连接的氢键受体和供体之间的相互作用的势能。这种势能为

有一些力场中，在势能函数中专门添加一个复杂的函数形式来描述氢键连接的原子之间的相互作用，这个函数一般都会将氢键的几何形状考虑到其中，因此通常会受到氢键受体原子、供体原子和氢原子坐标的影响。例如：YETI力场中描述氢键的形式的氢键势能不仅依赖氢与氢键受体原子之间的距离和氢原子与氢键供体受体原子形成的角度，还与氢键是否在受体原子中所在孤对电子对的方向上发生偏移有关。在GRID程序中则应用了另外一种依赖孤对电子对方向的形式。另外，还有一些力场并未包括精确的氢键势能项，而是依靠静电和范德华作用来重现氢键[8]。