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.

Lattice Inversion Modified Embedded Atom Method for FCC Metals

Lattice Inversion Modified Embedded Atom Method for Bcc Transition Metals

Development of Lattice Inversion Modified Embedded Atom Method and Its Applications

Embedded-atom-method interatomic potentials from lattice inversion

Development of Interatomic Potentials for FCC Metals Based on Lattice Inversion Method

Thermodynamic properties of fcc metals using reparameterized MEAM potentials

Chen'S Lattice Inversion Embedded-Atom Method for Ni-Al Alloy

Lattice-Inversion Embedded-Atom-Method Interatomic Potentials for Group-VA Transition Metals

A coarse-grained model for fcc metals based on hierarchical coupling between molecular dynamics and isothermal dissipative particle dynamics

MEAM interatomic potential for thermodynamic and mechanical properties of lithium allotropes

New Methods for Acquiring Interatomic Potentials of Complex Material and the Applications to Study of Intermetallic Compounds

Efficient Modeling of Composition-Dependent Lattice Constants of Multicomponent Metallic Solid Solutions Based on High-Throughput Experiment

Interatomic Potentials Between Distinct Atoms from First-Principles Calculation and Lattice-Inversion Method

Atomistic Simulation Study of the FCC and BCC Crystal-Melt Interface Stresses

Modified Möbius inversion transform and interatomic pair potentials in the bcc-metals Mo and Cr

Design and additive manufacturing of a modified face-centered cubic lattice with enhanced energy absorption capability

A bond counting model for accurate prediction of lattice parameter of bcc solid solution alloys

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

A predictive model of impurity diffusion coefficients in face-centered-cubic metallic systems based on machine-learning

Developing a second nearest-neighbor modified embedded atom method interatomic potential for lithium

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