skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: An O(N) and parallel approach to integral problems by a kernel-independent fast multipole method: Application to polarization and magnetization of interacting particles

Abstract

Large classes of materials systems in physics and engineering are governed by magnetic and electrostatic interactions. Continuum or mesoscale descriptions of such systems can be cast in terms of integral equations, whose direct computational evaluation requires O(N2) operations, where N is the number of unknowns. Such a scaling, which arises from the many-body nature of the relevant Green's function, has precluded wide-spread adoption of integral methods for solution of large-scale scientific and engineering problems. In this work, a parallel computational approach is presented that relies on using scalable open source libraries and utilizes a kernel-independent Fast Multipole Method (FMM) to evaluate the integrals in O(N) operations, with O(N) memory cost, thereby substantially improving the scalability and efficiency of computational integral methods. We demonstrate the accuracy, efficiency, and scalability of our approach in the context of two examples. In the first, we solve a boundary value problem for a ferroelectric/ferromagnetic volume in free space. In the second, we solve an electrostatic problem involving polarizable dielectric bodies in an unbounded dielectric medium. Lastly, the results from these test cases show that our proposed parallel approach, which is built on a kernel-independent FMM, can enable highly efficient and accurate simulations and allow formore » considerable flexibility in a broad range of applications.« less

Authors:
 [1];  [2]; ORCiD logo [3];  [4];  [5];  [6];  [7]; ORCiD logo [8]
  1. Argonne National Lab. (ANL),Lemont, IL (United States)
  2. Univ. of Chicago, Chicago, IL (United States)
  3. Argonne National Lab. (ANL), Lemont, IL (United States)
  4. Stanford Univ., Stanford, CA (United States)
  5. Argonne National Lab. (ANL), Lemont, IL (United States); Univ. of Chicago, Chicago, IL (United States); KCG Holdings, Inc. (United States)
  6. Univ. of Chicago, Chicago, IL (United States); Univ. Nacional de Colombia-Medellin, Medellin (Columbia)
  7. Argonne National Lab. (ANL),Lemont, IL (United States); Univ. of Chicago, Chicago, IL (United States)
  8. Argonne National Lab. (ANL),Lemont, IL (United States); Northwestern-Argonne Institute for Science and Engineering, Evanston, IL (United States)
Publication Date:
Research Org.:
Argonne National Lab. (ANL), Argonne, IL (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC-22). Materials Sciences and Engineering Division; USDOE
OSTI Identifier:
1339573
Alternate Identifier(s):
OSTI ID: 1420630
Grant/Contract Number:  
AC02-06CH11357
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Chemical Physics
Additional Journal Information:
Journal Volume: 145; Journal Issue: 6; Journal ID: ISSN 0021-9606
Publisher:
American Institute of Physics (AIP)
Country of Publication:
United States
Language:
English
Subject:
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Citation Formats

