用Python分析方阱位能函数与Lennard-Jones 位能函数的近似关系

在我们研究微观世界分子间的相互作用时，分子间的相互作用力是无法避开的话题。然而，分子间作用力十分复杂，大部分作用力原理的研究尚停留在定性描述阶段。不过，我们可以将分子间作用力分为斥力和引力两大类，以方便我们的研究。故人们为了研究分子间的斥力和引力，便构造了经验的位能函数，以便更好地分析分子间的作用力。位能函数模型十分重要，目前已构成了现代统计力学计算和计算机分子模拟的基础。

下面我们讨论两种最为常见的位能函数，即方阱位能函数与Lennard-Jones 位能函数。并分析它们之间的关联，这有助于我们更加深入而形象地理解分子间作用力的效果。

1.方阱位能函数

方阱(square-well) 位能函数是一个最简单的同时考虑分子间排斥和吸引两种相互作用的位能函数。该模型假设分子间吸引力仅在分子间距 r 的一定范围内存在，且为常数。其位能函数如下： 

方阱位能模型将分子看成直径为σ的硬球，而吸引力仅在σ<r≤λσ的范围内存在，且为不随r变化的常数-ε，好像分子周围存在一个方阱，阱宽为σ(λ-1)，阱深为ε。 

 图1 方阱位能函数

方阱(square-well) 位能函数数学处理简单又具有3个可调参数σ、ε和λ，因此在理论计算和分子模拟中仍得到广泛应用。

2.Lennard-Jones 位能函数

最早由数学家琼斯于1924年提出。由于其解析形式简单而被广泛使用，特别是用来描述惰性气体分子间相互作用尤为精确。其位能函数如下：