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.

Rapidly Convergent Expansion Method For Calculating The Effective Conductivity Of Three-Dimensional Lattices Of Symmetric Inclusions

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

A fast numerical method for the conductivity of heterogeneous media with Dirichlet boundary conditions based on discrete sine–cosine transforms

Research on the Calculation Method of Convection Effects on Three-Dimensional Dendritic Growth Based on LBM

On the computation of lattice sums without translational invariance

New Iterative Method for Three-Dimensional Eddy-Current Problems

Effective conductivity of a random suspension of highly conducting spherical particles

Rapidly convergent representations for Green' functions for Laplace' equation

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

The Simulation Of Conductivity Of Binary Mixtures With Lattices Boltzmann Method

Series expansion for the effective conductivity of a periodic dilute composite with thermal resistance at the two-phase interface

The electrostatic potential of a periodic lattice

A robust lattice Boltzmann scheme for high-throughput predicting effective thermal conductivity of reinforced composites

Effective Conductivities of Heterogeneous Media Containing Multiple Inclusions with Various Spatial Distributions

A Discrete Sine‐Cosine Transforms Galerkin Method for the Conductivity of Heterogeneous Materials With Mixed Dirichlet/Neumann Boundary Conditions

3D modelling and sensitivity in DC resistivity using charge density

Estimation of Conductivity of Binary Mixtures Using the Lattice Boltzmann Method

Calculation of Effective Characteristics of a 2D Composite with Rhombic Voids Using an Inhomogeneous Cell Model

Effective Conductivity of Spiral and other Radial Symmetric Assemblages

A generalized differential scheme for the effective conductivity of percolating microinhomogeneous materials with the Hall effect

Convergence of direct space exchange lattice sums in polymer band structure calculations