Dielectric Properties of Nanoconfined Water from Ab Initio Thermopotentiostat Molecular Dynamics

Dielectric profile at the Pt(111)/water interface

The dielectric constant, although a simplified concept when considering atomic scales, enters many mean-field, electrochemical interface models and constant potential models as an important parameter. Here, we use ab initio and machine-learned molecular dynamics to scrutinize the behavior of the electronic contribution to ɛr(z) as a function of distance z from a Pt(111) surface. We show that the resulting dielectric profile can largely be explained as a sum of the metallic response and the density-scaled water response at the interface. A slight enhancement of the dielectric response close to the surface can be explained by elongated, strongly polarizable orbitals induced by metal/water bonding. In spite of this enhancement, our results suggest the presence of a region with a very low dielectric constant close to the surface (where the orientational dielectric response does not kick in yet), even for water in contact with hydrophilic metallic interfaces. This region will restrict the double layer capacitance to relatively low values even at potentials where dielectric saturation does not play a role yet. This finding has implications on possible interpretations of double layer capacitances, the dependence of surface electric fields on the ion size, and on electrochemical kinetics.

Revealing the Reaction Pathway of Anodic Hydrogen Evolution at Magnesium Surfaces in Aqueous Electrolytes

Aqueous metal corrosion is a major economic concern in modern society. A phenomenon that has puzzled generations of scientists in this field is the so-called anomalous hydrogen evolution: the violent dissolution of magnesium under electron-deficient (anodic) conditions, accompanied by strong hydrogen evolution and a key mechanism hampering Mg technology. Experimental studies have indicated the presence of univalent Mg+ in solution, but these findings have been largely ignored because they defy our common chemical understanding and evaded direct experimental observation. Using recent advances in the ab initio description of solid-liquid electrochemical interfaces under controlled potential conditions, we describe the full reaction path of Mg atom dissolution from a kinked Mg surface under anodic conditions. Our study reveals the formation of a solvated [Mg2+(OH)-]+ ion complex, challenging the conventional assumption of Mg2+ ion formation. This insight provides an intuitive explanation for the postulated presence of (Coulombically) univalent Mg+ ions, and the absence of protective oxide/hydroxide layers normally formed under anodic/oxidizing conditions. The discovery of this unexpected and unconventional reaction mechanism is crucial for identifying new strategies for corrosion prevention and can be transferred to other metals.

A new method to calculate broadband dielectric spectra of solvents from molecular dynamics simulations demonstrated with polarizable force fields

Simulating the dielectric spectra of solvents requires the nuanced definition of inter- and intra-molecular forces. Non-polarizable force fields, while thoroughly benchmarked for dielectric applications, do not capture all the spectral features of solvents, such as water. Conversely, polarizable force fields have been largely untested in the context of dielectric spectroscopy but include charge and dipole fluctuations that contribute to intermolecular interactions. We benchmark non-polarizable force fields and the polarizable force fields AMOEBA03 and HIPPO for liquid water and find that the polarizable force fields can capture all the experimentally observed spectral features with varying degrees of accuracy. However, the non-polarizable force fields miss at least one peak. To diagnose this deficiency, we decompose the liquid water spectra from polarizable force fields at multiple temperatures into static and induced dipole contributions and find that the peak originates from induced dipole contributions. Broadening our inquiry to other solvents parameterized with the AMOEBA09 force field, we demonstrate good agreement between the experimental and simulated dielectric spectra of methanol and formamide. To produce these spectra, we develop a new computational approach to calculate the dielectric spectrum via the fluctuation dissipation theorem. This method minimizes the error in both the low and high frequency portions of the spectrum, improving the overall accuracy of the simulated spectrum and broadening the computed frequency range.

Nonlinear Effects of Hydrophobic Confinement on the Electronic Structure and Dielectric Response of Water

Fundamental studies of the dielectrics of confined water are critical to understand the ion transport across biological and synthetic nanochannels. The relevance of these fundamental studies, however, surmounts the difficulty of probing water's dielectric constant as a function of a fine variation in confinement. In this work, we explore the computational efficiency of machine learning potentials to derive the confinement effects on the dielectric constant, polarization, and dipole moment of water. Our simulations predict an enhancement of the axial dielectric constant of water under extreme confinement, arising from either the formation of ferroelectric structures of ordered water or larger dipole fluctuations facilitated by the disruption of water's H-bond network. Our study highlights the impact of hydrophobic nanoconfinement on the dielectric constant and on the ionic and electronic structure of water molecules, pointing to the importance of geometric flexibility and electronic polarizability to properly model confinement effects on water.

Deprotonation of formic, acetic acids and bicarbonate ion in slit silica nanopores at infinite dilution and in the presence of electrolytes

