Skip to main content
construction release_alert
Scholars@Duke will be undergoing maintenance April 11-15. Some features may be unavailable during this time.
cancel
Journal cover image

A multipole-based algorithm for efficient calculation of forces and potentials in macroscopic periodic assemblies of particles

Publication ,  Journal Article
Lambert, CG; Darden, TA; Board, JA
Published in: Journal of Computational Physics
January 1, 1996

A new and efficient algorithm based on multipole techniques is presented which calculates the electrostatic forces and potentials in macroscopic periodic assemblies of particles. The fast multipole algorithm (FMA) can be used to compute forces with in the n-particle unit cell in O(n) time. For the cubic lattice, forces due to a 3k X 3k X 3k lattice of images of the unit cell, containing 33k n particles, can be computed in O(nk2 + k3 log k) time to arbitrary precision. The algorithm was readily added onto an existing FMA implementation, and computational results are presented. Accurate electrostatic computations were done on a 38 X 38 X 38 region of 100000-particle unit cells, giving a volume of 28 quadrillion particles at less than a twofold cost over computing the forces and potentials in the unit cell alone. In practice, a k = 4 ⋯ 6 simulation approximates the true infinite lattice Ewald sum forces (including the shape-dependent dipole correction) to high accuracy, taking 25-30 % more time to compute than only the unit cell. The method extends to noncubic unit cell shapes, and noncubic macroscopic shapes. Simple code modifications allowed computation of forces within macroscopic spheres and ellipsoids, and within near-infinite square, circular, and elliptical surfaces formed of unit cubes replicated along two of the three axes. In addition to efficient periodic simulations, the method provides a powerful tool to study limiting behavior of various finite crystal shapes, as well as surface phenomena in molecular dynamics simulations. © 1996 Academic Press, Inc.

Duke Scholars

Published In

Journal of Computational Physics

DOI

ISSN

0021-9991

Publication Date

January 1, 1996

Volume

126

Issue

2

Start / End Page

274 / 285

Related Subject Headings

  • Applied Mathematics
  • 51 Physical sciences
  • 49 Mathematical sciences
  • 40 Engineering
  • 09 Engineering
  • 02 Physical Sciences
  • 01 Mathematical Sciences
 

Citation

APA
Chicago
ICMJE
MLA
NLM
Lambert, C. G., Darden, T. A., & Board, J. A. (1996). A multipole-based algorithm for efficient calculation of forces and potentials in macroscopic periodic assemblies of particles. Journal of Computational Physics, 126(2), 274–285. https://doi.org/10.1006/jcph.1996.0137
Lambert, C. G., T. A. Darden, and J. A. Board. “A multipole-based algorithm for efficient calculation of forces and potentials in macroscopic periodic assemblies of particles.” Journal of Computational Physics 126, no. 2 (January 1, 1996): 274–85. https://doi.org/10.1006/jcph.1996.0137.
Lambert CG, Darden TA, Board JA. A multipole-based algorithm for efficient calculation of forces and potentials in macroscopic periodic assemblies of particles. Journal of Computational Physics. 1996 Jan 1;126(2):274–85.
Lambert, C. G., et al. “A multipole-based algorithm for efficient calculation of forces and potentials in macroscopic periodic assemblies of particles.” Journal of Computational Physics, vol. 126, no. 2, Jan. 1996, pp. 274–85. Scopus, doi:10.1006/jcph.1996.0137.
Lambert CG, Darden TA, Board JA. A multipole-based algorithm for efficient calculation of forces and potentials in macroscopic periodic assemblies of particles. Journal of Computational Physics. 1996 Jan 1;126(2):274–285.
Journal cover image

Published In

Journal of Computational Physics

DOI

ISSN

0021-9991

Publication Date

January 1, 1996

Volume

126

Issue

2

Start / End Page

274 / 285

Related Subject Headings

  • Applied Mathematics
  • 51 Physical sciences
  • 49 Mathematical sciences
  • 40 Engineering
  • 09 Engineering
  • 02 Physical Sciences
  • 01 Mathematical Sciences