Fast Updating Multipole Coulombic Potential Calculation
Thomas Hoft, Bradley Alpert
We present a numerical method to efficiently and accurately re-compute the Coulomb potential of a large ensemble of charged particles after a subset of the particles undergoes a change of position. Errors are bounded even after a large number of such shifts, making it practical for use in Monte Carlo Markov chain methods in molecular dynamics, computational astrophysics, computational chemistry, and other applications. The method uses truncated multipole expansions of the potential and a tree decomposition of the computational domain to reduce the computational complexity. Computational costs scale logarithmically in the size of the problem. Scaling, accuracy, and efficiency are confirmed with numerical experiments. The new method outperforms a direct calculation for moderate problem sizes.
and Alpert, B.
Fast Updating Multipole Coulombic Potential Calculation, Siam Journal on Scientific Computing, [online], https://doi.org/10.1137/16M1096189, https://tsapps.nist.gov/publication/get_pdf.cfm?pub_id=921841
(Accessed March 1, 2024)