A simple efficient algorithm for molecular simulations of constant potential electrodes

Ranisha S. Sitlapersad,Anthony R. Thornton,Wouter K. den Otter

DOI: https://doi.org/10.1063/5.0171502

IF: 4.304

2024-01-18

The Journal of Chemical Physics

Abstract:Increasingly, society requires high power, high energy storage devices for applications ranging from electric vehicles to buffers on the electric grid. Supercapacitors are a promising contribution to meeting these demands, though there still remain unsolved practical problems. Molecular dynamics simulations can shed light on the relevant molecular level processes in electric double layer capacitors, but these simulations are computationally very demanding. Our focus here is on the algorithmic complexity of the constant potential method (CPM), which uses dedicated electrostatics solvers to maintain a fixed potential difference between two conducting electrodes. We show how any standard electrostatics solver—capable of calculating the energies and forces on all atoms—can be used to implement CPM with a minimum of coding. As an example, we compare our generalized implementation of CPM, based on invocations of the particle–particle–particle–mesh routine of the Large-scale Atomic/Molecular Massively Parallel Simulator, with a traditional implementation based on a dedicated re-implementation of Ewald summation. Both methods yield comparable results on four test systems, with the former achieving a substantial gain in speed and improved scalability. The step from dedicated electrostatic solvers to generic routines is made possible by noting that CPM’s traditional narrow Gaussian point-spread of atomic charges on the electrodes effectively endows point-like atoms with chemical hardness, i.e., an intra-atomic energy quadratic in the charge.

chemistry, physical,physics, atomic, molecular & chemical

What problem does this paper attempt to address?

This paper focuses on how to efficiently handle the problem of constant potential electrodes in molecular dynamics simulations, particularly for electrochemical systems in high-power, high-energy storage devices such as supercapacitors. Traditional methods, such as Constant Potential Method (CPM), use specialized electrostatic solvers to maintain a fixed potential difference between the two electrodes, but this approach is computationally expensive. The paper proposes a generalized Constant Potential Method (GCPM) that allows the use of any standard electrostatic solver, such as Particle-Particle-Particle-Mesh (PPPM) method, to achieve CPM without the need for rewriting dedicated code. This approach not only improves speed, scalability, but also simplifies the code. By treating the electrode atoms as point charges with "chemical hardness," i.e., energy penalty when the atomic charge changes, GCPM avoids the complexity caused by narrow Gaussian charge distribution in traditional CPM. The paper demonstrates the superiority of GCPM over traditional CPM based on the Ewald summation in four test systems, achieving significant speed improvements and better computational complexity. Additionally, the authors discuss the implementation of GCPM in LAMMPS software and provide numerical examples to demonstrate its efficiency. In summary, this paper aims to address the computational efficiency issues in handling constant potential electrodes in molecular simulations and proposes a new algorithm that can simulate key processes in electrochemical systems faster and more flexibly.

A simple efficient algorithm for molecular simulations of constant potential electrodes

ELECTRODE: An electrochemistry package for atomistic simulations

Fully Periodic, Computationally Efficient Constant Potential Molecular Dynamics Simulations of Ionic Liquid Supercapacitors

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

Charging and discharging a supercapacitor in molecular simulations

Reliability of Constant Charge Method for Molecular Dynamics Simulations on EDLCs in Nanometer and Sub-Nanometer Spaces

Fast electrostatic solvers for kinetic Monte Carlo simulations

Mass-zero constrained molecular dynamics for electrode charges in simulations of electrochemical systems

Grand-canonical molecular dynamics simulations powered by a hybrid 4D nonequilibrium MD/MC method: Implementation in LAMMPS and applications to electrolyte solutions

Grand Canonical Ensemble Approaches in CP2K for Modeling Electrochemistry at Constant Electrode Potentials

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Constant-pH Simulations with the Polarizable Atomic Multipole AMOEBA Force Field

Simulating the charging mechanism of a realistic nanoporous carbon-based supercapacitor using a fully polarizable model

Molecular insights on optimizing nanoporous carbon-based supercapacitors with various electrolytes

PMC-IZ: A Simple Algorithm for the Electrostatics Calculation in Slab Geometric Molecular Dynamics Simulations

Advanced Potential Energy Surfaces for Molecular Simulation

Accuracy of the Explicit Energy-Conserving Particle-in-Cell Method for Under-resolved Simulations of Capacitively Coupled Plasma Discharges

Charge equilibration model with shielded long-range Coulomb for reactive molecular dynamics simulations

A Seamless Grid-Based Interface for Mean-Field QM/MM Coupled with Efficient Solvation Free Energy Calculations

Direct Implicit and Explicit Energy-Conserving Particle-in-Cell Methods for Modeling of Capacitively-Coupled Plasma Devices

Scalable molecular simulation of electrolyte solutions with quantum chemical accuracy