Accelerating QM/MM simulations of electrochemical interfaces through machine learning of electronic charge densities

Abstract

A crucial aspect in the simulation of electrochemical interfaces consists in treating the distribution of electronic charge of electrode materials that are put in contact with an electrolyte solution. Recently, it has been shown how a machine-learning method that specifically targets the electronic charge density, also known as SALTED, can be used to predict the long-range response of metal electrodes in model electrochemical cells. In this work, we provide a full integration of SALTED with MetalWalls, a program for performing classical simulations of electrochemical systems. We do so by deriving a spherical harmonics extension of the Ewald summation method, which allows us to efficiently compute the electric field originated by the predicted electrode charge distribution. We show how to use this method to drive the molecular dynamics of an aqueous electrolyte solution under the quantum electric field of a gold electrode, which is matched to the accuracy of density-functional theory. Notably, we find that the resulting atomic forces present a small error of the order of 1 meV/Å, demonstrating the great effectiveness of adopting an electron-density path in predicting the electrostatics of the system. Upon running the data-driven dynamics over about 3 ns, we observe qualitative differences in the interfacial distribution of the electrolyte with respect to the results of a classical simulation. By greatly accelerating quantum-mechanics/molecular-mechanics approaches applied to electrochemical systems, our method opens the door to nanoseconds timescales in the accurate atomistic description of the electrical double layer.

10 pages, 6 figures

References in corpus (11)

• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• E(3)-Equivariant Graph Neural Networks for Data-Efficient and Accurate Interatomic Potentials
• A Fourth-Generation High-Dimensional Neural Network Potential with Accurate Electrostatics Including Non-local Charge Transfer
• Incorporating long-range physics in atomic-scale machine learning
• Charge fluctuations from molecular simulations in the constant-potential ensemble
• Simulating electrochemical systems by combining the finite field method with a constant potential electrode
• A semiclassical Thomas-Fermi model to tune the metallicity of electrodes in molecular simulations
• Effect of the metallicity on the capacitance of gold - aqueous sodium chloride interfaces
• Predicting the Charge Density Response in Metal Electrodes
• PiNNwall: Heterogeneous Electrode Models from Integrating Machine Learning and Atomistic Simulation
• Impact of quantum-chemical metrics on the machine learning prediction of electron density

Cited by in corpus (4)

• Machine Learning Potential for Electrochemical Interfaces with Hybrid Representation of Dielectric Response
• Machine Learning and Data-Driven Methods in Computational Surface and Interface Science
• Learning the Electrostatic Response of the Electron Density through a Symmetry-Adapted Vector Field Model
• Predictive Free Energy Simulations Through Hierarchical Distillation of Quantum Hamiltonians