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.

Nonequilibrium molecular dynamics for accelerated computation of ion–ion correlated conductivity beyond Nernst–Einstein limitation

Capturing the Nuclear Quantum Effects in Molecular Dynamics for Lattice Thermal Conductivity Calculations: Using Ice As Example

Size Dependent Lithium-Ion Conductivity of Solid Electrolytes in Machine Learning Molcular Dynamics Simulations

Equivalence of the Equilibrium and the Nonequilibrium Molecular Dynamics Methods for Thermal Conductivity Calculations: from Bulk to Nanowire Silicon

Thermal conductivity of solid electrolytes with ab initio accuracy

Correlations from ion-pairing and the Nernst-Einstein equation

High-throughput investigation of stability and Li diffusion of doped solid electrolytes via neural network potential without configurational knowledge

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

Simulating Diffusion Properties of Solid‐State Electrolytes via a Neural Network Potential: Performance and Training Scheme

Analysis of diffusion in solid state electrolytes through MD-simulations, improvement of the Li-ion conductivity in β-Li3PS4 as an example

Simulating diffusion properties of solid-state electrolytes via a neural network potential: Performance and training scheme

Thermal conductivity of Li$_3$PS$_4$ solid electrolytes with ab initio accuracy

Scalable molecular simulation of electrolyte solutions with quantum chemical accuracy

Accurate description of ion migration in solid-state ion conductors from machine-learning molecular dynamics

Transformer Enables Ion Transport Behavior Evolution and Conductivity Regulation for Solid Electrolyte

Defect-driven anomalous transport in fast-ion conducting solid electrolytes

Atomic cluster expansion force field based thermal property material design with density functional theory level accuracy in non-equilibrium molecular dynamics calculations over sub-million atoms

Correlation between mechanical properties and ionic conductivity of sodium superionic conductors: a relative density-dependent relationship

Which Interfaces Matter Most? A Non-monotonic Relationship between Grain Boundary Defect Chemistry and Ionic Conductivity in a Concentrated Solid Electrolyte

Lattice Dynamics in the NASICON NaZr2(PO4)3 Solid Electrolyte from Temperature-Dependent Neutron Diffraction, NMR, and Ab Initio Computational Studies