An evaluation of the capacitive behavior of supercapacitors as a function of the radius of cations using simulations with a constant potential method

Theoretical description of the electrical double layer for a mixture of n ionic species with arbitrary size and charge asymmetries. I. Spherical geometry

In this work, we propose a theoretical finite element description of the ionic profiles of a general mixture of n species of spherical charged particles dissolved in an implicit solvent, with arbitrary size and charge asymmetries, neutralizing a spherical macroion. This approach aims to close the gap between the nano- and the micro-scales in macroion solutions, taking into account the ion correlations and ionic excluded volume effects consistently. When these last two features are neglected, the classical non-linear Poisson-Boltzmann theory for n ionic species-with different ionic closest approach distances to the colloidal surface-is recovered as a limit case. As a proof of concept, we study the electrical double layer of an electroneutral mixture of oppositely charged colloids and small microions, with an asymmetry 1:333 in size and 1:10 in valence, in salt-free and added salt environments. Our theoretical approach displays a good agreement regarding the ionic profiles, the integrated charge, and the mean electrostatic potential obtained from molecular dynamics simulations with explicit-sized microions. Although the non-linear Poisson-Boltzmann colloid-colloid and colloid-microion profiles differ notably from those obtained via molecular dynamics simulations with explicit small-sized ions, the associated mean electrostatic potential agrees well with the corresponding explicit microion simulations.

Molecular dynamics of preferential adsorption in mixed alkali–halide electrolytes at graphene electrodes

Understanding the microscopic behavior of aqueous electrolyte solutions in contact with graphene and related carbon surfaces is important in electrochemical technologies, such as capacitive deionization or supercapacitors. In this work, we focus on preferential adsorption of ions in mixed alkali–halide electrolytes containing different fractions of Li+/Na+ or Li+/K+ and/or Na+/K+ cations with Cl− anions dissolved in water. We performed molecular dynamics simulations of the solutions in contact with both neutral and positively and negatively charged graphene surfaces under ambient conditions, using the effectively polarizable force field. The simulations show that large ions are often intuitively attracted to oppositely charged electrodes. In contrast, the adsorption behavior of small ions tends to be counterintuitive. In mixed-cation solutions, one of the cations always supports the adsorption of the other cation, while the other cation weakens the adsorption of the first cation. In mixed-cation solutions containing large and small cations simultaneously, adsorption of the larger cations varies dramatically with the electrode charge in an intuitive way, while adsorption of the smaller cations changes oppositely, i.e., in a counterintuitive way. For (Li/K)Cl mixed-cation solutions, these effects allow the control of Li+ adsorption by varying the electrode charge, whereas, for LiCl single-salt solutions, Li+ adsorption is nearly independent of the electrode charge. We rationalize this cation–cation lever effect as a result of a competition between three driving forces: (i) direct graphene–ion interactions, (ii) the strong tendency of the solutions to saturate the network of non-covalent intermolecular bonds, and (iii) the tendency to suppress local charge accumulation in any region larger than typical interparticle distances. We analyze the driving forces in detail using a general method for intermolecular bonding based on spatial distribution functions and different contributions to the total charge density profiles. The analysis helps to predict whether an ion is more affected by each of the three driving forces, depending on the strength of the ion solvation shells and the compatibility between the contributions of the charge density profiles due to the ion and water molecules. This approach is general and can also be applied to other solutions under different thermodynamic conditions.