计算范德华力中的诱导力

以二聚体A－B为例。

首先应该计算二聚体的诱导偶极矩，以及固有偶极矩，方法参考如何计算诱导偶极矩。注意单位为Debye。计算过程可以考虑使用Excel表格的计算功能。

A的固有偶极(P_ax,P_ay,P_az)，B的固有偶极矩为(P_bx,P_by,P_bz)；A的诱导偶极矩为(P'_ax,P'_ay,P'_az)，B的诱导偶极矩为(P'_bx,P'_by,P'_bz)；

A,B电偶极子的正电中心、负电中心坐标（A⁺_x,A⁺_y,A⁺_z），（A^-_x,A^-_y,A^-_z），（B⁺_x,B⁺_y,B⁺_z），（B^-_x,B^-_y,B^-_z）计算方法参考分子的电偶极矩、四极矩的计算

先将偶极子还原成正负点电荷的模型，偶极子间相互作用还原为4个点电荷作用。这样做可以使得计算的公式简单，偶极子间的作用公式比较复杂，但实际上是等价的

例如A的固有偶极（Xa,Ya,Za），偶极子正负电荷中心的坐标：（A⁺_x,A⁺_y,A⁺_z）、（A^-_x,A^-_y,A^-_z），二者之间的距离L_A，单位为埃。

统一单位：

偶极矩：1 a.u. = 2.542 Debye = 8.478×10-30 C·m
长度：1 埃＝1.889726 a.u.

则换算出来的正负点电荷为（统一为原子单位a.u.）：

Q=P/L
=(|(Pax,Pay,Paz)|/2.542)/(L<sub>A</sub>＊1.889726)

这样得到＋Q、－Q点电荷对，分别置于坐标（A⁺_x,A⁺_y,A⁺_z）与（A^-_x,A^-_y,A^-_z）

根据诱导力的定义，将(P_ax,P_ay,P_az)与(P'_bx,P'_by,P'_bz)按上面的方式，转换为4个点电荷相互作用（力、能量）； 并将(P_bx,P_by,P_bz)与(P'_ax,P'_ay,P'_az)按上面的方式，转换为4个点电荷相互作用（力、能量）；

二者相加，即得到二聚体之间的诱导力、能量。