Comparison of the performances of different computational methods to calculate the electrochemical stability of selected ionic liquids

The electrochemical stability windows (ESW) of selected ionic liquids have been calculated by comparing different computational approaches previously suggested in the literature. The molecular systems under study are based on di-alkyl imidazolium and tetra-alkyl ammonium cations coupled with two different imide anions (namely, bis-fluorosulfonyl imide and bis-trifluoromethyl sulfonyl imide), for which an experimental investigation of the ESW is available. Thermodynamic oxidation and reduction potentials have here been estimated by different models based on calculations either on single ions or on ionic couples. Various Density Functional Theory (DFT) functionals (MP2, B3LYP, B3LYP including a polarizable medium and empirical dispersion forces) were exploited. Both vertical and adiabatic transitions between the starting states and the oxidized or reduced states were considered. The approach based on calculations on ionic couples is not able to reproduce the experimental data, whatever the used DFT functional. The best quantitative agreement is obtained by calculations on single ions when the MP2 functional in vacuum is considered and the transitions between differently charged states are vertical (purely electronic without the relaxation of the structure). The B3LYP functional underestimates the ESW. The inclusion of a polar medium excessively widens the ESW, while a large shrinkage of the ESW is obtained by adopting an adiabatic transition scheme instead of a vertical transition one.