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

Andrea Grisafi,Mathieu Salanne

2024-05-13

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.

Materials Science

What problem does this paper attempt to address?

The purpose of this paper is to address the efficiency and accuracy issues in calculating the electron charge distribution on simulated electrochemical interfaces. Traditional quantum mechanics/molecular mechanics (QM/MM) methods often limit large-scale simulations of electrochemical systems due to high computational costs. In this paper, the authors propose a machine learning (ML) approach, particularly using the SALTED algorithm, to predict the electron charge density of metal electrodes and integrate it with the MetalWalls program for classical simulations of electrochemical systems. The authors develop a spherical harmonic expansion-based Ewald summation method in the paper, which effectively calculates the electric field generated by the predicted charge distribution and accurately calculates the electrostatic forces on electrolyte atoms. This method avoids the traditional partial charge-based approaches and demonstrates the effectiveness of predicting the electrostatics of the system by predicting the electron density pathway. Through data-driven molecular dynamics simulations of approximately 3 nanoseconds, the study reveals significant differences between the distribution at the electrolyte interface and the results from pure classical simulations. The main contribution of the paper is the acceleration of QM/MM methods in electrochemical systems, making it possible to achieve accurate atomic-level descriptions of the electric double layer at the nanosecond timescale. This is crucial for understanding the interactions between metal electrodes and liquid electrolytes, and has potential applications in electrochemistry, electrocatalysis, and energy storage.

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

MetalWalls: Simulating electrochemical interfaces between polarizable electrolytes and metallic electrodes

Scalable molecular simulation of electrolyte solutions with quantum chemical accuracy

MetalWalls: A Classical Molecular Dynamics Software Dedicated to the Simulation of Electrochemical Systems

Predicting the Charge Density Response in Metal Electrodes

Machine Learning Potential for Electrochemical Interfaces with Hybrid Representation of Dielectric Response

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

Advancing Multiscale Molecular Modeling with Machine Learning-Derived Electrostatics

Electrostatic Embedding Machine Learning for Ground and Excited State Molecular Dynamics of Solvated Molecules

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

Artificial Intelligence and QM/MM with a Polarizable Reactive Force Field for Next-Generation Electrocatalysts

Machine Learning in QM/MM Molecular Dynamics Simulations of Condensed-Phase Systems

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

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

An Electrostatically Embedded QM/MM Scheme for Electrified Interfaces

Development and application of hybrid AIMD/cDFT simulations for atomic-to-mesoscale chemistry

Analytic gradients for the electrostatic embedding QM/MM in periodic boundary conditions using particle-mesh Ewald sums and electrostatic potential fitted charge operators

Molecular Simulation of Electrode-Solution Interfaces

Advancement of the Homogeneous Background Method for the Computational Simulation of Electrochemical Interfaces

Hybrid Quantum Mechanical, Molecular Mechanical, and Machine Learning Potential for Computing Aqueous-Phase Adsorption Free Energies on Metal Surfaces

Molecular modelling and simulation of electrolyte solutions, biomolecules, and wetting of component surfaces