Open Access Article
Imgon Hwang
*ab,
Mantas Leketasab,
Ryan Bragg
c,
Sittipong Kaewmorakotab,
John M. Griffin
c and
Robert A. W. Dryfe*ab
aDepartment of Chemistry, University of Manchester, Oxford Road, Manchester, M13 9PL, UK. E-mail: imgon.hwang@manchester.ac.uk; robert.dryfe@manchester.ac.uk
bHenry Royce Institute, University of Manchester, Oxford Road, Manchester, M13 9PL, UK
cDepartment of Chemistry, Lancaster University, Lancaster, LA1 4YB, UK
First published on 18th February 2026
Electric double-layer capacitors (EDLCs) are promising energy storage devices due to their high power density and long cycle life, but their lower energy density compared to batteries remains a major limitation. Highly concentrated electrolytes present a safer and more sustainable alternative to conventional organic electrolytes, while retaining the ability to operate at high-voltage. Since many water-in-salt (WIS) electrolytes are based on fluorinated anions, the length of the fluorinated chain plays a critical role in determining electrolyte properties. In this study, we investigate the impact of anion structure by comparing three fluorinated imide-based salts, i.e. LiTFSI (lithium bis(trifluoromethanesulfonyl)imide), LiPFSI (lithium bis(pentafluoroethanesulfonyl)imide), and LiNFSI (lithium bis(nonafluorobutanesulfonyl)imide), for supercapacitor applications. Electrochemical characterization, Raman spectroscopy, and solid-state nuclear magnetic resonance (NMR) reveal that LiPFSI, with its higher proportion of double-donor single-acceptor (DDA) hydrogen bonding profiles, forms a highly coordinated anionic network at an optimal concentration of 5 mol kg−1. Solid-state NMR spectroscopy shows that LiPFSI preferentially occupies the in-pore environment of nanoporous carbon, achieving a favourable ion to solvent balance that supports efficient charge storage and enhanced capacitance. In contrast, although LiNFSI can form dense ion agglomerates at high concentrations, the bulky anion and associated higher viscosity restrict ion transport and pore accessibility. These findings highlight the importance of anion selection and concentration, and pore-level electrolyte distribution in tailoring high-concentrated electrolytes for next-generation high-voltage, high-capacity supercapacitors.
An important recent development with aqueous electrolytes has been the development of highly concentrated water-in-salt (WIS) electrolytes, which significantly expanded the electrochemical stability window (ESW) towards 3 V, i.e. far exceeding the 1.23 V limit of conventional aqueous electrolytes.11 This innovation, although reported in the context of Li+-ion batteries, has been used to narrow the performance gap between aqueous and organic-based supercapacitors, making WIS electrolytes attractive for next-generation energy storage applications.11–23 Among various electrolyte candidates, imide-based anions such as lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) have demonstrated excellent thermal and electrochemical stability, in turn prompting research into their use in high-voltage aqueous supercapacitors.11 These unique properties have driven extensive exploration of imide-based WIS electrolytes for high-voltage energy storage applications.23–31 Despite the success of WIS electrolytes, most prior studies have focused on single salt systems at extremely high concentrations (typically >20 mol kg−1, hereafter denoted m),31 where water molecules are fully coordinated and free solvent is nearly absent. However, recent studies have revealed that many relevant physicochemical changes—such as ion pairing, clustering, and non-ideal transport behavior—begin to manifest themselves at intermediate concentrations, well below the canonical WIS threshold.32,33 For example, Suh et al. demonstrated that LiTFSI and NaClO4 electrolytes reach maximum ionic conductivity near 5–6 m, with further increases in concentration leading to enhanced ion aggregation and reduced mobility due to stronger cation–anion interaction.34 These findings suggest that the optimal properties of concentrated electrolytes may emerge in the 5 m regime, where conductivity is maximized, marking a transition where ion–ion interactions begin to dominate over ion–solvent interactions, even though water is still the majority component.34
This study introduces a comparative analysis of three fluorinated imide-based anions, LiTFSI, lithium bis(pentafluoroethanesulfonyl)imide (LiPFSI), and lithium bis(nonafluorobutanesulfonyl)imide (LiNFSI), to elucidate how variations in fluorinated chain length influence ion solvation, electrolyte–carbon interactions, and charge storage mechanisms. A combination of Raman spectroscopy and solid-state nuclear magnetic resonance (NMR) spectroscopy is employed to provide multi-scale understanding. Raman spectroscopy is used to investigate the solvation structure and hydrogen-bonding interactions in the electrolyte, providing insights into how different anions influence the local ionic environment,35,36 while solid-state NMR offers molecular-level insights into electrolyte behavior within nanoporous carbon electrodes, clarifying anion-dependent ion accessibility and charge storage mechanisms.37,38
To the best of our knowledge, this is the first comparative study of LiTFSI, LiPFSI, and LiNFSI electrolytes at intermediate concentrations (∼5 M), combining the methods mentioned above to link bulk solvation structure and pore-level ion partitioning with supercapacitor performance. Unlike prior work that predominantly focused on single salts at ultra-high concentrations, this work probes the structure–property relationships across different anions, integrating analysis from bulk electrolyte characterization to in-pore molecular behavior and device-level electrochemical performance. This approach provides insight into how molecular–scale interactions translate into ion transport, electrolyte partitioning, and charge storage efficiency, thereby informing the rational design of high energy density aqueous supercapacitors.
:
1
:
1 mass ratio were dispersed in iso-propyl alcohol (IPA, Fisher Chemical, UK). The slurry was stirred for at least 2 hours until most of the solvent evaporated to achieve a non-sticky, mouldable consistency. The flattened film was rolled into 50 µm thick sheets using an Electric Cold Roller Press (model MSK-HRP-MR100DC MTI, KJ Group, USA). The sheets were dried in an oven at 70 °C for one hour, then transferred to a vacuum oven and dried overnight at 100 °C. Dried sheets were punched into discs of 12 mm diameter (approximately 7 mg). Electrodes using YP-50F carbon (Kuraray, Japan) were prepared following the same procedure.
The electrodes were prepared manually, and some variation in mass loading between samples is therefore unavoidable. Based on repeated electrode mass measurements, the mass loading variation was estimated to be within approximately ±20%, which is typical for manually prepared porous carbon electrodes. To minimize the influence of this variability, electrochemical measurements were consistently performed using material taken from the central region of the electrode coatings, where thickness and composition are most uniform.
Raman spectra were acquired using a Renishaw inVia microscope equipped with a 532 nm excitation laser (Renishaw, UK), operating at a power of 0.274 mW, with a grating of 1800 lines per mm and a 50× long working distance objective. Deconvolution of the Raman spectra was carried out using a Gaussian function, allowing peak positions to vary freely to achieve the best fit with the experimental data. Deconvolution was performed with Microcal Origin 2022b (OriginLab Co., Northampton, MA, USA), applying 5-point Savitzky–Golay smoothing to determine the number and positions of sub-bands. Prior to analysis, the spectra were baseline-corrected. The bands were then deconvoluted into the most appropriate number of Gaussian components by monitoring both the residual sum of squares and the χ2 value as functions of the number of Gaussian components.39 X-ray diffraction (XRD) analysis was performed on free-standing activated carbon samples using a Philips X'pert PRO diffractometer with Cu Kα radiation (λ = 0.154 nm) and operating at 40 kV and 30 mA.
The wetting properties of the concentrated electrolyte were investigated using a Theta Optical Tensiometer (Biolin Scientific, Finland) with OneAttension software (version 2.3). In this experiment, the droplet size was controlled at a fixed volume of 1.5 µL. Measurements were performed statically in a humidity chamber to minimize liquid evaporation. A non-porous carbon substrate, the basal plane of highly ordered pyrolytic graphite (HOPG, ZYA grade, Scanwel, USA), was used as a proxy for the activated carbon material. Images were captured after the electrolyte was placed on the HOPG (t = 30 s) at a rate of 10.2 Hz. The contact angles were analyzed using OneAttension software, based on the Young–Laplace equation.40,41 Contact angle measurements were performed with more than ten independent droplets for each electrolyte.
The 3-electrode electrochemical investigations were performed in Poly(etheretherketone) (PEEK) Swagelok-type cells by sandwiching a porous glass microfiber membrane (Whatman, thickness = 0.20 mm) between two electrodes and graphite foil current collectors. For the electrolyte stability investigations, three-electrode cells were built with a working AC electrode (∼7 mg), a greater mass of AC (∼30 mg) as a counter electrode and a silver-wire pseudo-reference electrode (AgQRE37,42). The AgQRE was calibrated using a standard redox couple (10 mM ferrocene solution in an ionic liquid), and its stability was verified subsequently by monitoring the open-circuit potential directly in the corresponding electrolytes over the duration of the electrochemical measurements. All the cells were filled with the 1, 5 or 10 m LiTFSI, LiPFSI, or LiNFSI electrolyte solutions. The electrochemical measurements in three-electrode cells were conducted by cyclic voltammetry (CV) at 1, 10, and 100 mV s−1 and the galvanostatic charge and discharge (GCD) method at 1, 2, 5, 10 A g−1 with the working electrode swept between 0 and +0.5 V (positively) and 0 and −0.5 V (negatively) vs. the Ag wire quasi-reference electrode. Electrochemical impedance spectroscopy (EIS) measurements were performed at 0 V (two-electrode configuration) and at 0 V (vs. AgQRE, three-electrode configuration), using a 5 mV amplitude ac modulation across a frequency range of 100 kHz to 10 mHz.
All electrochemical measurements (CV, GCD, and EIS) were performed using independently prepared electrodes and electrolyte batches, showing good reproducibility (variation ±5%). Electrochemical performance is reported as gravimetric capacitance (F g−1) to account for any minor variations in electrode mass.
As well as bulk interactions, understanding of capacitance requires some knowledge of electrolyte–substrate interactions, hence the wettability of the electrolytes on highly ordered pyrolytic graphite (HOPG) was evaluated using contact angle measurements shown in Fig. S3. At a concentration of 5 m for each electrolyte, the solution contact angles in air were approximately 23(±2.8)° for LiTFSI, 24(±3.6)° for LiPFSI, and 18(±5.3)° for LiNFSI. These values are significantly lower than those of pure water on HOPG (62°),45 reflecting the affinity of the fluorinated anion for the graphite surface. LiNFSI shows the lowest contact angle: increased wettability is significant because it reflects stronger adsorption of the electrolyte on the electrode, which can increase ion accessibility and improve electrochemical properties at the electrode–electrolyte interface.
To achieve deeper insights into how the fluorinated anions affect the conductivity and electrochemical behavior of the electrolytes, we employed Raman spectroscopy to probe the chemical environment of water molecules in the electrolytes. Fig. 1 shows the OH stretching modes of LiTFSI, LiPFSI, and LiNFSI aqueous solutions at concentrations of 1 m, 5 m, and 10 m. The broad Raman band was deconvoluted into three Gaussian peaks corresponding to double-donor double-acceptor (DDAA), double-donor single-acceptor (DDA), and single-donor single-acceptor (DA) hydrogen bonding motifs. The deconvolution procedure and representative fitting results for the OH stretching region are provided in Fig. S4 and Table S1. This approach is commonly used to analyze changes in hydrogen-bonding networks that arise from variations in electrolyte composition and concentration.46,47
![]() | ||
| Fig. 1 Raman OH stretching vibration of LiTFSI, LiPFSI, and LiNFSI at the (a) 1 m, (b) 5 m, and (c) 10 m concentrations. | ||
At a concentration of 1 m, the Raman spectra of all three electrolytes exhibit similar OH vibration characteristics that closely resemble those of pure water (Fig. 1a). This similarity suggests that at low concentrations, the solvation structures of the electrolytes are largely comparable, with lithium cations (Li+) being well-solvated by water molecules and the anions (TFSI−, PFSI−, NFSI−) remaining largely dissociated in solutions. As the concentration increases to 5 m, distinct differences in solvation and ion association begin to emerge. Deconvolution of the OH stretching region reveals that the proportion of strongly hydrogen-bonded water molecules (DDAA structures) varies among the electrolytes in the order LiNFSI > LiTFSI > LiPFSI, while the content of moderately hydrogen-bonded species (DDA structures) follows the reverse trend (Fig. 1b). These variations reflect how each anion modulates the local solvation environment, indicating that the identity of the anion plays a critical role in determining the balance between free and associated Li+ species. This evolving solvation behaviour suggests that, at intermediate concentrations, anion-dependent ion association plays a dominant role in shaping electrolyte structure. In the case of LiPFSI, the lower DDAA and higher DDA content imply the formation of a more compact and aggregated anionic network,48 which likely limits the availability of free, fully solvated Li+ ions.49 Conversely, the higher DDAA content observed for LiNFSI indicates the presence of a larger fraction of bulk-like water molecules that are less directly involved in Li+ solvation. This behavior is consistent with stronger ion–ion interactions and the formation of contact ion pairs or aggregates, which reduce the availability of water molecules for cation coordination. LiTFSI exhibits intermediate characteristics, reflecting a balance between ion dissociation and aggregation arising from the competing effect of charge delocalization and coordination flexibility of the TFSI− anion. Together, these distinctions highlight a progressive evolution of solvation with concentration, driven by enhanced ion–ion interactions and reduced effective solvation of Li+.
Additional insight into these interactions is provided by Raman spectra of LiNFSI in aqueous solution, which reveal a concentration-independent S–N–S stretching mode centred at approximately 732 cm−1 across the full investigated concentration range (1–10 m, see Fig. S5). This spectral invariance suggests a stable solvation structure for the NFSI− anion and indicates that increasing salt concentration does not significantly perturb its local coordination environment. Notably, a distinct shoulder at ∼752 cm−1 persists across all concentrations, including the most dilute regimes (normalized Raman spectra as shown in Fig. S5d). The persistence of this higher-frequency feature under dilute conditions rules out its origin being due to bulk aggregation or concentration-induced clustering. Instead, its spectral position and stability suggest the presence of contact ion pairs (CIP) or chelated Li+–NFSI− complexes, coexisting alongside more weakly bound or solvent-separated species.50,51 Similar shoulders have been identified in other aqueous salt systems, where they have been conclusively assigned to strongly coordinated contact ion pairs, observable even at low concentrations.50,51 This behaviour stands in marked contrast to that observed for the structurally related LiTFSI and LiPFSI salts. In those systems, the S–N–S stretching mode typically appears near ∼740–745 cm−1, but exhibits a progressive blue shift (toward higher frequencies) as salt concentration increases. This shift has been widely interpreted as a spectroscopic signature of increasing Li+–anion interaction strength, caused by the formation of contact ion pairs and higher-order aggregates (e.g., triple ions or clusters) due to reduced solvent availability and enhanced electrostatic correlation at higher concentrations.52 Namely, in contrast to LiTFSI and LiPFSI, where increased concentration induces stronger pairing, LiNFSI appears to form stable CIPs even in dilute solutions, contributing to persistent transport limitations across the full concentration range. Notably, the presence of stable contact ion pairs in LiNFSI even at low concentration, as suggested by Raman analysis, provides a potential microscopic explanation for its higher viscosity and lower ionic conductivity. Persistent ion pairing and aggregation increase the effective hydrodynamic size of charge carriers and reduce the population of free Li+ ions, thereby simultaneously enhancing viscous resistance and suppressing ionic mobility, consistent with molecular dynamics and scattering studies of imide-based water-in-salt electrolytes.48
These differences in solvation structures across concentrations highlight distinct trends in ion pairing, aggregation, and water interactions, all of which are crucial determinants of electrolyte performance in energy storage applications.48 The unique behavior of LiPFSI, characterized by a higher DDA content at 5 m, is indicative of a moderately structured hydrogen-bonding network that promotes enhanced ion transport and improved pore accessibility (vide infra). Conversely, LiNFSI, with its larger anion size and comparatively weaker hydrogen-bonding networks, exhibits more pronounced transport limitations at higher concentrations due to its solvation dynamics and ion aggregation behavior.
To experimentally assess the electrochemical performance of high-concentration electrolytes containing different anionic species, symmetric two-electrode coin cells were assembled and evaluated using cyclic voltammetry (CV) and galvanostatic charge–discharge (GCD) techniques. The electrodes were fabricated from commercially available activated carbon (YEC-8A), featuring a surface area of 2390 m2 g−1 and a pore size distribution with a peak maximum at approximately 0.87 nm.37 As shown in Fig. S6, the CVs obtained at a scan rate of 10 mV s−1 reveal that all three electrolytes achieve maximum capacitance near 5 m concentration: the capacitance values are presented in Fig. S7a. At higher scan rates, the current response is increasingly dominated by surface-controlled capacitance, thereby minimizing differences in overall capacitance due to reduced diffusion effects.53,54 Under these conditions, the capacitance largely reflects the electrolyte's ability to support rapid ion transport—a behaviour consistent with the observed ion conductivity and viscosity trends (Fig. S1 and S2). The CV responses shown in Fig. 2a–c, recorded at a scan rate of 1 mV s−1, reveal that pore transport effects become increasingly important, leading to more pronounced differences in capacitance across the electrolytes (Fig. S7b).54 Under these slow conditions, sufficient time is available for ions to penetrate the nanoporous carbon electrode, making diffusion-controlled capacitance within the pores the dominant charge storage mechanism.37 At this scan rate, LiTFSI and LiPFSI exhibit comparable capacitance values (Fig. S7b). Similarly, contact angle measurements on HOPG suggest nearly identical wettability between these electrolytes (see discussion above and Fig. S3). In contrast, LiNFSI displays significantly lower capacitance at 1 mV s−1. Interestingly, it also shows the lowest contact angle (18(±5.3)° at 5 m), indicating high apparent surface wettability on HOPG: this reflects the trade-off between the higher adsorption of the NFSI and reduced pore accessibility because of its larger size/lower viscosity. Note that the differences between the idealized (planar, basal plane HOPG) and capacitive electrode (activated carbon, mixed with binder and graphene) materials makes direct comparison difficult. We emphasis the notable reduction in contact angle relative to pure water (62°)45 – highlighting the substantial enhancement in wetting behavior induced by the presence of these salts. This observation is particularly important because good wetting is a prerequisite for electrolyte access to nanoscale pores. This behavior is consistent with porous electrode theories, which suggests that during longer charge times, ions first occupy larger interparticle voids before diffusing into smaller pores.55
Capacitance values derived from GCD measurements (Fig. 2d–f and S8–S10) align well with the CV results, reaffirming that peak capacitance is achieved at a concentration of 5 m. This consistency across techniques suggests that the high bulk ionic conductivity observed at 5 m (Fig. S1b) effectively translates to efficient in-pore ion transport, supporting enhanced charge storage. In other words, the favourable bulk transport properties at this concentration appear to facilitate ion access into the porous electrode structure, enabling maximum utilization of the available surface area. This correlation highlights the importance of optimizing not just ionic concentration, but also the interplay between ion mobility and pore accessibility. At 1 m and 5 m, capacitance increases with the applied voltage window, potentially due to enhanced charge storage and minor contributions from faradaic side reactions, such as electrode oxidation. However, at 10 m, this trend plateaus despite the higher ionic concentration, likely due to increased viscosity and reduced mobility, which hinder further charge accumulation. Although both CV and GCD characterize charge storage, they do so under different control modes. CV, being a potentiostatic technique, allows rapid evaluation of surface and diffusion-controlled processes across a range of time scales, depending on the scan rate. In contrast, GCD, being galvanostatic in nature, subjects the system to a constant current, revealing how the electrolyte responds to sustained ion flow and longer charge duration. These complementary approaches offer a more comprehensive understanding of the underlying electrochemical mechanisms.56 Finally, while higher voltage windows can improve capacitance, they also increase the risk of electrolyte decomposition or carbon oxidation, especially in more dilute electrolytes, e.g. the 1 m solution. Raman spectroscopy (Fig. 1) indicates that these dilute solutions have lower DDA content and a less structured hydrogen-bonding network, which may increase the availability of free water molecules and facilitate side reactions—ultimately compromising electrochemical stability.
The coulombic efficiency (CE), η = (discharge time/charge time) × 100, calculated from the GCD results (Fig. 2g–i), shows a direct correlation with the CV analysis. At 1 m, the CE is generally low, indicating significant energy losses during charge/discharge cycles (Fig. 2g). This inefficiency is further supported by the asymmetric shape of the discharge curves in both the CV (Fig. S6a) and GCD results (Fig. S8), suggesting the occurrence of parasitic reactions. These findings highlight the limitations of dilute electrolytes in maintaining efficient and reversible charge storage processes. Raman spectroscopy provides additional insight into this behavior. In dilute electrolytes, the weaker hydrogen-bonding networks, which are evidenced by a lower DDA structure, may contribute to increased side reactions and reduced CE. In contrast, at higher concentrations (5 m and 10 m, as shown in Fig. S9 and S10), the stronger solvation structures, with higher DDA content, correlate with enhanced electrochemical stability and improved CE (Fig. 2h and i), consistent with earlier observations.
At a high concentration of 10 m, the electrochemical behavior of LiTFSI and LiPFSI differs from that at lower concentrations. The capacitance values remain nearly constant despite increasing voltage, in line with the ESW observed in the CV data. This stability suggests effective suppression of electrolyte decomposition and electrode degradation, likely due to modified solvation structures that reduce the number of free solvent molecules, as confirmed by Raman spectroscopy (Fig. 1). The high CE observed at 10 m further reflects the minimal side reactions and more efficient charge storage processes under these conditions. In contrast, LiNFSI exhibits markedly different behavior at high concentrations. GCD results reveal a significant drop in charge/discharge efficiency, which is attributed to the high viscosity of the electrolyte, impairing ion transport and increases internal resistance, as is evident from the increased ohmic drop and decreased capacitance (Fig. S10). Additionally, the bulky NFSI− anion, which is non-spherical and conformationally flexible, possesses a larger effective spatial extent than TFSI− and PFSI−. For comparison, structural analyzes of ionic liquids have shown that the effective longer molecular dimension of the TFSI− anion is approximately 0.79 nm,57 a value already comparable to the dominant pore size (∼0.87 nm) of the YEC-8A activated carbon used in this study. The longer perfluoroalkyl chain of NFSI−, combined with strong ion pairing and aggregation at high concentrations, is therefore expected to further increase the effective size of the charge-carrying species, limiting access to narrow micropores and contributing to the observed capacitance loss (see further discussion below). Similarly, X-ray scattering data suggests that the characteristic d-spacing of the aggregates formed by these ions ranges from 0.7 to 1.1 nm, as the anion size increases.48 It is therefore possible that the optimal capacitance of PFSI− anion reflects the match between the modal average pore size and the anion size. The presence of residual water may trigger parasitic faradaic reactions such as water electrolysis. This behavior is consistent with the CV data, which show a sharp increase in capacitance beyond 1.8 V, indicating faradaic contributions rather than purely capacitive charging. These factors collectively limit the practical application of LiNFSI in high-voltage energy storage systems. Long-term cycling stability results (Fig. S11) further reinforce these findings.
To further validate these findings, electrochemical impedance spectroscopy (EIS) was performed. The Nyquist (Fig. S12) and Bode plots (Fig. S13) obtained from EIS are used to provide qualitative insights into the accessibility within the porous electrode structure and elucidate ion transport behavior at the electrode–electrolyte interface as a function of frequency.58 In the Nyquist plot, the high-frequency region reveals that at 1 m concentration (Fig. S12a), LiTFSI exhibits a lower equivalent series resistance (ESR, real-axis intercept) and charge transfer resistance (semicircle diameter) compared to LiPFSI. However, this trend reverses at higher concentrations: at 5 m (Fig. S12b) and 10 m (Fig. S12c), LiPFSI demonstrates lower ESR and charge transfer resistances compared to LiTFSI. These trends are consistent with the CV and GCD results, reinforcing the observed improvements in electrochemical performance at higher concentrations. Notably, LiNFSI consistently shows the highest ESR and charge transfer resistance across all concentrations, reflecting its sluggish ion transport properties. The Bode plot (Fig. S13) provides additional insights into the frequency dependent capacitive behavior of the electrolytes. As the frequency decreases, the phase angle shifts negatively, with ideal supercapacitor behavior approaching −90°.59 At 5 m, LiTFSI reaches a phase angle of −45° more rapidly (at a higher frequency) than LiPFSI, suggesting more efficient ion transport and faster double-layer formation. However, this trend reverses at 10 m, where LiPFSI surpasses LiTFSI, indicating a shift in ion transport dynamics.60 The frequency corresponding to the −45° phase angle serves as a key indicator of capacitive performance, highlighting the balance between surface and pore-level ion transport.59 The observed reversal in the phase angle trend at 10 m suggests that the reduced solvation of LiPFSI facilitates deeper ion penetration into the micropores, enhancing diffusion-controlled charge storage. This observation aligns with previous reports in ionic liquid-based systems,61,62 reinforcing the conclusion that capacitance and overall electrochemical performance are closely governed by the pore characteristics of the electrode and the hydration structure of the WIS electrolyte.
While the two-electrode setup effectively captures overall energy storage behaviour, it is limited in isolating the intrinsic electrochemical characteristics of individual electrodes. To address this, a three-electrode Swagelok-type cell was also employed in this study, allowing for a more accurate and detailed comparison of electrode behaviour. This configuration enables clearer differentiation of the charge storage mechanisms and more precise evaluation of electrolyte-dependent performance differences.
Fig. 3a–c presents the CV results for electrolyte solutions containing different anions, measured across varying concentrations. The CV experiments were conducted at a low scan rate of 1 mV s−1 to accurately detect faradaic reactions such as electrolyte decomposition. For these experiments, the measurements were performed by sweeping the voltage from 0 V (vs. AgQRE) in both positive and negative directions. The results showed that all electrolytes, except for 10 m LiNFSI, exhibited a characteristic rectangular CV shape indicative of EDLC behaviour across all tested concentrations. However, as the scan rate increased (Fig. S13 and S14 at 10 mV s−1 and 100 mV s−1, respectively), the rectangular CV shape gradually distorted. This phenomenon occurs because, at higher scan rates, the ability of ions within the electrolyte to reach the electrode surface becomes increasingly limited at higher currents. In other words, at high scan rates, ions do not have sufficient time to migrate and interact with the electrode surface, leading to the increasing resistance observed in the CV. In a typical EDLC system, the capacitance tends to decrease at high scan rates, as ion transport becomes the rate-limiting factor. However, for 10 m LiNFSI, the opposite trend was observed. At a high scan rate of 100 mV s−1 (Fig. S15), the CV retained a nearly rectangular shape and exhibited the highest capacitance. Conversely, as the scan rate decreased to 1 mV s−1 (Fig. 3c), the CV curve transitioned to almost linear I–V curve, with a sharp decline in capacitance (specific capacitance values calculated from CV areas can be found in Fig. S16). This anomalous behavior can be attributed to the high viscosity and low ionic conductivity of the 10 m LiNFSI electrolyte (Fig. S1). At high scan rates, the dominant charge storage mechanism appears to be limited to interfacial ion rearrangement near the outer surface of the nanoporous electrode. In such cases, ions do not have time to penetrate far into electrode pores, but fast rearrangement of ions at the electrode/electrolyte interface can still support capacitive behaviour. Interestingly, this surface-dominated charge storage remains efficient for 10 m LiNFSI at high scan rates, despite the system's slow ion mobility. At low scan rates, more time is available for ions to access the internal porous network, which ideally should lead to higher capacitance. However, in the case of 10 m LiNFSI, this deeper pore accessibility is severely limited. The large effective size of the NFSI− anion, together with the high viscosity of electrolyte, restricts ion mobility and prevents effective ion transport into sub-nanometre pores. For comparison, the dominant pore size of the YEC-8A activated carbon (0.87 nm) is slightly larger than the effective molecular dimension of the TFSI− anion (∼0.79 nm).57 The longer fluorinated chain of NFSI−, combined with strong ion pairing and aggregation at high concentration, is therefore expected to further increase the effective size of the charge-carrying species, making ion entry into narrow micropores sterically hindered and energetically unfavourable. Under these conditions, crowding, ion layering, and electrostatic confinement effects dominate, leading to limited pore utilization and suppressed capacitance despite the low scan rate.
This behaviour aligns with prior observations in ionic liquids and highly concentrated electrolytes, where interfacial ion dynamics dominate at high scan rates, enabling efficient charge storage near the electrode surface. At lower scan rates, however, charge storage increasingly depends on ion transport into the internal porous structure of the electrode, which becomes rate-limiting in systems with sluggish ionic mobility. In such cases, especially when pore sizes approach the dimensions of solvated ions, charge storage is further hindered by confinement effects, including ion layering, crowding, and overscreening at electrified interfaces.63 Thus, for 10 m LiNFSI, the charge storage mechanism transitions from surface-limited to diffusion-limited behaviour as scan rate decreases, underscoring the importance of both ion properties and pore accessibility in highly concentrated system.
To further elucidate the impact of pore architecture under conditions where internal pore transport dominates, we compared two carbon electrodes: YEC-8A (average pore size: ∼0.87 nm) and YP-50F (1.1–1.8 nm).64,65 Fig. S16 shows the CV measurements collected at a scan rate of 1 mV s−1 for various LiNFSI concentrations using a 3-electrode configuration. While a direct comparison of specific capacitance between different activated carbons is complicated due to variations in morphology and surface functionality, we can still assess the influence of pore size, particularly when using longer chain ions (NFSI−). At 1 m and 5 m (Fig. S17a and b), YP-50F exhibited slightly higher capacitance than YEC-8A, likely due to its broader pore structure enabling greater ion accessibility (specific capacitance values calculated from CV areas can be found in Fig. S18). Despite these differences, at these concentrations both carbon electrodes performed without significant transport limitations. However, at a 10 m concentration, the performance diverged clearly. As previously discussed, YEC-8A showed severely distorted CVs and diminished capacitance, indicative of hindered ion transport. In contrast, YP-50F maintained a more rectangular CV shape and higher specific capacitance (Fig. S17c and S18), suggesting that its wider pores mitigated the steric and transport challenges posed by bulky NFSI− anions. This comparison emphasizes that ion transport is influenced not only by electrolyte properties but also by the match between pore size and ion dimensions. It is important to note, however, that larger pores are not universally better, e.g., prior studies suggested optimal double-layer capacitance occurs near ∼0.85 nm in aqueous systems.66 Nonetheless, under the extreme conditions of 10 m LiNFSI, i.e., high viscosity, low mobility, and ion clustering, such narrow pores become largely inaccessible. These findings highlight the need for co-optimization of electrolyte chemistry and pore architecture to achieve efficient charge storage under high concentration conditions.
To further validate the influence of electrolyte composition on charge storage behaviour across a range of concentrations, we calculated specific capacitance values from the GCD curves recorded at current densities from 0.1 A g−1 to 10 A g−1 for the positive electrode (0–0.5 V vs. AgQRE) and are presented in Fig. 3d–f. These values corroborate CV trends, following the order LiPFSI > LiTFSI > LiNFSI, with a maximum at 5 m. While GCD-based capacitance values are lower than those from CV (Fig. S19), this discrepancy is expected due to contributions from parasitic currents at voltage extremes, which are more prominent in CV due to its potentiostatic nature.56,67 The capacitance values calculated for the negative electrode (−0.5 to 0 V vs. AgQRE) are also shown in Fig. S19. Although Li+ adsorption occurs predominantly at the negative electrode, the identity of the counterion significantly affects the overall charge screening environment. In highly concentrated systems, such as 10 m LiNFSI, interfacial ion organization is strongly modulated by ion–ion correlations. According to the Kornyshev model and related work on ionic liquids, phenomena such as “overscreening” or “crowding” can occur at the electrode interface, leading to anion-induced changes in double-layer structure and capacitance. These effects are particularly relevant under high current densities, where ion layering and local structuring play an important role in determining device performance.68–70 Fig. 3g–i presents the Nyquist plots measured at 0 V (vs. AgQRE). Notably, all three Li-based electrolyte solutions containing different anions exhibit the steepest slopes at a concentration of 5 m, which can be interpreted as the most ideal EDLC behavior (slope closest to vertical). The observed decrease in slope at low frequencies, despite the measurement being at 0 V where faradaic reactions are not expected, may arise from other interfacial processes such as ion diffusion limitations or subtle surface adsorption phenomena at the electrode/electrolyte interface.71 Such variations in slope may also be influenced by the porous structure of the electrode, interfacial resistance, and ion diffusion processes.
To further investigate the retention of electrolyte species on the electrode surface after electrochemical cycling, post-mortem Raman spectroscopy and X-ray diffraction (XRD) analyses were conducted. Raman spectra indicate that LiTFSI and LiPFSI were almost undetected on the carbon surface after electrochemical experiments (Fig. S20), suggesting their efficient desorption and removal from the electrode. In contrast, LiNFSI could still be detected, even at a 5 m concentration, with an increasing Raman signal at higher concentrations. XRD analysis further supports these findings, showing the presence of residual LiNFSI species on the electrode surface post-cycling, while no crystalline phases related to LiTFSI or LiPFSI were observed (Fig. S21). Although some retention may result from the absence of a washing step, the consistent difference across identically treated samples suggests that LiNFSI exhibits a stronger affinity for the carbon surface as also evidenced by its lower contact angle and improved wettability. While EIS techniques are undoubtedly a powerful tool for elucidating the reaction mechanisms of supercapacitors, in highly concentrated electrolytes, diffusion processes within the electrolyte may predominantly determine the overall behavior, even if redox-active species capable of inducing faradaic reactions are present. Moreover, the behavior of ions adsorbed on the electrode surface often remains relatively stable across much of the frequency spectrum, with noticeable changes typically occurring only in the very low-frequency region.72 Therefore, to further investigate the extent to which the electrolyte penetrates the nanopores of electrode, we conducted additional NMR analysis. This approach allowed us to assess whether the electrolyte effectively infiltrates the micropores of the electrode, or if ion transport is restricted at specific pore sizes, thereby providing a more precise interpretation of the low-frequency behavior observed in the EIS analysis.
To gain insight into the wetting of the carbon surface, magic-angle spinning (MAS) NMR experiments were performed on carbon soaked with LiPFSI. Measurements were carried out by soaking a defined mass of YEC-8a carbon (15 mg) with 20 µL of electrolyte within a MAS rotor. The volume of electrolyte used is approximately 1.3 times larger than the pore volume of the carbon sample. The sealed rotor was then spun at 5 kHz and 1H and 19F MAS NMR spectra were recorded. Fig. 4a shows 1H MAS NMR spectra of carbon soaked with LiPFSI at concentrations between 1 m and 10 m. In each spectrum, two main resonances are observed. Peaks at higher (more positive) chemical shift values correspond to species in excess electrolyte outside the carbon particles (ex-pore). Peaks at lower chemical shifts (more negative) correspond to species in electrolyte adsorbed within the carbon pores (in-pore). The origin of the shift difference between in-pore and ex-pore environments is magnetic field-induced ring currents associated with the delocalized electrons within the carbon surface.73,74
As the concentration increases from 1 m to 10 m, the overall intensity of the spectrum reduces, reflecting the lower water content of the sample. Within each spectrum, the intensity of the ex-pore resonance at ∼4 ppm reduces relative to the intensity of the in-pore resonance at −0.5 ppm. This shows that the water preferentially occupies the in-pore environment. As the total amount of water reduces with increasing concentration, the ex-pore environment is depleted in preference to the in-pore environment. Because the ex-pore resonance overlaps slightly with the in-pore resonance, the intensity of the in-pore resonance also reduces by a small amount with increasing concentration. The magnitude of the ring current shift between the in-pore and ex-pore resonances is approximately 5 ppm, which is consistent with values measured in other studies.38,75,76 For the 10 m concentration, the ring current shift is slightly reduced due to a small reduction in the chemical shift of the ex-pore resonance. This suggests that hydrogen bonding in the ex-pore water is reduced, which is consistent with the relatively low concentration of water remaining in the ex-pore environment which will, on average, see a more LiPFSI-rich environment. Deconvolutions of the 1 m and 10 m spectra (Fig. S22) show evidence of multiple in-pore environments which likely reflects a range of pore environments present in the disordered carbon structure.
Fig. 4b shows 19F MAS NMR spectra of same samples. In each spectrum, two sets of resonances are observed corresponding to the CF2 (−118–−121 ppm) and CF3 (−80–−85 ppm) groups in the PFSI side chains. In these spectra, the overall intensity follows the opposite trend to the 1H spectra, reflecting the increasing PFSI content of the sample as the electrolyte concentration increases. For the 1 m concentration, only in-pore resonances are observed at −85 and −121 ppm, showing that the PFSI− ions preferentially occupy the in-pore environment and a negligible proportion remain in the ex-pore environment. This indicates that the ex-pore H2O observed in the 1H MAS NMR spectrum is essentially free water that is depleted of LiPFSI. As the electrolyte concentration increases, ex-pore resonances at −80 and −118 ppm emerge and increase in intensity. This suggests that the emergence of the ex-pore resonance arises from pore saturation, where, above a certain loading and concentration, there is no space remaining in the in-pore environment and the LiPFSI is forced to occupy the ex-pore environment. This trend is also mirrored by the 7Li MAS NMR spectra of the same samples (Fig. S22), which shows that the partitioning of the Li cations between the in-pore and ex-pore environments broadly follows that of the PFSI anions. In the 5 m spectrum, the intensities of the in-pore resonances appear higher than in the 10 m spectrum; however, deconvolutions reveal the in-pore ion populations for both concentrations are very similar (Table S2), further evidencing saturation of the in-pore environment. Deconvolutions also enable the local ion
:
water molar ratios within the in-pore and ex-pore environments to be determined (Table S2).38 These can be normalized by the ionic molar ratios of the neat solutions, where a value below 1 shows that the local concentration is reduced relative to the neat solution, and a value above 1 shows that the local concentration is increased relative to the neat solution. As shown in Fig. 4c, the normalized in-pore ionic ratio reduces with increasing concentration, while the normalized ex-pore concentration increases with increasing concentration. The crossing point appears to be at a concentration of approximately 7 m, which corresponds to where the H2O and LiPFSI would be expected to be equally distributed across the system.
Overall, the NMR results show that LiPFSI preferentially occupies the in-pore environment, and only occupies the ex-pore environment when the pores become saturated with high-concentration loadings. This is broadly in line with the HOPG contact angle measurements which demonstrate the high affinity of LiPFSI for the carbon surface, although the pronounced differences in the in-pore electrolyte composition viewed by NMR highlight the subtleties of the relative affinities of LiPFSI and H2O for the carbon surface in a way which is not evident from the bulk contact angle measurement. Indeed, the carbon-LiPFSI system exhibits a degree of hydrophobicity which appears to depend on the bulk electrolyte concentration. At 1 m concentration when all of the LiPFSI is inside the pores, it can accommodate a smaller fraction of the H2O than when in the bulk, resulting in some H2O residing outside the pores. At 10 m concentration, almost all of the H2O is accommodated inside the pores, although the total amount of H2O present here is much lower than in the 1 m sample. In the 5 m sample, the in-pore environment contains the highest combined amount of LiPFSI and H2O of the samples studied. It should be noted that in a functioning capacitor, the electrolyte is present in excess relative to the carbon, so there will be ex-pore H2O and LiPFSI in all cases. However, the in-pore electrolyte composition in the absence of an applied potential is expected to be broadly comparable to the observations made here. This may indicate that the maximum capacitance observed at ∼5 m concentration corresponds to where the in-pore electrolyte has an optimum balance between the number of ions present to store charge, and the amount of H2O present to facilitate movement of the ions through the porous network during charging.
These results provide practical guidelines for designing high-voltage aqueous supercapacitors: selecting mid-concentration electrolytes and matching them with suitable carbon pore sizes optimizes ion transport and charge storage, highlighting the importance of co-designing electrolyte chemistry and electrode architecture.
These findings highlight that carbon pore accessibility, alongside ion size and solvation behavior, is a critical factor in determining performance. Future design of electrolytes for supercapacitor applications should therefore consider not only bulk physicochemical properties, but also how ions interact with and distribute within the nanoporous electrode structure. This integrated understanding is essential for advancing high-concentration aqueous electrolytes tailored specifically for next generation supercapacitor technologies.
Supplementary information (SI): data on bulk solution properties, contact angles of the electrolytes, further Raman, XRD and NMR data, along with more detailed electrochemical data. See DOI: https://doi.org/10.1039/d5ta09457g.
| This journal is © The Royal Society of Chemistry 2026 |