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.

Calculation of electronic properties of amorphous alloys

Self-consistent supercell approach to alloys with local environment effects

Efficient Calculation of the Lattice Thermal Conductivity by Atomistic Simulations with Ab Initio Accuracy

Simple theory of electronic-structure calculations for amorphous transition-metal alloys

Exafs Measurements and Reverse Monte Carlo Modeling of Atomic Structure in Amorphous Ni80p20 Alloys

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

First-principle and molecular dynamics calculations for physical properties of Ni–Sn alloy system

An Embedded Atom Method for Alloy Nanoparticles

Gaussian-basis LDA and GGA Calculations for Alkali-Metal Equations of State

Computing Elastic Tensors of Amorphous Materials from First-Principles

Prediction of NMR, X-ray and Mössbauer experimental results for amorphous Li-Si alloys using a novel DFTB model

Experimental and molecular dynamics simulation study of structure of liquid and amorphous Ni62Nb38 alloy

First Principles Calculation of Elastic Moduli of Early-Late Transition Metal Alloys

Atomic structure and physical properties of Ni–Nb amorphous alloys determined by an n-body potential

Electronic structure and self energies of randomly substituted solids using density functional theory and model calculations

Structure and Density of Transition Metal-Metalloid Amorphous Alloys

Molecular Dynamic Simulations and Atomic Structures of Amorphous Materials

Development of Lattice Inversion Modified Embedded Atom Method and Its Applications

Multiscale Energy Density Algorithm and Application to Surface Structure of Ni Matrix of Superalloy

First principles residual resistivity using locally self-consistent multiple scattering method

The ab initio amorphous materials database: Empowering machine learning to decode diffusivity