Abstract:The constant potential molecular dynamics simulation method proposed by Siepmann and Sprik and reformulated later by Reed (SR-CPM) has been widely employed to investigate the metallic electrolyte/electrode interfaces, especially for conducting nanochannels with complex connectivity, *e.g.*, carbide-derived carbon or graphene-assembled membrane. This work makes substantial extensions of this seminal SR-CPM approach. First, we introduce two numerical techniques to determine electrode atom charges with an order of magnitude improvement in computational efficiency compared with those widely employed methods. The first numerical technique dramatically accelerates the to calculation of the Ewald interaction matrix $\mathbf{E}$, which takes advantage of the existing highly optimised electrostatic codes. The second technique introduces a new preconditioning technique in the conjugate gradient method to considerably increase the computational efficiency of a linear equation system that determines electrode atomic charges. Our improved SR-CPM implemented in the LAMMPS package can handle extra-large systems, *e.g.*, over 8.1 million electrode atoms. Second, after demonstrating the importance of the electroneutrality constraint, we propose a two-step method to enforce electroneutrality in the following post-treatment step, applicable for matrix and iterative techniques. Third, we propose a solid theoretical analysis for the adjustable parameter $\alpha_i$ (namely the atomic Hubbard-U $U_i^0$), which is arbitrarily selected in many SR-CPM simulation practices. We proposed that the optimised $\alpha_i$ or $U_i^0$ should compensate for the electrical potential/energy discrepancy between the discrete atomistic model and the continuum limit. The analytical and optimal ${\alpha}_i^0$ values are derived for a series of 2D materials.

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

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

Embedded-atom-method interatomic potentials from lattice inversion

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

Analysis of semi-empirical interatomic potentials appropriate for simulation of crystalline and liquid Al and Cu

Long-range empirical potential model: extension to hexagonal close-packed metals

Development of New Interatomic Potentials Appropriate for Crystalline and Liquid Iron

Composition-dependent interatomic potentials: A systematic approach to modelling multicomponent alloys

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

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

Modified Analytic Embedded Atom Method Potential for Chromium

A Response Embedded Atom Method of Interatomic Potentials

Constraints on phase stability, defect energies, and elastic constants of metals described by EAM-type potentials

Meta-Atom Molecular Dynamics for Studying Material Property Dependent Deformation Mechanisms of Alloys

Thermodynamic properties of fcc metals using reparameterized MEAM potentials

Analysis of the EAM and MEAM Potentials for Modeling Localized States of the Ni3Al and Pt3Al Crystals

Quantum corrected cell model for an anharmonic generalized Lennard-Jones solid

Molecular Dynamic Simulation of Liquid-Vapor Coexistence of Metals Modeled Using Modified Empirical Pair Potentials

Development of suitable interatomic potentials for simulation of liquid and amorphous Cu–Zr alloys

Modified embedded-atom method potential for high-temperature crystal-melt properties of Ti–Ni alloys and its application to phase field simulation of solidification

Accurate Path-Integral Molecular Dynamics Calculation of Aluminum with Improved Empirical Ionic Potentials