Abstract:Atomistic modeling of electrified interfaces remains a major issue for detailed insights in electrocatalysis, corrosion, electrodeposition, batteries, and related devices such as pseudocapacitors. In these domains, the use of grand-canonical density functional theory (GC-DFT) in combination with implicit solvation models has become popular. GC-DFT can be conveniently applied not only to metallic surfaces but also to semiconducting oxides and sulfides and is, furthermore, sufficiently robust to achieve a consistent description of reaction pathways. However, the accuracy of implicit solvation models for solvation effects at interfaces is in general unknown. One promising way to overcome the limitations of implicit solvents is going toward hybrid quantum mechanical (QM)/molecular mechanics (MM) models. For capturing the electrochemical potential dependence, the key quantity is the capacitance, i.e., the relation between the surface charge and the electrochemical potential. In order to retrieve the electrochemical potential from a QM/MM hybrid scheme, an electrostatic embedding is required. Furthermore, the charge of the surface and of the solvent regions has to be strictly opposite in order to consistently simulate charge-neutral unit cells in MM and in QM. To achieve such a QM/MM scheme, we present the implementation of electrostatic embedding in the VASP code. This scheme is broadly applicable to any neutral or charged solid/liquid interface. Here, we demonstrate its use in the context of GC-DFT for the hydrogen evolution reaction (HER) over a noble-metal-free electrocatalyst, MoS2. We investigate the effect of electrostatic embedding compared to the implicit solvent model for three contrasting active sites on MoS2: (i) the sulfur vacancy defect, which is rather apolar; (ii) a Mo antisite defect, where the active site is a surface bound highly polar OH group; and (iii) a reconstructed edge site, which is generally believed to be responsible for most of the catalytic activity. According to our results, the electrostatic embedding leads to almost indistinguishable results compared to the implicit solvent for the apolar system but has a significant effect on polar sites. This demonstrates the reliability of the hybrid QM/MM, electrostatically embedded solvation model for electrified interfaces.

Modeling solid-liquid interface reactions with next generation extended Lagrangian quantum-based molecular dynamics

Quantum mechanical and molecular mechanical simulation approaches bridging length and time scales for simulation of interface reactions in realistic environments

Extended Lagrangian Born-Oppenheimer Molecular Dynamics with DFT+U

Electronic Structure Modeling of Electrochemical Reactions at Electrode/Electrolyte Interfaces in Lithium Ion Batteries

First-principles molecular dynamics simulations at solid-liquid interfaces with a continuum solvent

Bridging the gap between mesoscopic and molecular models of solid/liquid interfaces out-of-equilibrium

(Invited) Study of Electrochemical Interfaces Using DFT and the Classical Liquid Theory Hybrid Simulations

Nonadiabatic Dynamics of Molecules Interacting with Metal Surfaces: A Quantum-Classical Approach Based on Langevin Dynamics and the Hierarchical Equations of Motion

Microenvironment effects from first principles multiscale modeling of electrochemical CO2 reduction

Hybrid Combination of Quantum Mechanics with Quantum-Based Polarizable Reactive Force Field for Large Scale Full Solvent Simulations of Electrocatalysis

Molecular Simulation of Electrode-Solution Interfaces

Modeling Operando Electrochemical CO 2 Reduction

Reaction intermediates during operando electrocatalysis identified from full solvent quantum mechanics molecular dynamics

Catalytic flow with a coupled Finite Difference -- Lattice Boltzmann scheme

Simulation of potential and species distribution in a Li||Bi liquid metal battery using coupled meshes

An Electrostatically Embedded QM/MM Scheme for Electrified Interfaces

(Invited) Oxide-Water Interfaces from Ab Initio Based Deep Neural Network Simulations

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

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

Simulation of lithium hydroxide decomposition using deep potential molecular dynamics

Modeling of a non-aqueous Li-O2 battery incorporating synergistic reaction mechanisms, microstructure, and species transport in the porous electrode