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.

A note on accurate pressure calculations of Coulomb systems with periodic boundary conditions

Pressure of Coulomb systems with volume-dependent long-range potentials

Pressure Calculation in Polar and Charged Systems using Ewald Summation: Results for the Extended Simple Point Charge Model of Water

Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory

Force and Heat Current Formulas for Many-Body Potentials in Molecular Dynamics Simulations with Applications to Thermal Conductivity Calculations

Lekner summations and Ewald summations for quasi-two dimensional systems

Electrostatics of solvated systems in periodic boundary conditions

A Note on Pore Fluid Pressure Ratios in the Coulomb Wedge Theory

From the Vlasov-Poisson-Boltzmann System to the Vlasov-Poisson-Landau System with Coulomb Potential

Methods of Pressure for Fluid Systems with Periodic Boundary Conditions in Molecular Dynamics

Rigorous Error Bounds for Ewald Summation of Electrostatics at Planar Interfaces

Computation of the High Temperature Coulomb Density Matrix in Periodic Boundary Conditions

Self-consistent pressure-dependent on-site Coulomb correction for zero-temperature equations of state of $f$-electron metals

The HI Content and Extent of Low Surface Brightness Galaxies - Could LSB Galaxies be Responsible for Damped Ly-alpha Absorption?

Accurate boundary-integral formulations for the calculation of electrostatic forces with an implicit-solvent model

Some remarks on Coulombic effects in and scattering and the determination of

Symmetries of the Kac-Peterson Modular Matrices of Affine Algebras

Non-Ewald methods for evaluating the electrostatic interactions of charge systems: similarity and difference

Finite-Difference Approximations and Local Algorithm for the Poisson and Poisson-Boltzmann Electrostatics

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