Treecode algorithm for pairwise electrostatic interactions with solvent-solute polarization

An O(N log N) treecode algorithm is presented for computing pairwise interactions of electrostatic free energy for reaction potentials with polarization effects due to the macroscopic solvent. A multipole expansion for a cluster is used to account for particles inside the cluster, where a spatial difference is applied to obtain the expansion coefficients of the polarization function. Numerical tests are performed to illustrate the accuracy and efficiency of the approach. The algorithm is significant in speeding up generalized Born methods for biomolecular simulations under the framework of macroscopically treating solvents.