反应场算法

在反应场（reaction field）算法中，作用在一个偶极矩i上的力被分为近程 相互作用和长程相互作用两个部分。近程相互作用部分包括了反应球内（即截断半径rc内）所有偶极矩对偶极矩i的相互作用，长程相互作用部分是反应球外（即截断半径rc以外）的所有偶极矩对偶极矩i的相互作用。其中，反应球外的偶极矩被近似为介电常数ε的连续介质，并在反应球内形成一个反应场。由连 续介质电动力学可以知道，作用在偶极矩i上的反应场的大小正比于反应球的矩。

式中，μj为反应球内的偶极矩，求和遍及反应球内所有偶极矩，也包括偶极矩i。偶极矩i与反应场的作用能的贡献为Ei• μj，被加到反应球内的相互作用上。

反应场算法的缺点是反应球内的偶极子数变化时，能量和受力的变化是不连 续的，但可以通过一个开关函数来解决不连续问题。