Richard M.
Fogarty
a,
Richard P.
Matthews
ab,
Patricia A.
Hunt
*ac and
Kevin R. J.
Lovelock
*d
aDepartment of Chemistry, Imperial College London, UK
bSchool of Health, Sport and Bioscience, University of East London, UK
cSchool of Chemical and Physical Sciences, Victoria University of Wellington, New Zealand. E-mail: patricia.hunt@vuw.ac.nz
dDepartment of Chemistry, University of Reading, UK. E-mail: k.r.j.lovelock@reading.ac.uk
First published on 21st March 2025
The electronic structure of ionic liquids (ILs) is a key factor in their chemical reactivity. Experimental techniques provide insight into IL electronic structure (e.g., X-ray photoelectron spectroscopy, XPS), but are impractical for screening large numbers of potential ILs. Computational screening offers an alternative approach, but current ab initio calculation methods (ion-pairs or large calculations with periodic boundaries) are not suitable for screening. We establish that a simple and computationally low-cost method, lone-ions evaluated at the B3LYP-D3(BJ)/6-311+G(d,p) level employing a generalised solvation model SMD (solvation model based on density), captures IL liquid-phase density-of-states (DoS) with good accuracy by validating against XPS data for a wide range of ILs. The additivity of the results from individual lone-ion calculations provides a significant advantage, enabling predictions of the DoS for a large number of ILs and delivering a significant step towards the computational screening of ILs for many applications.
The possibility, out of the potentially vast number of ILs, that an ideal IL exists for a particular application is an appealing prospect. The challenge is to identify the optimal IL. Large-scale experimental screening is daunting, making computational screening hugely advantageous. The use of expensive and technically demanding ab initio methods or DFT molecular dynamics to provide reliable computational results is well established for ILs.24 Computational methodologies for screening significant quantities of ILs must avoid such expensive, technically demanding calculations, but still capture sufficient information, in particular the solvation environment of the cations and anions, to be representative of the IL. Therefore, critical in the context of large-scale screening is the minimum level of computation required to reliably deliver information on liquid-phase ions.
To have high confidence in the predictions produced by any computational method, experimental validation is essential. The most established validation procedure for the electronic structure of ILs is a visual match of the calculated IL DoS with the experimental valence electronic structure measured using non-resonant valence X-ray photoelectron spectroscopy (XPS, e.g., hν = 1486 eV) and/or ultraviolet photoelectron spectroscopy (UPS, e.g., hν = 21.2 eV).25–43 Gelius-weighted IL DoS, where photoionisation cross-sections are included,44 can give more confidence in the computed spectra.33–35,37–41 However, the reliance on a purely visual validation of computed spectra with experimental spectra is unsatisfactory and will not work for many ILs, given that anion contributions dominate the experimental photoelectron spectra, meaning that cation contributions cannot be observed and validated.
A more quantitative validation of DoS calculations is now possible using recently reported experimental XPS data for 60 ILs, which includes values for Ei(IL) and the energy difference between the cation and anion highest occupied orbitals within the IL environment, ΔEB(ion HOFO) = EB(cation HOFO) – EB(anion HOFO) (where EB = binding energy).21,22 In particular the comparison of calculated and experimental ΔEB(HOFO) values allows us to validate the ability of the calculated DoS to the capture liquid-phase ion–ion solvation environment of both the cation and anion.
There is currently no established method of calculating the IL DoS that is inexpensive and technically undemanding. Calculations of ion-pairs in the gas-phase (ion-pair-GP), which are relatively inexpensive but reasonably technically demanding, perform acceptably, although small EB shifts are required to obtain a good visual match of the IL DoS with experimental valence XP spectra.30,32–35,37 Summation of individual lone-ions-GP DoS to produce an IL DoS are inexpensive and technically undemanding, but do not adequately capture the bulk IL DoS.25–29,32,37
IL DoS can potentially be captured within a liquid environment using a generalised solvation model SMD (solvation model based on density),45,46 which involves placing the substrate/solute in a cavity and representing the surrounding liquid as a charge distribution on the cavity surface. The SMD model can account for electrostatic effects on electronic structure, but not specific covalent interactions such as ion–ion coordination within the first solvation sphere. Recently, anion–cation interactions in ILs were established to be primarily electrostatic and non-specific using a combination of EB(core) from XPS and very expensive and technically demanding ab initio molecular dynamics calculations.47 These results indicate that the SMD might work well for IL DoS.
In the current manuscript we report on an efficient and accurate computational method for evaluating IL DoS. The methodology is validated against a range of experimental data and more complex calculations. Two new quantitative (non-visual) validation methods are presented, using recently published Ei(IL) and ΔEB(ion HOFO) experimental data for 39 ILs21,22 to establish the accuracy of the computed IL DoS (Fig. 1c and d). The validated method is then applied to a wide range of IL systems, demonstrating the ability to screen and provide predictive information for 560 ILs (Fig. 2). The ions studied cover a wide chemical space, including six imidazolium cations, a phosphonium cation, two ammonium cations, and 35 anions; see ESI,† Section S1 for a full list including chemical structures.
![]() | ||
Fig. 1 (a) Comparison of the experimental non-resonant valence XP spectrum (recorded at hν = 161.0 eV) with calculated Gelius-weighted DoS (calculated at hν = 160 eV) for varying model systems for [C4C1Im][SCN], vertically offset for clarity. The yellow and pink filled rectangles are to guide the eye to the EB regions where anionic and cationic contributions respectively dominated the experimental valence XP spectrum. The vertical dashed grey lines represent the calculated cation HOFO energies and anion HOFO energies for the different systems, demonstrating ΔEB(ion HOFO) for each system. The experimental valence XP spectrum was charge corrected as given in ref. 22. (b) Comparison of calculated Gelius-weighted DoS (calculated at hν = 160 eV) for two model systems (dimer-GP and lone-ions-SMD) for [C4C1Im][SCN]. A single shift was applied to each of the calculated Gelius-weighted DoS; the EB shift value for each calculated Gelius-weighted DoS is given on the extreme right-hand side of the figure. Correlations of experimental (taken from ref. 22) and calculated descriptors of IL electronic structure: (c) ΔEB(ion HOFO) and (d) ionisation energy, Ei(IL). Calculations are for lone-ions-SMD (39 ILs in total). |
![]() | ||
Fig. 2 Predictions using lone-ions-SMD. (a) IL ionisation energies, Ei(IL,calc.) for lone-ions-SMD. (b) ΔEB(ion HOFO,calc.) = EB(cation HOFO,calc.) − EB(anion HOFO,calc.) for lone-ions-SMD. When ΔEB(ion HOFO,calc.) is positive the anion is the IL HOMO (red), when ΔEB(ion HOFO,calc.) is negative the cation is the IL HOMO (blue) and when ΔEB(ion HOFO,calc.) is zero either the cation and the anion could be the IL HOMO (white). Ei(IL,calc.) for lone-ions-SMD used to produce these graphs are given in ESI,† Tables S12 and S13. |
Optimisation was carried out without symmetry constraints. The SCF convergence criteria were tightened from the Gaussian 09 defaults to 10−9 on the density matrix and 10−7 on the energy matrix (scf = conver = 9). The numerical integration grid was a pruned grid with 99 radial shells and 590 angular points per shell (int = ultrafine). Frequency analysis was carried out for all optimised structures, which are confirmed as minima by the absence of imaginary modes.
The SMD (solvation model based on density) was used.46 Unless otherwise stated, [C4C1Im][SCN] SMD parameters were used for [C4C1Im][SCN] IL calculations while [C4C1Im][PF6] parameters were used for calculations for other ILs (very slightly different SMD parameters were used for the [Bi2Cl8]2− anion as detailed in ref. 61 and 64). Some testing of the relative permittivity (εr) was carried out using default parameters within G09 and G16. The SMD parameters employed are shown in Table S3 (ESI†).
Gelius-weighted IL DoS were computed (lone-ions-GP, lone-ions-SMD, ion-pair-GP, ion-pair-SMD, dimer-GP and dimer-SMD, where dimer = two ion-pairs); for lone-ions-SMD the IL DoS was produced via summation of individual cation and anion lone-ion-SMD DoS. A single EB shift is applied to align the computed HOMO (lowest EB, i.e., right most peak in Fig. 1a and b) of the DoS with the lowest EB XPS peak.
Further computational details can be found in the ESI,† Section S2.
Further validation of the ability of the lone-ions-SMD methodology to recover the liquid-phase inter-ion solvation is evidenced by excellent visual matches of calculated DoS and three different experimental XPS data types. First, good visual matches are obtained at different hν for lone-ions-SMD DoS and experimental valence XPS for both [C4C1Im][SCN] and [C8C1Im]Cl (ESI,† Fig. S2 and S3); the relative peak areas of XP spectra vary depending on the wavelength of the light, and the use of Gelius-weighting allows this variation to be recovered by the lone-ions-SMD DoS calculations. Second, computed lone-ion-SMD DoS show superb visual matches to experimental lone-ion valence XPS data for 12 anions and four cations (ESI,† Section S4), validating the calculations; for example, all four cyano-containing anions give excellent visual matches to the calculated lone-ion DoS to experimental data (ESI,† Fig. S5). Third, computed lone-ion-SMD partial DoS for one atomic orbital show good visual matches to experimental valence resonant XPS data for six anions and five cations (ESI,† Section S5), again validating the calculations; for example, for [HSO4]− the O 2p lone-ion-SMD partial DoS visually matches very well to the experimental oxygen valence resonant XPS data (ESI,† Fig. S8f).
Quantitative comparisons using two different parameters/metrics show that the lone-ions-SMD method recovers the liquid-phase inter-ion solvation very well. First, the quantitative comparison of ΔEB(ion HOFO,exp.)22versus ΔEB(ion HOFO,lone-ions-SMD) gives an excellent linear correlation (R2 = 0.91, Fig. 1c). For example, for [C4C1Im][SCN] lone-ions-SMD gave ΔEB(ion HOFO,calc.) = 2.2 eV, which is the same (within error limits) as ΔEB(ion HOFO,exp.) = 1.9 eV ± 0.5 eV.21 In contrast, ΔEB(ion HOFO,lone-ions-GP) gives a poor linear correlation (R2 = 0.36, ESI,† Fig. S10a).
In the second quantitative assessment, comparison of Ei(exp.) taken from ref. 22 with Ei(calc.) (i.e., the calculated HOMO EB) for lone-ions-SMD gives a good linear correlation with R2 = 0.89 (Fig. 1d). In contrast, Ei(lone-ions-GP) gives a poor correlation with R2 = 0.50 (ESI,† Fig. S11a). These results also show that lone-ions-SMD calculations are sufficient to capture the relative difference in Ei between different ions. The y-intercept = −2.2 eV of Fig. 1d does not go through zero demonstrating an approximately constant error relative to Ei(exp.). For [C4C1Im][SCN], Ei(dimer-SMD), which is expected to give accurate results, is 1.63 eV lower than Ei(exp.). This result suggests that both the approximately constant offset between Ei(calc.) and Ei(exp.) and the need for a single EB shift to align computed HOMO (top of the DoS) with the lowest EB XPS peak is due to the DFT method, rather than the SMD method not capturing the liquid-phase solvation environment.65–67 Our results are consistent with previous studies where strong linear correlations (but with non-zero intercepts) were observed between DFT and experimental Ei for a wide range of functionals and molecules.68,69
Our calculations appear to have picked up a flaw in the experimental analysis for one IL, demonstrating that calculations can provide a valuable check for experimental results. An outlier for the correlation in ΔEB(ion HOFO) is the IL [C2C1Im][FAP] ([C2C1Im]+ = 1-ethyl-3-methylimidazolium and [FAP]− = tris(pentafluoroethyl)trifluorophosphate), ΔEB(ion HOFO,lone-ions-SMD) = −0.70 eV compared to ΔEB(ion HOFO,exp.) = −3.8 eV (Fig. 1c). On close inspection of the experimental results published in ref. 22 we expect that there was a small anion contribution at a similar EB to the cation contribution, which was very challenging to observe experimentally due to the overlap of the cationic and anionic experimental contributions, leading to an incorrect assignment of the anion EB(HOFO).
Lone-ions-SMD captures the liquid-phase solvation for dianion-containing ILs very well. Two of the ILs studied contained dianions, [P6,6,6,14]2[ZnCl4] ([P6,6,6,14]+ = tetradecyl(trihexyl)phosphonium) and [C8C1Im]2[ZnCl4]. These results provide evidence that no special treatment (beyond lone-ions-SMD) is required for ILs containing dianions.
A sum of non-specific electrostatic inter-ion interactions, as present in lone-ions-SMD, have the same effect on IL orbital energies as found in an explicit larger cluster, i.e., a dimer. Knowledge of the relative impact of generalised, non-specific solvent interactions, versus specific (covalent or ionic) ion-pair interactions is highly desirable. For [C4C1Im][SCN] the computed DoS for the dimer-GP and sum of lone-ions-SMD are visually near identical (Fig. 1b). Moreover, the individual ion Ei(HOFO) and in the IL Ei(HOMO) shift required to visually match the experimental XP spectrum is near identical, +1.95 eV and +1.94 eV respectively.
For the most important input parameter in the SMD model, the relative permittivity εr, ∼11 is a good general value for predictions of screening, with an Ei(HOFO) error of < ± 0.2 eV. For ILs 7.5 < εr < 20 is the normal range of experimental εr values.70–73 A single value for εr that can reasonably represent any IL is highly desirable to simplify the survey procedure. Within the SMD model the dependence of MO energies on εr has been evaluated (ESI,† Section S7 and Fig. S12). The impact on ΔEB(ion HOFO) is shown to scale as (εr − 1)/εr (ESI,† Fig. S13). A simple estimate for the error incurred using a generalised representative εr parameter can be determined by comparing the change in lone-ions-SMD Ei(cation HOFO) over a range of εr. The change in Ei(cation HOFO) lone-ions-SMD for εr = 7.5 to 20, for a dialkylimidazolium cation is 0.4 eV (ESI,† Fig. S13), with the average Ei(cation HOFO, lone-ions-SMD) occurring at εr ∼ 11.
Based on the results presented thus far, we have established that an excellent compromise between accuracy and computational effort is achieved using lone-ions-SMD calculations. Lone-ions-SMD calculations are the simplest and cheapest model (studied here) that accurately captures the IL DoS. We have shown that lone-ions-SMD matches well to experimental data both qualitatively (visual matches to experimental valence XP spectra and RXP spectra) and quantitatively (linear correlations for Ei(IL) and ΔEB(ion HOFO)). Using lone-ions-SMD compared to e.g., five ion-pairs (where the effect of the counterions would be close to fully captured) causes minimal loss of output quality with respect to the valence levels, and avoids the necessity of dynamically sampling a wide range of cluster conformers.
Relative to the lone-ions-SMD, calculations involving ion-pairs and dimers carry a significantly greater computational cost and require more skill to generate and evaluate, and thus are too expensive and time-consuming for screening many ILs. An attempt has been made to quantify the computational cost saving for different model systems (ESI,† Table S4). As more ions are added to a system, a larger number of possible conformers (relative arrangement of the cations and anions) needs to be evaluated. Broadly, the time saved increases as more ILs are screened. For example, to screen 1000 ILs made up from a pool of 20 cations and 50 anions, the lone-ions calculations have ∼430 times lower computational cost than ion-pair calculations and ∼7900 times lower computational cost than dimer calculations (ESI,† Table S4). These estimates do not include the user time required to select the ion-pair/dimer conformers or run a dynamic sampling procedure, which is very dependent on the experience of the user.
The additivity of the individual cation and anion lone-ion calculation results provides a significant advantage when screening IL combinations, and means low-cost predictions of DoS can be achieved for a large number of ILs. We have made predictions based on 16 cations and 35 anions, which can form 560 different ILs; Ei(calc.) (Fig. 2a) and ΔEB(ion HOFO,calc.) (Fig. 2b).
Using the validated lone-ions-SMD method, we predict that: (i) the most readily ionised group of ILs is [cation]2[FeCl4] with Ei(calc.) = 4.9 eV (Fig. 2a, predicted IL DoS ESI,† Fig. S15a); (ii) the least readily ionised IL is [N2,2,1,0][SbF6] with Ei(calc.) = 10.1 eV (Fig. 2a, predicted IL DoS ESI,† Fig. S15c). This second result is consistent with the use of ammonium cations and highly fluorinated anions as supporting electrolytes in electrochemistry.74,75
ΔEB(ion HOFO) (Fig. 2b) can be used to identify the ion HOFO that is the IL HOMO (ESI,† Section S8). For ILs it might reasonably be assumed that the anion provides the IL HOMO (i.e., the anion will be ionised before the cation). However, this assumption is not always true. Using the data collected here, the anion is predicted to be the HOMO for 474 ILs and the cation is predicted to be the HOMO for 86 ILs (ESI,† Fig. S14). Thus, 15% of the ILs studied here are predicted to have a cation-based IL HOMO. This conclusion might seem counterintuitive; for ∼15% of ILs the process [cation]+ → [cation]2+ + e− is more favourable than the process [anion]− → anion + e−. For example, for [C4C1Im][SbF6] ΔEB(ion HOFO,calc.) = −2.8 eV (Fig. 2b) is a remarkably large negative value (predicted IL DoS ESI,† Fig. S15b). We have previously shown that atoms of a high electronegativity (within an anion) increase many of the MO EB, including the HOMO.22 Anions containing F show a particularly large effect,22,76 as clearly evidenced in Fig. 2b. The one anion that is not halide-containing (and particularly has no F) that gives a large EB(anion HOFO) is [B(CN)4]− (tetracyanoborate).
For ILs with ΔEB(ion HOFO,calc.) ≈ 0 eV, which are represented by the white/very pale colored boxes on Fig. 2b (e.g., [CnC1Im][NTf2] where n = 2 to 8 and [NTf2]− = bis[(trifluoromethane)sulfonyl]imide), it is possible that environmental effects could change the relative ordering of the IL ion HOFOs and ultimately affect the IL properties and reactivity. Such ILs could be expected to be more strongly influenced by the presence of solutes such as ion contaminants or water, H-bonding additives or changing IL organisation at an electrode surface.
It is expected that IL DoS will be strongly linked to reactivity. Consequently, further investigation is needed to determine whether this lone-ion-SMD calculation method, and therefore IL DoS, can be used to predict other IL properties and reactivity.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d5cp00214a |
This journal is © the Owner Societies 2025 |