Finite field formalism for bulk electrolyte solutions

Impedance of nanocapacitors from molecular simulations to understand the dynamics of confined electrolytes

Nanoelectrochemical devices have become a promising candidate technology across various applications, including sensing and energy storage, and provide new platforms for studying fundamental properties of electrode/electrolyte interfaces. In this work, we employ constant-potential molecular dynamics simulations to investigate the impedance of gold-aqueous electrolyte nanocapacitors, exploiting a recently introduced fluctuation-dissipation relation. In particular, we relate the frequency-dependent impedance of these nanocapacitors to the complex conductivity of the bulk electrolyte in different regimes, and use this connection to design simple but accurate equivalent circuit models. We show that the electrode/electrolyte interfacial contribution is essentially capacitive and that the electrolyte response is bulk-like even when the interelectrode distance is only a few nanometers, provided that the latter is sufficiently large compared to the Debye screening length. We extensively compare our simulation results with spectroscopy experiments and predictions from analytical theories. In contrast to experiments, direct access in simulations to the ionic and solvent contributions to the polarization allows us to highlight their significant and persistent anticorrelation and to investigate the microscopic origin of the timescales observed in the impedance spectrum. This work opens avenues for the molecular interpretation of impedance measurements, and offers valuable contributions for future developments of accurate coarse-grained representations of confined electrolytes.

Computational dielectric spectroscopy on solid-solution interface by time-dependent voltage applied molecular dynamics simulation

A frequency-dependent dielectric constant characterizes the dielectric response of a medium and also represents the time scale of system's collective dynamics. Although it is valuable not only academically but also practically for developing advanced devices, getting the value of a solution at the interface with a solid or electrode surface is challenging both experimentally and computationally. Here, we propose a computational method that imitates the dielectric spectroscopy and AC impedance measurement. It combines a time-dependent voltage applied molecular dynamics simulation with an equivalent circuit representation of a system composed of a solution confined between two identical electrodes. It gives the frequency-dependent dielectric constants of the bulk solution and the interface simultaneously. Unlike the conventional method, it does not require computation of a dipole autocorrelation function and its Fourier transformation. Application of the method on a system of water confined between polarizable Pt electrodes gives the static dielectric constant and the relaxation time of the bulk water in good agreement with previous simulation results and experimental values. In addition, it gives a much smaller static dielectric constant at the interface, consistent with previous observations. The outline of the dielectric dispersion curve of the interface seems similar to that of the bulk, but the relaxation time is several times faster.

The limit of macroscopic homogeneous ice nucleation at the nanoscale

Nucleation in small volumes of water has garnered renewed interest due to the relevance of pore condensation and freezing under conditions of low partial pressures of water, such as in the upper troposphere. Molecular simulations can in principle provide insight on this process at the molecular scale that is challenging to achieve experimentally. However, there are discrepancies in the literature as to whether the rate in confined systems is enhanced or suppressed relative to bulk water at the same temperature and pressure. In this study, we investigate the extent to which the size of the critical nucleus and the rate at which it grows in thin films of water are affected by the thickness of the film. Our results suggest that nucleation remains bulk-like in films that are barely large enough accommodate a critical nucleus. This conclusion seems robust to the presence of physical confining boundaries. We also discuss the difficulties in unambiguously determining homogeneous nucleation rates in nanoscale systems, owing to the challenges in defining the volume. Our results suggest any impact on a film's thickness on the rate is largely inconsequential for present day experiments.

Electrical noise in electrolytes: a theoretical perspective

Seemingly unrelated experiments such as electrolyte transport through nanotubes, nano-scale electrochemistry, NMR relaxometry and surface force balance measurements, all probe electrical fluctuations: of the electric current, the charge and polarization, the field gradient (for quadrupolar nuclei) and the coupled mass/charge densities. The fluctuations of such various observables arise from the same underlying microscopic dynamics of the ions and solvent molecules. In principle, the relevant length and time scales of these dynamics are encoded in the dynamic structure factors. However, modelling the latter for frequencies and wavevectors spanning many orders of magnitude remains a great challenge to interpret the experiments in terms of physical processes such as solvation dynamics, diffusion, electrostatic and hydrodynamic interactions between ions, interactions with solid surfaces, etc. Here, we highlight the central role of the charge-charge dynamic structure factor in the fluctuations of electrical observables in electrolytes and offer a unifying perspective over a variety of complementary experiments. We further analyze this quantity in the special case of an aqueous NaCl electrolyte, using simulations with explicit ions and an explicit or implicit solvent. We discuss the ability of the standard Poisson-Nernst-Planck theory to capture the simulation results, and how the predictions can be improved. We finally discuss the contributions of ions and water to the total charge fluctuations. This work illustrates an ongoing effort towards a comprehensive understanding of electrical fluctuations in bulk and confined electrolytes, in order to enable experimentalists to decipher the microscopic properties encoded in the measured electrical noise.

