Abstract:The multilevel summation method (MSM) offers an efficient algorithm utilizing convolution for evaluating long-range forces arising in molecular dynamics simulations. Shifting the balance of computation and communication, MSM provides key advantages over the ubiquitous particle–mesh Ewald (PME) method, offering better scaling on parallel computers and permitting more modeling flexibility, with support for periodic systems as does PME but also for semiperiodic and nonperiodic systems. The version of MSM available in the simulation program NAMD is described, and its performance and accuracy are compared with the PME method. The accuracy feasible for MSM in practical applications reproduces PME results for water property calculations of density, diffusion constant, dielectric constant, surface tension, radial distribution function, and distance-dependent Kirkwood factor, even though the numerical accuracy of PME is higher than that of MSM. Excellent agreement between MSM and PME is found also for interface potentials of air–water and membrane–water interfaces, where long-range Coulombic interactions are crucial. Applications demonstrate also the suitability of MSM for systems with semiperiodic and nonperiodic boundaries. For this purpose, simulations have been performed with periodic boundaries along directions parallel to a membrane surface but not along the surface normal, yielding membrane pore formation induced by an imbalance of charge across the membrane. Using a similar semiperiodic boundary condition, ion conduction through a graphene nanopore driven by an ion gradient has been simulated. Furthermore, proteins have been simulated inside a single spherical water droplet. Finally, parallel scalability results show the ability of MSM to outperform PME when scaling a system of modest size (less than 100 K atoms) to over a thousand processors, demonstrating the suitability of MSM for large-scale parallel simulation.

A Simple and Efficient Lattice Summation Method for Metallic Electrodes in Constant Potential Molecular Dynamics Simulation

ELECTRODE: An electrochemistry package for atomistic simulations

A Modified Elastic Constant Calculation Method for Triclinic-Like Special Quasi-Random Structures: Application to Pd-M (M = Cu, Ag) Solid Solutions

Mass-zero constrained molecular dynamics for electrode charges in simulations of electrochemical systems

Efficient Calculation of the Lattice Thermal Conductivity by Atomistic Simulations with Ab Initio Accuracy

Constant Potential and Constrained Charge Ensembles for Simulations of Conductive Electrodes

Polarizable Potentials For Metals: The Density Readjusting Embedded Atom Method (DR-EAM)

Simulations of Coulomb systems with slab geometry using an efficient 3d Ewald summation method

Multilevel Summation Method for Electrostatic Force Evaluation

An accurate and efficient self-consistent approach for calculating electron transport through molecular electronic devices: including the corrections of electrodes

Numerical Solution of 3D Poisson-Nernst-Planck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment

Incorporating surface polarization effects into large-scale coarse-grained Molecular Dynamics simulation

Scalable molecular simulation of electrolyte solutions with quantum chemical accuracy

Charge equilibration model with shielded long-range Coulomb for reactive molecular dynamics simulations

Bridging the Gap Between Simulated and Experimental Ionic Conductivities in Lithium Superionic Conductors

Prediction of Li intercalation voltages in rechargeable battery cathode materials: effects of exchange-correlation functional, van der Waals interactions, and Hubbard $U$

Multiscale Theory in the Molecular Simulation of Electrolyte Solutions

Marcus-Hush-Chidsey Kinetics at Electrode-Electrolyte Interfaces

Generalized correction to embedded-atom potentials for modeling equilibrium and non-equilibrium properties of metals

Efficient Parallel Implementation of Ewald Summation in Molecular Dynamics Simulations on Multi-Core Platforms

MetalWalls: Simulating electrochemical interfaces between polarizable electrolytes and metallic electrodes