HSMA: An O(N) electrostatics package implemented in LAMMPS
arXiv:2104.05260 · doi:10.1016/j.cpc.2022.108332
Abstract
We implement two recently developed fast Coulomb solvers, HSMA3D [J. Chem. Phys. 149 (8) (2018) 084111] and HSMA2D [J. Chem. Phys. 152 (13) (2020) 134109], into a new user package HSMA for molecular dynamics simulation engine LAMMPS. The HSMA package is designed for efficient and accurate modeling of electrostatic interactions in 3D and 2D periodic systems with dielectric effects at the O(N) cost. The implementation is hybrid MPI and OpenMP parallelized and compatible with existing LAMMPS functionalities. The vectorization technique following AVX512 instructions is adopted for acceleration. To establish the validity of our implementation, we have presented extensive comparisons to the widely used particle-particle particle-mesh (PPPM) algorithm in LAMMPS and other dielectric solvers. With the proper choice of algorithm parameters and parallelization setup, the package enables calculations of electrostatic interactions that outperform the standard PPPM in speed for a wide range of particle numbers.
References in corpus (21)
- Steric Effects in Electrolytes: A Modified Poisson-Boltzmann Equation
- Comments on PPPM, FMM, and the Ewald Method for Large Periodic Coulombic Systems
- How to mesh up Ewald sums (II): An accurate error estimate for the P3M algorithm
- Electrostatics in Periodic Slab Geometries I
- Dielectric effects in the self-assembly of binary colloidal aggregates
- Onset of DNA Aggregation in Presence of Monovalent and Multivalent Counterions
- MMM2D: A fast and accurate summation method for electrostatic interactions in 2D slab geometries
- Simulation of Charged Systems in Heterogeneous Dielectric Media via a True Energy Functional
- Efficient and accurate simulation of dynamic dielectric objects
- Effect of Image Forces on Polyelectrolyte Adsorption at a Charged Surface
- Electrolytes between dielectric charged surfaces: Simulations and theory
- A variational formulation of electrostatics in a medium with spatially varying dielectric permittivity
- Simulations of Polyelectrolyte Adsorption to a Dielectric Like-Charged Surface
- Dielectric modulation of ion transport near interfaces
- Efficient solutions of self-consistent mean field equations for dewetting and electrostatics in nonuniform liquids
- A fast spectral method for electrostatics in doubly-periodic slit channels
- Asymmetric electrolytes near structured dielectric interfaces
- A method to compute periodic sums
- An O(N) and parallel approach to integral problems by a kernel-independent fast multipole method: Application to polarization and magnetization of interacting particles
- Accurate and Efficient Numerical Simulation of Dielectrically Anisotropic Particles
- A random batch Ewald method for particle systems with Coulomb interactions