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.

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

A Response Embedded Atom Method of Interatomic Potentials

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

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

MetalWalls: Simulating electrochemical interfaces between polarizable electrolytes and metallic electrodes

Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

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

Embedded-atom-method interatomic potentials from lattice inversion

Electrostatic interactions in atomistic and machine-learned potentials for polar materials

A Polarizable Valence Electron Density Based Force Field for High-Energy Interactions between Atoms and Molecules

Immersed boundary method for dynamic simulation of polarizable colloids of arbitrary shape in explicit ion electrolytes

Advanced Potential Energy Surfaces for Molecular Simulation

A fast, convenient, polarizable electrostatic model for molecular dynamics

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

Symmetrized Drude Oscillator Force Fields Improve Numerical Performance of Polarizable Molecular Dynamics

Reconciling simulated melting and ground-state properties of metals with a modified embedded-atom method potential

Development of Nonbonded Models for Metal Cations Using the Electronic Continuum Correction

Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields

The electrostatic interaction of an external charged system with a metal surface: a simplified density functional theory approach

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

Predicting the Charge Density Response in Metal Electrodes