RESUMO
Electrochemistry is central to many applications, ranging from biology to energy science. Studies now involve a wide range of techniques, both experimental and theoretical. Modeling and simulations methods, such as density functional theory or molecular dynamics, provide key information on the structural and dynamic properties of the systems. Of particular importance are polarization effects of the electrode/electrolyte interface, which are difficult to simulate accurately. Here, we show how these electrostatic interactions are taken into account in the framework of the Ewald summation method. We discuss, in particular, the formal setup for calculations that enforce periodic boundary conditions in two directions, a geometry that more closely reflects the characteristics of typical electrolyte/electrode systems and presents some differences with respect to the more common case of periodic boundary conditions in three dimensions. These formal developments are implemented and tested in MetalWalls, a molecular dynamics software that captures the polarization of the electrolyte and allows the simulation of electrodes maintained at a constant potential. We also discuss the technical aspects involved in the calculation of two sets of coupled degrees of freedom, namely the induced dipoles and the electrode charges. We validate the implementation, first on simple systems, then on the well-known interface between graphite electrodes and a room-temperature ionic liquid. We finally illustrate the capabilities of MetalWalls by studying the adsorption of a complex functionalized electrolyte on a graphite electrode.
RESUMO
In situ NMR spectroscopy is a powerful technique to investigate charge storage mechanisms in carbon-based supercapacitors thanks to its ability to distinguish ionic and molecular species adsorbed in the porous electrodes from those in the bulk electrolyte. The NMR peak corresponding to the adsorbed species shows a clear change of chemical shift as the applied potential difference is varied. This variation in chemical shift is thought to originate from a combination of ion reorganisation in the pores and changes in ring current shifts due to the changes of electronic density in the carbon. While previous Density Functional Theory calculations suggested that the electronic density has a large effect, the relative contributions of these two effects is challenging to untangle. Here, we use mesoscopic simulations to simulate NMR spectra and investigate the relative importance of ion reorganisation and ring currents on the resulting chemical shift. The model is able to predict chemical shifts in good agreement with NMR experiments and indicates that the ring currents are the dominant contribution. A thorough analysis of a specific electrode/electrolyte combination for which detailed NMR experiments have been reported allows us to confirm that local ion reorganisation has a very limited effect but the relative quantities of ions in pores of different sizes, which can change upon charging/discharging, can lead to a significant effect. Our findings suggest that in situ NMR spectra of supercapacitors may provide insights into the electronic structure of carbon materials in the future.