Jiang, Xikai, Li, Jiyuan, Zhao, Xujun, Qin, Jian, Karpeev, Dmitry, Hernandez-Ortiz, Juan, de Pablo, Juan J., and Heinonen, Olle. An O(N) and parallel approach to integral problems by a kernel-independent fast multipole method: Application to polarization and magnetization of interacting particles. United States: N. p., 2016. Web. doi:10.1063/1.4960436.
Jiang, Xikai, Li, Jiyuan, Zhao, Xujun, Qin, Jian, Karpeev, Dmitry, Hernandez-Ortiz, Juan, de Pablo, Juan J., & Heinonen, Olle. An O(N) and parallel approach to integral problems by a kernel-independent fast multipole method: Application to polarization and magnetization of interacting particles. United States. doi:10.1063/1.4960436.
Jiang, Xikai, Li, Jiyuan, Zhao, Xujun, Qin, Jian, Karpeev, Dmitry, Hernandez-Ortiz, Juan, de Pablo, Juan J., and Heinonen, Olle. Wed . "An O(N) and parallel approach to integral problems by a kernel-independent fast multipole method: Application to polarization and magnetization of interacting particles". United States. doi:10.1063/1.4960436. https://www.osti.gov/servlets/purl/1339573.
@article{osti_1339573,
title = {An O(N) and parallel approach to integral problems by a kernel-independent fast multipole method: Application to polarization and magnetization of interacting particles},
author = {Jiang, Xikai and Li, Jiyuan and Zhao, Xujun and Qin, Jian and Karpeev, Dmitry and Hernandez-Ortiz, Juan and de Pablo, Juan J. and Heinonen, Olle},
abstractNote = {Large classes of materials systems in physics and engineering are governed by magnetic and electrostatic interactions. Continuum or mesoscale descriptions of such systems can be cast in terms of integral equations, whose direct computational evaluation requires O(N2) operations, where N is the number of unknowns. Such a scaling, which arises from the many-body nature of the relevant Green's function, has precluded wide-spread adoption of integral methods for solution of large-scale scientific and engineering problems. In this work, a parallel computational approach is presented that relies on using scalable open source libraries and utilizes a kernel-independent Fast Multipole Method (FMM) to evaluate the integrals in O(N) operations, with O(N) memory cost, thereby substantially improving the scalability and efficiency of computational integral methods. We demonstrate the accuracy, efficiency, and scalability of our approach in the context of two examples. In the first, we solve a boundary value problem for a ferroelectric/ferromagnetic volume in free space. In the second, we solve an electrostatic problem involving polarizable dielectric bodies in an unbounded dielectric medium. Lastly, the results from these test cases show that our proposed parallel approach, which is built on a kernel-independent FMM, can enable highly efficient and accurate simulations and allow for considerable flexibility in a broad range of applications.},
doi = {10.1063/1.4960436},
journal = {Journal of Chemical Physics},
number = 6,
volume = 145,
place = {United States},
year = {2016},
month = {8}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 2 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs
journal, January 1998


Biomolecular electrostatics using a fast multipole BEM on up to 512 gpus and a billion unknowns
journal, June 2011

  • Yokota, Rio; Bardhan, Jaydeep P.; Knepley, Matthew G.
  • Computer Physics Communications, Vol. 182, Issue 6
  • DOI: 10.1016/j.cpc.2011.02.013

Fast boundary methods for magnetostatic interactions in micromagnetics
conference, January 2003

  • Forster, H.; Schrefl, T.; Fidler, J.
  • Digests of INTERMAG 2003. International Magnetics Conference, Digest of INTERMAG 2003. International Magnetics Conference (Cat. No.03CH37401)
  • DOI: 10.1109/INTMAG.2003.1230308

An iterative, fast, linear-scaling method for computing induced charges on arbitrary dielectric boundaries
journal, April 2010

  • Tyagi, Sandeep; Süzen, Mehmet; Sega, Marcello
  • The Journal of Chemical Physics, Vol. 132, Issue 15
  • DOI: 10.1063/1.3376011

The black-box fast multipole method
journal, December 2009


Comparison of scalable fast methods for long-range interactions
journal, December 2013


Compression of boundary element matrix in micromagnetic simulations
journal, April 2009

  • Knittel, A.; Franchin, M.; Bordignon, G.
  • Journal of Applied Physics, Vol. 105, Issue 7
  • DOI: 10.1063/1.3072032

Particle mesh Ewald: An N ⋅log( N ) method for Ewald sums in large systems
journal, June 1993

  • Darden, Tom; York, Darrin; Pedersen, Lee
  • The Journal of Chemical Physics, Vol. 98, Issue 12
  • DOI: 10.1063/1.464397

Hybrid method for computing demagnetizing fields
journal, March 1990

  • Fredkin, D. R.; Koehler, T. R.
  • IEEE Transactions on Magnetics, Vol. 26, Issue 2
  • DOI: 10.1109/20.106342

Comparison of efficient techniques for the simulation of dielectric objects in electrolytes
journal, June 2015


Boundary Integral Methods for Multicomponent Fluids and Multiphase Materials
journal, May 2001

  • Hou, T. Y.; Lowengrub, J. S.; Shelley, M. J.
  • Journal of Computational Physics, Vol. 169, Issue 2
  • DOI: 10.1006/jcph.2000.6626

A theory of interactions between polarizable dielectric spheres
journal, May 2016


A regularized boundary integral equation method for elastodynamic crack problems
journal, January 1989

  • Nishimura, N.; Kobayashi, S.
  • Computational Mechanics, Vol. 4, Issue 4
  • DOI: 10.1007/BF00301390

The fast multipole boundary element method for potential problems: A tutorial
journal, May 2006


Read and write processes, and head technology for perpendicular recording
journal, March 2009

  • Gao, Kai Zhong; Heinonen, Olle; Chen, Yonghua
  • Journal of Magnetism and Magnetic Materials, Vol. 321, Issue 6
  • DOI: 10.1016/j.jmmm.2008.05.025

Direct observation of particle interactions and clustering in charged granular streams
journal, July 2015

  • Lee, Victor; Waitukaitis, Scott R.; Miskin, Marc Z.
  • Nature Physics, Vol. 11, Issue 9
  • DOI: 10.1038/nphys3396

Task-Based FMM for Multicore Architectures
journal, January 2014

  • Agullo, Emmanuel; Bramas, Bérenger; Coulaud, Olivier
  • SIAM Journal on Scientific Computing, Vol. 36, Issue 1
  • DOI: 10.1137/130915662

Geometrical validity of curvilinear finite elements
journal, January 2013


libMesh : a C++ library for parallel adaptive mesh refinement/coarsening simulations
journal, November 2006

  • Kirk, Benjamin S.; Peterson, John W.; Stogner, Roy H.
  • Engineering with Computers, Vol. 22, Issue 3-4
  • DOI: 10.1007/s00366-006-0049-3

A fast solution method for three-dimensional many-particle problems of linear elasticity
journal, August 1998


Computing induced charges in inhomogeneous dielectric media: Application in a Monte Carlo simulation of complex ionic systems
journal, April 2004


Comparison of three-dimensional Poisson solution methods for particle-based simulation and inhomogeneous dielectrics
journal, July 2012


Magnetization reversal in circularly exchange-biased ferromagnetic disks
journal, January 2009


Phase separation versus supersolid behavior in frustrated antiferromagnets
journal, May 2011

  • Albuquerque, A. Fabricio; Laflorencie, Nicolas; Picon, Jean-David
  • Physical Review B, Vol. 83, Issue 17
  • DOI: 10.1103/PhysRevB.83.174421

GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
journal, July 1986

  • Saad, Youcef; Schultz, Martin H.
  • SIAM Journal on Scientific and Statistical Computing, Vol. 7, Issue 3
  • DOI: 10.1137/0907058

Efficient and accurate simulation of dynamic dielectric objects
journal, February 2014

  • Barros, Kipton; Sinkovits, Daniel; Luijten, Erik
  • The Journal of Chemical Physics, Vol. 140, Issue 6
  • DOI: 10.1063/1.4863451

Electrostatic potential inside ionic solutions confined by dielectrics: a variational approach
journal, January 2001

  • Allen, Rosalind; Hansen, Jean-Pierre; Melchionna, Simone
  • Physical Chemistry Chemical Physics, Vol. 3, Issue 19
  • DOI: 10.1039/b105176h

Fast directional multilevel summation for oscillatory kernels based on Chebyshev interpolation
journal, February 2012

  • Messner, Matthias; Schanz, Martin; Darve, Eric
  • Journal of Computational Physics, Vol. 231, Issue 4
  • DOI: 10.1016/j.jcp.2011.09.027

The fast multipole boundary element method for molecular electrostatics: An optimal approach for large systems
journal, July 1995

  • Bharadwaj, Ranganathan; Windemuth, Andreas; Sridharan, S.
  • Journal of Computational Chemistry, Vol. 16, Issue 7
  • DOI: 10.1002/jcc.540160707

Order N algorithm for computation of electrostatic interactions in biomolecular systems
journal, December 2006

  • Lu, B.; Cheng, X.; Huang, J.
  • Proceedings of the National Academy of Sciences, Vol. 103, Issue 51
  • DOI: 10.1073/pnas.0605166103

A 10000 particle molecular dynamics model with long range forces
journal, September 1973


Spectral Approximation of the Free-Space Heat Kernel
journal, July 2000

  • Greengard, Leslie; Lin, Patrick
  • Applied and Computational Harmonic Analysis, Vol. 9, Issue 1
  • DOI: 10.1006/acha.2000.0310

Numerical methods for the stray-field calculation: A comparison of recently developed algorithms
journal, January 2013


Direct Observation of Unconventional Topological Spin Structure in Coupled Magnetic Discs
journal, February 2012


Phase diagram of complex fluids using an efficient integral equation method
journal, November 2005

  • Charpentier, I.; Jakse, N.
  • The Journal of Chemical Physics, Vol. 123, Issue 20
  • DOI: 10.1063/1.2117010

Dielectric Effects in the Self-Assembly of Binary Colloidal Aggregates
journal, July 2014


The finite element method with penalty
journal, May 1973


Demagnetizing factors for rectangular ferromagnetic prisms
journal, March 1998

  • Aharoni, Amikam
  • Journal of Applied Physics, Vol. 83, Issue 6
  • DOI: 10.1063/1.367113