Abstract:Atomic description of electrochemical systems requires reactive interaction potential to explicitly describe the chemistry between atoms and molecules and the evolving charge distribution and polarization effects. Calculating Coulomb electrostatic interactions and polarization effects requires a better estimate of the partial charge distribution in molecular systems. However, models such as reactive force fields and charge equilibration (QEq) include Coulomb interactions up to a short-distance cutoff for better computational speeds. Ignoring long-distance electrostatic interaction affects the ability to describe electrochemistry in large systems. We studied the long-range Coulomb effects among charged particles and extended the QEq method to include long-range effects. By this extension, we anticipate a proper account of Coulomb interactions in reactive molecular dynamics simulations. We validate the approach by computing charges on a series of metal-organic frameworks and some simple systems. Results are compared to regular QEq and quantum mechanics calculations. The study shows slightly overestimated charge values in the regular QEq approach. Moreover, our method was combined with Ewald summation to compute forces and evaluate the long-range effects of simple capacitor configurations. There were noticeable differences between the calculated charges with/without long-range Coulomb interactions. The difference, which may have originated from the long-range influence on the capacitor ions, makes the Ewald method a better descriptor of Coulomb electrostatics for charged electrodes. The approach explored in this study enabled the atomic description of electrochemical systems with realistic electrolyte thickness while accounting for the electrostatic effects of charged electrodes throughout the dielectric layer in devices like batteries and emerging solid-state memory.

Local Simulation Algorithms for Coulomb Interaction

Auxiliary field simulation and Coulomb's Law

Monte Carlo simulation of electrostatic interactions in inhomogeneous dielectric media: Correct sampling for the local lattice simulation algorithm

A new approach for efficient simulation of Coulomb interactions in ionic fluids

A multiscale hybrid Maxwellian-Monte-Carlo Coulomb collision algorithm for particle simulations

A RANDOM BATCH EWALD METHOD FOR PARTICLE SYSTEMS WITH COULOMB INTERACTIONS

Local molecular field theory for effective attractions between like charged objects in systems with strong Coulomb interactions

Super-Scalable Molecular Dynamics Algorithm

Learning the Quantum Centroid Force Correction in Molecular Systems: A Localized Approach

Electrostatic interactions in water: nonlocal electrostatic approach

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

Fast electrostatic solvers for kinetic Monte Carlo simulations

Guidelines for Free-Energy Calculations Involving Charge Changes

Surface-resolved Dynamic Simulation of Charged Non-Spherical Particles

A GPU-based Large-Scale Monte Carlo Simulation Method for Systems with Long-Range Interactions.

Local Structure-Preserving Relaxation Method for Charged Systems on Unstructured Meshes

Self-consistent Coulomb interactions for machine learning interatomic potentials

Local Structure-Preserving Relaxation Method for Equilibrium of Charged Systems on Unstructured Meshes

Monte Carlo generation of localised particle trajectories

Grand-canonical molecular dynamics simulations powered by a hybrid 4D nonequilibrium MD/MC method: Implementation in LAMMPS and applications to electrolyte solutions