Ionic fluctuations in finite volumes: fractional noise and hyperuniformity

Observing finite regions of a bigger system is a common aim, from microscopy to molecular simulations. In the latter especially, there is ongoing interest in predicting thermodynamic properties from tracking fluctuations in finite observation volumes. However, kinetic properties have received little attention, especially not in ionic solutions, where electrostatic interactions play a decisive role. Here, we probe ionic fluctuations in finite volumes with Brownian dynamics and build an analytical framework that reproduces our simulation results and is broadly applicable to other systems with pairwise interactions. Particle number and charge correlations exhibit a rich phenomenology with time, characterized by a diversity of timescales. The noise spectrum of both quantities decays as 1/f3/2, where f is the frequency. This signature of fractional noise shows the universality of 1/f3/2 scalings when observing diffusing particles in finite domains. The hyperuniform behaviour of charge fluctuations, namely that correlations scale with the area of the observation volume, is preserved in time. Correlations even become proportional to the box perimeter at sufficiently long times. Our results pave the way to understand fluctuations in more complex systems, from nanopores to single-particle electrochemistry.

Simulating dielectric spectra: A demonstration of the direct electric field method and a new model for the nonlinear dielectric response

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.

A theory for the stabilization of polar crystal surfaces by a liquid environment

Polar crystal surfaces play an important role in the functionality of many materials and have been studied extensively over many decades. In this article, a theoretical framework is presented that extends existing theories by placing the surrounding solution environment on an equal footing with the crystal itself; this is advantageous, e.g., when considering processes such as crystal growth from solution. By considering the polar crystal as a stack of parallel plate capacitors immersed in a solution environment, the equilibrium adsorbed surface charge density is derived by minimizing the free energy of the system. In analogy to the well-known diverging surface energy of a polar crystal surface at zero temperature, for a crystal in solution it is shown that the "polar catastrophe" manifests as a diverging free energy cost to perturb the system from equilibrium. Going further than existing theories, the present formulation predicts that fluctuations in the adsorbed surface charge density become increasingly suppressed with increasing crystal thickness. We also show how, in the slab geometry often employed in both theoretical and computational studies of interfaces, an electric displacement field emerges as an electrostatic boundary condition, the origins of which are rooted in the slab geometry itself, rather than the use of periodic boundary conditions. This aspect of the work provides a firmer theoretical basis for the recent observation that standard "slab corrections" fail to correctly describe, even qualitatively, polar crystal surfaces in solution.

Can molecular simulations reliably compare homogeneous and heterogeneous ice nucleation?

In principle, the answer to the posed titular question is undoubtedly "yes." But in practice, requisite reference data for homogeneous systems have been obtained with a treatment of intermolecular interactions that is different from that typically employed for heterogeneous systems. In this article, we assess the impact of the choice of truncation scheme when comparing water in homogeneous and inhomogeneous environments. Specifically, we use explicit free energy calculations and a simple mean field analysis to demonstrate that using the "cut-and-shift" version of the Lennard-Jones potential (common to most simple point charge models of water) results in a systematic increase in the melting temperature of ice I_h. In addition, by drawing an analogy between a change in cutoff and a change in pressure, we use existing literature data for homogeneous ice nucleation at negative pressures to suggest that enhancements due to heterogeneous nucleation may have been overestimated by several orders of magnitude.

Transport coefficients of gel electrolytes: A molecular dynamics simulation study

The responses of gel electrolytes to stimuli make them useful in applications such as sensors and actuators. However, few studies have explored their transport properties from a molecular viewpoint. We studied the transport coefficients of gel electrolytes based on perfluorinated sulfonic acid using molecular dynamics simulations. The transport coefficients for electric and pressure fields, namely, the ionic conductivity, Darcy permeability, and cross coupling constant, were calculated based on Kubo’s linear response theory from the corresponding velocity correlation functions and mean square displacements. The effects of the water content of the gel electrolyte and those of the monovalent cationic species were also analyzed. The calculated transport coefficients qualitatively agree with the reported experimental results. The role of the cross coupling constants in determining the functional efficiency of gel electrolytes as pressure sensors or electroactive actuators is discussed.

Dielectric response of thin water films: a thermodynamic perspective

