Structure and dynamics of electrical double layers in organic electrolytes

The organic electrolyte of tetraethylammonium tetrafluoroborate (TEABF(4)) in the aprotic solvent of acetonitrile (ACN) is widely used in electrochemical systems such as electrochemical capacitors. In this paper, we examine the solvation of TEA(+) and BF(4)(-) in ACN, and the structure, capacitance, and dynamics of the electrical double layers (EDLs) in the TEABF(4-)ACN electrolyte using molecular dynamics simulations complemented with quantum density functional theory calculations. The solvation of TEA(+) and BF(4)(-) ions is found to be much weaker than that of small inorganic ions in aqueous solutions, and the ACN molecules in the solvation shell of both types of ions show only weak packing and orientational ordering. These solvation characteristics are caused by the large size, charge delocalization, and irregular shape (in the case of TEA(+) cation) of the ions. Near neutral electrodes, the double-layer structure in the organic electrolyte exhibits a rich organization: the solvent shows strong layering and orientational ordering, ions are significantly contact-adsorbed on the electrode, and alternating layers of cations/anions penetrate ca. 1.1 nm into the bulk electrolyte. The significant contact adsorption of ions and the alternating layering of cation/anion are new features found for EDLs in organic electrolytes. These features essentially originate from the fact that van der Waals interactions between organic ions and the electrode are strong and the partial desolvation of these ions occurs easily, as a result of the large size of the organic ions. Near charged electrodes, distinct counter-ion concentration peaks form, and the ion distribution cannot be described by the Helmholtz model or the Helmholtz + Poisson-Boltzmann model. This is because the number of counter-ions adsorbed on the electrode exceeds the number of electrons on the electrode, and the electrode is over-screened in parts of the EDL. The computed capacitances of the EDLs are in good agreement with that inferred from experimental measurements. Both the rotations (ACN only) and translations of interfacial ACN and ions are found to slow down as the electrode is electrified. We also observe an asymmetrical dependence of these motions on the sign of the electrode charge. The rotation/diffusion of ACN and the diffusion of ions in the region beyond the first ACN or ion layer differ only weakly from those in the bulk.