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.

Computational study of the effective three-ion interaction potentials in liquid metals with high density of electron gas

Calculation of an effective three-ionic interactions potential in metallic hydrogen

Probing ion-ion and electron-ion correlations in liquid metals within the quantum hypernetted chain approximation

Ground-state properties of metallic solids from ab initio coupled-cluster theory

Nature of Non-Adiabatic Electron-Ion Forces in Liquid Metals

Manyparticle interactions and local structure of the metallic hydrogen at zero pressure

First-principles molecular dynamics of liquid alkali metals based on the quantal hyper-netted chain theory

Electron correlation effects on the dielectric function of liquid metals

Nonperturbative ab initio approach for calculating the electrical conductivity of a liquid metal

Insights into the Electron-Electron Interaction from Quantum Monte Carlo Calculations

Effective potential of electron-electron interaction in the semiinfinite electron gas

Highly accurate electronic structure of metallic solids from coupled-cluster theory with nonperturbative triple excitations

Non-perturbative $ab$ $initio$ approach for calculating the electrical conductivity of a liquid metal

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

An ab initio study of the static, dynamic and electronic properties of some liquid 5d transition metals near melting

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

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

Ab initio evaluation of the electron-ion energy transfer rate in a non-equilibrium warm dense metal

Temperature evolution of the structure of liquid 3d transition metals: MD study

Theoretical investigation of atomic dynamics of molten metals at melting points

Accurate Surface and Finite Temperature Bulk Properties of Lithium Metal at Large Scales using Machine Learning Interaction Potentials