The surface of a polar liquid presents a special environment for the solvation and organization of charged solutes, which differ from bulk behaviors in important ways. These differences have motivated many attempts to understand electrostatic response at aqueous interfaces in terms of a spatially varying dielectric permittivity, typically concluding that the dielectric constant of interfacial water is significantly lower than in the bulk liquid. Such analyses, however, are complicated by the potentially nonlocal nature of dielectric response over the short length scales of interfacial heterogeneity. Here we circumvent this problem for thin water films by adopting a thermodynamic approach. Using molecular simulations, we calculate the solvent's contribution to the reversible work of charging a parallel plate capacitor. We find good agreement with a simple dielectric continuum model that assumes bulk dielectric permittivity all the way up to the liquid's boundary, even for very thin (∼1 nm) films. This comparison requires careful attention to the placement of dielectric boundaries between liquid and vapor, which also resolves apparent discrepancies with dielectric imaging experiments.

Free energy calculations from molecular simulations reveal that water's interfacial dielectric response is well-described by bulk properties.

Computational Amperometry of Nanoscale Capacitors in Molecular Simulations

In recent years, constant applied potential molecular dynamics has allowed researchers to study the structure and dynamics of the electrochemical double-layer of a large variety of nanoscale capacitors. Nevertheless, it has remained impossible to simulate polarized electrodes at fixed total charge. Here, we show that combining a constant potential electrode with a finite electric displacement fills this gap by allowing us to simulate open-circuit conditions. The method can be extended by applying an electric displacement ramp to perform computational amperometry experiments at different current intensities. As in experiments, the full capacitance of the system is obtained at low intensity, but this quantity decreases when the applied ramp becomes too fast with respect to the microscopic dynamics of the liquid.

Electric-field-based Poisson-Boltzmann theory: Treating mobile charge as polarization

Mobile charge in an electrolytic solution can in principle be represented as the divergence of ionic polarization. After adding explicit solvent polarization a finite volume of an electrolyte can then be treated as a composite nonuniform dielectric body. Writing the electrostatic interactions as an integral over electric-field energy density we show that the Poisson-Boltzmann functional in this formulation is convex and can be used to derive the equilibrium equations for electric potential and ion concentration by a variational procedure developed by Ericksen for dielectric continua [J. L. Ericksen, Arch. Rational Mech. Anal. 183, 299 (2007)AVRMAW0003-952710.1007/s00205-006-0042-4]. The Maxwell field equations are enforced by extending the set of variational parameters by a vector potential representing the dielectric displacement which is fully transverse in a dielectric system without embedded external charge. The electric-field energy density in this representation is a function of the vector potential and the sum of ionic and solvent polarization making the mutual screening explicit. Transverse polarization is accounted for by construction, lifting the restriction to longitudinal polarization inherent in the electrostatic potential based formulation of Poisson-Boltzmann mean field theory.

Macroscopic surface charges from microscopic simulations

Attaining accurate average structural properties in a molecular simulation should be considered a prerequisite if one aims to elicit meaningful insights into a system's behavior. For charged surfaces in contact with an electrolyte solution, an obvious example is the density profile of ions along the direction normal to the surface. Here, we demonstrate that, in the slab geometry typically used in simulations, imposing an electric displacement field D determines the integrated surface charge density of adsorbed ions at charged interfaces. This allows us to obtain macroscopic surface charge densities irrespective of the slab thickness used in our simulations. We also show that the commonly used Yeh-Berkowitz method and the "mirrored slab" geometry both impose vanishing integrated surface charge densities. We present results both for relatively simple rocksalt (1 1 1) interfaces and the more complex case of kaolinite's basal faces in contact with an aqueous electrolyte solution.

Dielectric response with short-ranged electrostatics

The dielectric nature of polar liquids underpins much of their ability to act as useful solvents, but its description is complicated by the long-ranged nature of dipolar interactions. This is particularly pronounced under the periodic boundary conditions commonly used in molecular simulations. In this article, the dielectric properties of a water model whose intermolecular electrostatic interactions are entirely short-ranged are investigated. This is done within the framework of local molecular-field theory (LMFT), which provides a well-controlled mean-field treatment of long-ranged electrostatics. This short-ranged model gives a remarkably good performance on a number of counts, and its apparent shortcomings are readily accounted for. These results not only lend support to LMFT as an approach for understanding solvation behavior, but also are relevant to those developing interaction potentials based on local descriptions of liquid structure.

Frontiers in molecular simulation of solvated ions, molecules and interfaces

This themed collection is a collection of articles on frontiers in molecular simulation of solvated ions, molecules and interfaces.