CALCULATION OF THE REACTION FIELD DUE TO OFF-CENTER POINT MULTIPOLES

Two general methods for calculating the reaction field generated by a set of off-center point multipoles in a spherical cavity are presented . The methods are a generalization of Kirkwood's original theory for a n arbitrary charge distribution. A polytensor formulation, similar to that previously developed for direct multipole interactions, serves to organize the computation and allows straightforward extension to high er derivatives of reaction field energy and gradients of the potential . The computation is reduced to calculation of the Cartesian derivativ es of biaxial harmonics. Recursive and explicit formulas are given for the calculation. As an example, the incorporation of reaction field e ffects in computation of induced dipole moments is discussed. The seco nd procedure, the central multipole method, scales linearly in calcula tion time with the size of the system. Methods to obtain derivatives a nalytically based on this method are also described. Our developments allow use of reaction field energy terms with atomic multipole-based e mpirical potential energy functions. Both methods show particular prom ise for use in simulation of heterogeneous systems, such as biopolymer s, where the remainder of the cavity can be filled explicitly with sol vent. (C) 1997 American Institute of Physics.