A major challenge in the molecular simulation of electric double layer capacitors (EDLCs) is the choice of an appropriate model for the electrode. Typically, in such simulations the electrode surface is modeled using a uniform fixed charge on each of
the electrode atoms, which ignores the electrode response to local charge fluctuations induced by charge fluctuations in the electrolyte. In this work, we evaluate and compare this Fixed Charge Method (FCM) with the more realistic Constant Potential Method (CPM), [Reed, et al., J. Chem. Phys., 126, 084704 (2007)], in which the electrode charges fluctuate in order to maintain constant electric potential in each electrode. For this comparison, we utilize a simplified LiClO$_4$-acetonitrile/graphite EDLC. At low potential difference ($DeltaPsile 2V$), the two methods yield essentially identical results for ion and solvent density profiles; however, significant differences appear at higher $DeltaPsi$. At $DeltaPsige 4V$, the CPM ion density profiles show significant enhancement (over FCM) of partially electrode solvated Li$^+$ ions very close to the electrode surface. The ability of the CPM electrode to respond to local charge fluctuations in the electrolyte is seen to significantly lower the energy (and barrier) for the approach of Li$^+$ ions to the electrode surface.
Structural transformations at interfaces are of profound fundamental interest as complex examples of phase transitions in low-dimensional systems. Despite decades of extensive research, no compelling evidence exists for structural transformations in
high-angle grain boundaries in elemental systems. Here we show that the critical impediment to observations of such phase transformations in atomistic modeling has been rooted in inadequate simulation methodology. The proposed new methodology allows variations in atomic density inside the grain boundary and reveals multiple grain boundary phases with different atomic structures. Reversible first-order transformations between such phases are observed by varying temperature or injecting point defects into the boundary region. Due to the presence of multiple metastable phases, grain boundaries can absorb significant amounts of point defects created inside the material by processes such as irradiation. We propose a novel mechanism of radiation damage healing in metals which may guide further improvements in radiation resistance of metallic materials through grain boundary engineering.