Dielectric effects and the coupled electrostatics between the nanoconfined and the internal/external aqueous media contribute to the observed deviations of chemistry within the nanoconfined environment when compared with unconfined systems. A systematic understanding has remained elusive, especially with respect to background salt concentration and boundary condition effects like the nanopore surface chemistry and the reference state used to calculate free energies. We utilize molecular dynamics simulations along with thermodynamic integration to determine the free energy difference associated with acid-base chemistry in 2 nm and 4 nm slit pores open to a bulk-like reservoir. pKa increases are predicted when confining acetic acid, formic acid, and bicarbonate in the slits at infinite dilution conditions. We find that confinement weakens the acids, and the modulation of outer pore surface dipole magnitudes can tune the pKa shift values, suggesting that purely "intrinsic" electrostatic effect on confinement may not exist. At sufficiently high salt concentrations, the dielectric/electrostatic effects on pKa values diminish due to charge screening effects. These discoveries enable future modifications of nanopore chemistries to achieve desirable properties for industrial applications.

Spin Crossover and Exchange Effects on Oxygen Evolution Reaction Catalyzed by Bimetallic Metal Organic Frameworks

Bimetallic metal-organic frameworks (BMOFs) have shown a superior oxygen evolution reaction (OER) performance, attributed to the synergistic effects of dual metal sites. However, the significant role of these dual-metal synergies in the OER is not yet fully understood. In this study, we employed density functional theory to systematically investigate the OER performance of NiAl- and NiFe-based BMOFs by examining all possible spin states of each intermediate across diverse external potentials and pH environments. We found that the spin state featuring a shallow hole trap state and Ni ions with a higher oxidation state serve as strong oxidizing agents, promoting the OER. An external potential-induced spin crossover was observed in each intermediate, resulting in significant changes in the overall reaction and activation energies due to altered energy levels. Combining the constant potential method and the electrochemical nudged elastic band method, we mapped the minimum free energy barriers of the OER under varied external potential and pH by considering the spin crossover effect for both NiAl and NiFe BMOFs. The results showed that NiFe exhibits better OER thermodynamics and kinetics, which is in good agreement with experimentally measured OER polarization curves and Tafel plots. Moreover, we found that the improved OER kinetics of NiFe not only is attributed to lower barriers but also is a result of improved electrical conductivity arising from the synergistic effects of Ni-Fe dual-metal sites. Specifically, replacing the second metal Al with Fe leads to two significant outcomes: a reduction in both the band gap and the effective hole mass compared to NiAl, and the initiation of super- and double-exchange interactions within the Ni-F-Fe chain, thereby enhancing electron transfer and hopping and leading to the improved OER kinetics.

On the Thermodynamic Equivalence of Grand Canonical, Infinite-Size, and Capacitor-Based Models in First-Principle Electrochemistry

First principles-based computational and theoretical methods are constantly evolving trying to overcome the many obstacles towards a comprehensive understanding of electrochemical processes on an atomistic level. One of the major challenges has been the determination of reaction energetics under a constant potential. Here, a theoretical framework was proposed applying standard electronic structure methods and extrapolating to the infinite-cell size limit where reactions do not alter the potential. Today, electronically grand canonical modifications to electronic structure methods, holding the potential constant by varying the number of electrons in a finite simulation cell, become increasingly popular. In this perspective, we show that these two schemes are thermodynamically equivalent. Further, we link these methods to capacitive models of the interface, in the limit that the capacitance of the charging components (whether continuum or atomistic) are equal and invariant along the reaction pathway. We benchmark the three approaches with an example of alkali cation adsorption on Pt(111) showing that all three approaches converge in the cases of Li, Na and K. For Cs, however, strong deviation from the ideal conditions leads to a spread in the respective results. We discuss the latter by highlighting the cases of broken equivalence and assumptions among the approaches.

Activation Energies of Heterogeneous Electrocatalysis: A Theoretical Perspective

Heterogeneous electrochemistry is important for various applications. However, currently, there is limited information about activation energies. In this invited review, we review the challenges associated with calculating these activation energies. Specifically, we highlight three key difficulties in atomistic modeling: liquid structure, electrode potential, and electrolyte ions, along with state-of-the-art methods to address them. We aim to inspire more studies in the field of activation energies to better understand and design heterogeneous electrocatalysts.

An implicit electrolyte model for plane wave density functional theory exhibiting nonlinear response and a nonlocal cavity definition

We have developed and implemented an implicit electrolyte model in the Vienna Ab initio Simulation Package (VASP) that includes nonlinear dielectric and ionic responses as well as a nonlocal definition of the cavities defining the spatial regions where these responses can occur. The implementation into the existing VASPsol code is numerically efficient and exhibits robust convergence, requiring computational effort only slightly higher than the original linear polarizable continuum model. The nonlinear + nonlocal model is able to reproduce the characteristic "double hump" shape observed experimentally for the differential capacitance of an electrified metal interface while preventing "leakage" of the electrolyte into regions of space too small to contain a single water molecule or solvated ion. The model also gives a reasonable prediction of molecular solvation free energies as well as the self-ionization free energy of water and the absolute electron chemical potential of the standard hydrogen electrode. All of this, combined with the additional ability to run constant potential density functional theory calculations, should enable the routine computation of activation barriers for electrocatalytic processes.

Controlled Electrochemical Barrier Calculations without Potential Control

