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.

A fully-polarizable KS-DFT/AMOEBA embedding scheme for plane wave basis sets through the MiMiC framework

Wavefunction-based electrostatic-embedding QM/MM using CFOUR through MiMiC

Tinktep: A Fully Self-Consistent, Mutually Polarizable Qm/Mm Approach Based On The Amoeba Force Field

Self-Consistent Solution Of Kohn-Sham Equation By Real-Space Finite-Difference Method

Multistate, polarizable QM/MM embedding scheme based on the direct reaction field method: Solvatochromic shifts, analytical gradients and optimization of conical intersections in solution

A polarizable CASSCF/MM approach using the interface between OpenMMPol library and CFour

Multistate, Polarizable QM/MM Embedding Scheme Based on the Direct Reaction Field Method: Solvatochromic Shifts, Analytical Gradients and Optimizations of Conical Intersections in Solution

Interfacing the Core-Shell or the Drude Polarizable Force Field With Car-Parrinello Molecular Dynamics for QM/MM Simulations

Seamless integration of GEM, a density based-force field, for QM/MM simulations via LICHEM, Psi4 and Tinker-HP

MiMiC: A High-Performance Framework for Multiscale Molecular Dynamics Simulations

A Quantum-Mechanics Molecular-Mechanics scheme for extended systems

A Polarizable QM/MM Model That Combines the State-Averaged CASSCF and AMOEBA Force Field for Photoreactions in Proteins

An Electrostatically Embedded QM/MM Scheme for Electrified Interfaces

The OpenMMPol Library for Polarizable QM/MM Calculations of Properties and Dynamics

PyDFT-QMMM: A modular, extensible software framework for DFT-based QM/MM molecular dynamics

Amber Free Energy Tools: Interoperable Software for Free Energy Simulations Using Generalized Quantum Mechanical/molecular Mechanical and Machine Learning Potentials.

Performance of the AMOEBA Water Model in the Vicinity of QM Solutes: A Diagnosis Using Energy Decomposition Analysis.

Completely Multipolar Model for Many-Body Water-Ion and Ion-Ion Interactions

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

Assessment of embedding schemes in a hybrid machine learning/classical potentials (ML/MM) scheme

Effect of the QM Size, Basis Set, and Polarization on QM/MM Interaction Energy Decomposition Analysis