Optimal particle-mesh algorithms

Abstract Optimal particle-mesh algorithms are regarded as those which, for a given computational cost, best represent the physics relevant to the evolution of a computer experiment. A method to determine for arbitrary interparticle force laws a best combination of charge assignment and potential solvers in both momentum and energy conserving schemes is presented. Explicit expressions for the errors in forces and harmonic amplitudes together with expressions for influence functions which minimize those errors are given. A comparison of optimized versions of the common energy and momentum conserving schemes is in Section III, and it is shown how the results of the comparison may be used in deciding which scheme to use for some particular purpose. The application of the error minimizing method to P3M schemes for ionic systems and to collisionless plasmas is discussed in Sections IV and V, respectively. In Section VI, it is shown how the method may be used to obtain finite difference equations for collisionless plasma simulations.