Search

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. [ABSTRACT FROM AUTHOR]

Controlling electrochemical reactivity requires a detailed understanding of the charging behavior and thermodynamics of the electrochemical interface. Experiments can independently probe the overall charge response of the electrochemical double layer by capacitance measurements and the thermodynamics of the inner layer with potential of maximum entropy measurements. Relating these properties by computational modeling of the electrochemical interface has so far been challenging due to the low accuracy of classical molecular dynamics (MD) for capacitance and the limited time and length scales of ab initio MD. Here, we combine large ensembles of long-time-scale classical MD simulations with charge response from electronic density functional theory to predict the potential-dependent capacitance of a family of ideal aqueous electrochemical interfaces with different peak capacitances. We show that while the potential of maximum capacitance varies, this entire family exhibits an electrode charge of maximum capacitance (CMC) between −2.9 and −2.2 μC/cm2, regardless of the details in the electronic response. Simulated heating of the same interfaces reveals that the entropy peaks at a charge of maximum entropy (CME) of −5.1 ± 0.6 μC/cm2, in agreement with experimental findings for metallic electrodes. The CME and CMC both indicate asymmetric response of interfacial water that is stronger for negatively charged electrodes, while the difference between CME and CMC illustrates the richness in behavior of even the ideal electrochemical interface. [ABSTRACT FROM AUTHOR]

In those MD simulations, the NVT ensemble was erroneously applied to all atoms, including the frozen electrode atoms, instead of the electrolyte atoms alone. Electrode charge densities (41 curves on the left from blue for -16 C/cm2 to red for +16 C/cm2) for each of the charge models and placement schemes at z = 0 in Fig. The potential of maximum capacitance (PMC) varies from -0.09 to -0.22 V from the PZC across the family, but the charge of maximum capacitance (CMC, black squares) is constant at -3.7 C/cm2 in (a)-(c) and -3.3 C/cm2 in (d). [Extracted from the article]

The microscopic dielectric response is a key quantity for electronic materials such as organic semiconductors. Calculations of this response for molecular crystals are currently either expensive or rely on extreme simplifications such as multipole expansions which lack microscopic detail. We present an alternate approach using a microscopic analogue of the Clausius-Mossotti equation, which constructs the dielectric response of a crystal from an eigenvalue decomposition of the dielectric response of individual molecules. This method can potentially be used to examine the effects of defects, disorder, and surfaces on the dielectric properties of molecular solids. [ABSTRACT FROM AUTHOR]

Continuum solvation models enable electronic structure calculations of systems in liquid environments, but because of the large number of empirical parameters, they are limited to the class of systems in their fit set (typically organic molecules). Here, we derive a solvation model with no empirical parameters for the dielectric response by taking the linear response limit of a classical density functional for molecular liquids. This model directly incorporates the nonlocal dielectric response of the liquid using an angular momentum expansion, and with a single fit parameter for dispersion contributions it predicts solvation energies of neutral molecules with a RMS error of 1.3 kcal/mol in water and 0.8 kcal/mol in chloroform and carbon tetrachloride. We show that this model is more accurate for strongly polar and charged systems than previous solvation models because of the parameter-free electric response, and demonstrate its suitability for ab initio solvation, including self-consistent solvation in quantum Monte Carlo calculations. [ABSTRACT FROM AUTHOR]

Reliable first-principles calculations of electrochemical processes require accurate prediction of the interfacial capacitance, a challenge for current computationally efficient continuum solvation methodologies. We develop a model for the double layer of a metallic electrode that reproduces the features of the experimental capacitance of Ag(100) in a non-adsorbing, aqueous electrolyte, including a broad hump in the capacitance near the potential of zero charge and a dip in the capacitance under conditions of low ionic strength. Using this model, we identify the necessary characteristics of a solvation model suitable for first-principles electrochemistry of metal surfaces in non-adsorbing, aqueous electrolytes: dielectric and ionic nonlinearity, and a dielectric-only region at the interface. The dielectric nonlinearity, caused by the saturation of dipole rotational response in water, creates the capacitance hump, while ionic nonlinearity, caused by the compactness of the diffuse layer, generates the capacitance dip seen at low ionic strength. We show that none of the previously developed solvation models simultaneously meet all these criteria. We design the nonlinear electrochemical soft-sphere solvation model which both captures the capacitance features observed experimentally and serves as a general-purpose continuum solvation model. [ABSTRACT FROM AUTHOR]

We demonstrate a method to compute the dielectric spectra of fluids in molecular dynamics (MD) by directly applying electric fields to the simulation. We obtain spectra from MD simulations with low magnitude electric fields (≈0.01 V/Å) in agreement with spectra from the fluctuation-dissipation method for water and acetonitrile. We examine this method's trade-off between noise at low field magnitudes and the nonlinearity of the response at higher field magnitudes. We then apply the Booth equation to describe the nonlinear response of both fluids at low frequency (0.1 GHz) and high field magnitude (up to 0.5 V/Å). We develop a model of the frequency-dependent nonlinear response by combining the Booth description of the static nonlinear dielectric response of fluids with the frequency-dependent linear dielectric response of the Debye model. We find good agreement between our model and the MD simulations of the nonlinear dielectric response for both acetonitrile and water.

Fully harnessing electrochemical interfaces for reactions requires a detailed understanding of solvent effects in the electrochemical double layer. Predicting the significant impact of solvents on entropic and electronic properties of electrochemical interfaces has remained an open challenge of computational electrochemistry. Using molecular dynamics simulations of silver-water and silver-acetonitrile interfaces, we show that switching the solvent changes the signs for both the charge of maximum capacitance (CMC) and charge of maximum entropy (CME). Contrasting the capacitance and CME behavior of these two interfaces, we demonstrate that the preferred orientation of the solvent molecule and the corresponding charge density determine the sign of the CMC and CME and, hence, the qualitatively different charge asymmetry of the electrochemical interface.

Ab initio modeling of electrochemical systems is becoming a key tool for understanding and predicting electrochemical behavior. Development and careful benchmarking of computational electrochemical methods are essential to ensure their accuracy. Here, using charging curves for an electrode in the presence of an inert aqueous electrolyte, we demonstrate that most continuum models, which are parameterized and benchmarked for molecules, anions, and cations in solution, undersolvate metal surfaces, and underestimate the surface charge as a function of applied potential. We examine features of the electrolyte and interface that are captured by these models and identify improvements necessary for realistic electrochemical calculations of metal surfaces. Finally, we reparameterize popular solvation models using the surface charge of Ag(100) as a function of voltage to find improved accuracy for metal surfaces without significant change in utility for molecular and ionic solvation.