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.