极化率和电荷模型的推导

极化力场最常使用分布原子极化。已经提出了几种方法来推导这些极化率。所有的方法都是基于拟合原子极化来重现实验或量子力学分子极化或量子力学静电势。根据诱导偶极子之间允许的相互作用水平，这些方法也可以分为两组:加性模型和相互作用模型。在加性方法中，允许极化位点响应外部电场，但不响应分子内其他位点上的永久和诱导多极。在米勒的方法中，原子的极化率只取决于杂化态和原子上的电子数。

原子极化是通过用点电荷或外场探测分子并计算响应静电势而产生的。响应势能是在分子周围的网格上定义的带电探针存在和不存在(真空)时的量子力学势能之差。模型参数(如原子极化率αi和高斯指数)拟合为静电势的响应。在计算过程中，允许高斯诱导偶极子和探针电荷相互作用，但在拟合过程中只使用响应势能的贡献，假设分子内极化是恒定的。这样，计算的极化率不依赖于永久多极模型的选择或分子内极化。与许多其他方法不同，永久静电模型可以在推导极化后选择。

感应偶极子、永久电荷和更高的多极矩决定了分子的总静电特性。由于这些原因，永久和分布的原子多极矩需要基于对静电势的拟合来推导，并且需要根据假设的极化模型进行一致的参数化。

详见文献，内容太多看不下去了，总体一句话，所有的拟合都是围绕静电势进行开展的工作，应用开发的方法很多，很多力场有自己的开发算法，所以不再进行描述。