The knowledge of electrochemical activation energies under applied potential conditions is a prerequisite for understanding catalytic activity at electrochemical interfaces. Here, we present a new set of methods that can compute electrochemical barriers with accuracy comparable to that of constant potential grand canonical approaches, without the explicit need for a potentiostat. Instead, we Legendre transform a set of constant charge, canonical reaction paths. Additional straightforward approximations offer the possibility to compute electrochemical barriers at a fraction of computational cost and complexity, and the analytical inclusion of geometric response highlights the importance of incorporating electronic as well as the geometric degrees of freedom when evaluating electrochemical barriers.

Electric Double Layer Effects in Electrocatalysis: Insights from Ab Initio Simulation and Hierarchical Continuum Modeling

Structures of the electric double layer (EDL) at electrocatalytic interfaces, which are modulated by the material properties, the electrolyte characteristics (e.g., the pH, the types and concentrations of ions), and the electrode potential, play crucial roles in the reaction kinetics. Understanding the EDL effects in electrocatalysis has attracted substantial research interest in recent years. However, the intrinsic relationships between the specific EDL structures and electrocatalytic kinetics remain poorly understood, especially on the atomic scale. In this Perspective, we briefly review the recent advances in deciphering the EDL effects mainly in hydrogen and oxygen electrocatalysis through a multiscale approach, spanning from the atomistic scale simulated by ab initio methods to the macroscale by a hierarchical approach. We highlight the importance of resolving the local reaction environment, especially the local hydrogen bond network, in understanding EDL effects. Finally, some of the remaining challenges are outlined, and an outlook for future developments in these exciting frontiers is provided.

Dielectric response of confined water films from a classical density functional theory perspective

We re-examine the problem of the dielectric response of highly polar liquids such as water in confinement between two walls using simple two-variable density functional theory involving number and polarisation densities. In the longitudinal polarisation case where a perturbing field is applied perpendicularly to the walls, we show that the notion of the local dielectric constant, although ill-defined at a microscopic level, makes sense when coarse-graining over the typical size of a particle is introduced. The approach makes it possible to study the effective dielectric response of thin liquid films of various thicknesses in connection with the recent experiments of Fumagalli et al., [Science, 2018, 360, 1339-1342], and to discuss the notion of the interfacial dielectric constant. We argue that the observed properties as a function of slab dimensions, in particular the very low dielectric constants of the order of 2-3 measured for thin slabs of ∼1 nm thickness do not highlight any special properties of water but can be recovered for a generic polar solvent having similar particle size and the same high dielectric constant. Regarding the transverse polarisation case where the perturbing field is parallel to the walls, the associated effective dielectric constant as a function of slab dimensions reaches bulk-like values at much shorter widths than in the longitudinal case. In both cases, we find an oscillatory behaviour for slab thicknesses in the one nanometer range due to packing effects.

Finding Infinities in Nanoconfined Geothermal Electrolyte Static Dielectric Properties and Implications on Ion Adsorption/Pairing

Infinities should naturally occur in the dielectric responses of ionic solutions relevant to many geochemical, energy storage, and electrochemical applications at a strictly zero frequency. Using molecular dynamics simulations cross-referenced with coarse-grained Monte Carlo models, using nanoslit pore models at hydrothermal conditions, and treating confined mobile charges as polarization, we demonstrate the far reaching consequences. The dielectric permittivity profile perpendicular to the slit (ϵ⊥(z)) increases, not decreases, with ionic concentration, unlike in the more widely studied megahertz-to-gigahertz frequency range. In confined electrolytes, the divergences in ϵ⊥(z) correctly describe crossovers between bulk- and surface-dominated dielectric behavior. Nanoconfinement at low ionic concentrations changes monovalent ion energetics by 1-2 kJ/mol, but no dielectric property studied so far is universally correlated to ion adsorption or ion-ion interactions. We caution that infinities signal violation of the "electrical insulator" dielectric assumption.

Grand Canonical Ensemble Modeling of Electrochemical Interfaces Made Simple

Grand canonical ensemble (GCE) modeling of electrochemical interfaces, in which the electrochemical potential is converged to a preset constant, is essential for understanding electrochemistry and electrocatalysis at the electrodes. However, it requires developing efficient and robust algorithms to perform practical and effective GCE modeling with density functional theory (DFT) calculations. Herein, we developed an efficient and robust fully converged constant-potential (FCP) algorithm based on Newton's method and a polynomial fitting to calculate the necessary derivative for DFT calculations. We demonstrated with the constant-potential geometry optimization and Born-Oppenheimer molecular dynamics (BOMD) calculations that our FCP algorithm is resistant to the numerical instability that plagues other algorithms, and it delivers efficient convergence to the preset electrochemical potential and renders accurate forces for updating the nuclear positions of an electronically open system, outperforming other algorithms. The implementation of our FCP algorithm enables flexibility in using various computational codes and versatility in performing advanced tasks including the constant-potential enhanced-sampling BOMD simulations that we showcased with the modeling of the electrochemical hydrogenation of CO, and it is thus expected to find a wide spectrum of applications in the modeling of chemistry at electrochemical interfaces.