High-energy-density energy storage technology requires the use of a new generation of anode materials, e.g. alkali metal anodes (Li, Na or K) with high theoretical capacity, and in the case of Na and K, also having a higher abundance in the earth’s crust. Such reactive anodes necessitate the discovery of compatible electrolyte materials to support the safe and long-term cycling of these new energy storage devices. Polymer electrolytes (PEs) are now accepted as an enabler of the ultimate solid-state high-performance batteries.1-4 Given the appropriate chemistry and materials properties, these materials can address the essential safety issue, i.e. prohibiting the formation of dendrites on highly reactive metal anodes through their optimal mechanical properties, in addition to formation of a potentially beneficial solid electrolyte interphase (SEI) layer, depending on their chemistry.5 Furthermore, PEs possess the additional advantages of high thermal stability, high electrochemical stability and ideal flexibility that help to improve battery safety and cycle life.
The earliest PEs for battery applications were proposed based on archetypal polyethylene oxide (PEO) or copolymers of PEO.6, 7 The polar ether functional group in the polymer backbone coordinates Li+ and drives its diffusion through the strongly coupled backbone motion. These PEO based systems have been normally restricted to using relatively low salt concentrations (e.g. EO: Li+ = 20) in order to obtain maximum ionic conductivity, 8 although the value is still below the minimum conductivity required for practical applications at moderate temperatures. They also suffer low Li+ transference numbers (TLi+ ~ 0.2) as the ‘free’ anions, such as bis(trifluoromethanesulfonyl)imide (TFSI–), dominate the conductivity whereas the Li+ motion is restricted by the strong EO-Li coordination.9-11 Although PEO based PEs have been commercialized by Bolloré group, their use is still limited, in particular to positive electrode material operating below 4V vs Li+/Li°. A variety of polymer design strategies have been proposed to improve PE performance, e.g. designing new anion chemistries with restricted motions,12 using polycarbonate groups to enhance dissociation of lithium salt, 11, 13, 14 using block copolymers and crosslinking of polymer chains to improve mechanical properties and ionic conductivities, etc..15, 16 Besides those, the search for new polymer materials with higher metal ion transference number and conductivity is ongoing.
One approach to maximise the transference number of the metal ions is to use polymers that have any other ions chemically tethered to backbones (also known as polyelectrolytes) to immobilize their motion. In the case of polyanionic materials, in which anions are chemically bonded to polymer backbones, cation conduction dominates and the cation transference number tends to unity.17 Thus in the case of a lithium salt polyelectrolyte, the ionic conductivity is completely derived from Li+ conduction which is highly desirable in lithium batteries in order to overcome problems related with anion concentration gradients. Despite progress such single ion conductors with any plasticizers/solvents still typically suffer very low ionic conductivities that hinder their applications. 18, 19
Recently, polycationic PEs have emerged as potential solid-state solvents for lithium salts, in particular the monomer is a polymerizable IL cation, showing good performance in Lithium batteries.20, 21 This material is termed a polymeric or polymerized ionic liquid (PolyIL),22 which inherits both the excellent electrochemical performance of ILs and good physical, thermal and mechanical properties of polymers.23, 24 Interestingly, these materials have been shown to support the dissolution of extremely high lithium salt concentrations, such that the ratio of Li+: polycation is greater than 1, without sacrificing ionic conductivity,21, 25, 26 due to a unique Li-anion-polycation co-coordination structure (Li+ ions are distributed around the polycation backbone by anion-bridging26). This abnormal concentration - conductivity positive correlation differentiates PolyILs significantly from traditional PEs such as PEO-based electrolytes,27 but is in line with the ‘polymer-in-salt’ electrolytes (PISE)28-30 proposed by C. Austen Angell decades ago. In such cases, ion mobility will begin to increase again while the glass transition temperature (Tg) decreases, once the salt is the dominant component in the PISE system.28 However, the PolyIL is unique in that this salt concentration-conductivity positive correlation is also applicable to the low salt concentration range.21, 26
The fast Li+ conduction in “polyIL-in-salt” should be closely related to its transport mechanism. A recent report using computational methods, suggested a vehicular Li+ diffusion mechanism in PolyIL electrolytes which, however, only focused on a low to medium salt concentration range (i.e. Li+ : cation repeat unit £ 0.4).31 It is likely that this mechanism is concentration dependent, as has been noted for liquid electrolytes, since the salt concentration will alter the Li+ coordination environment.32-36 Given the recent promising demonstration of polyIL-in-salt systems in Li metal batteries, understanding the transport mechanisms in these new PISE materials could open up the possibility of designing PEs with high conductivity and high target ion transference number. Furthermore, while several systems have been proposed for Li devices, materials focusing on other battery chemistries such as Na and K are sparse and would enable solid-state batteries based on these more abundant materials.
Thus, in this work we used a computational approach to investigate the transport of Na+ and K+ in PolyIL based PISEs and discovered a rapid diffusion of alkali metal ions through a structural diffusion mechanism. Under the premise of weak interaction between metal ions and anions, the salt-concentrated environment favors fast metal ion diffusion. We also investigated the solvation and diffusion of the multivalent Mg2+ in these PolyILs and discussed the challenges in achieving high Mg2+ conduction in such systems. Finally, as a proof of concept, we validated the computational findings by experimental measurement of conductivity as well as in application of a Na based PolyIL PISE in a sodium symmetric cell, which showed excellent stable cycling behavior at a capacity of 0.5mAh cm-2 per cycle for over 100 cycles.
Computational design optimal polyIL-in-salt systems
The PolyIL studied in this work is poly(diallyldimethylammonium)bis(fluorosulfonyl)imide (PDADMA FSI) (Supplementary Fig1e). Our previous experimental work has studied the use of this polymer with increasing amounts of LiFSI salt at polycation to Li+ mole ratios from 2:1 to 1:6. Careful experimental measurements were used to find the optimum conductivity.26 However, such an approach to electrolyte design is both time consuming and costly and here we sought a more efficient computational approach. Based on MD simulations in that earlier work we learned that the optimal salt concentration (1:1.5) determined in PDADMA FSI/Li FSI allows the most lithium ions to be homogeneously distributed around polycationic backbones by co-coordinating with anions. At higher Li+ ratio, excessive Li+ form molten salt-like aggregates. Although in silico, these previous simulations suggested that Li+ diffusion could be further enhanced with these aggregates, this could not be realised in our experiments as those Li+-rich domains eventually transform into the less conductive crystalline salt-like new phase. Nevertheless, based on these understandings, we hypothesize here that a close-to-optimal salt concentration could at least be identified from a PolyIL PISE composition that gives the ‘saturated’ metal ion distribution around polymer chains, and can be determined via quantifying three different anion coordination states: (1) the polycation-anion single-coordination state (P+-A–, Fig 1a) that can be seen in the case of neat polyIL or with low salt concentrations; (2) the metal ion-anion single-coordination state (Me+-A–, Fig 1c) as seen in the molten salt-like regime; and (3) the polycation-anion-metal ion co-coordination state (P+-A–-Me+, Fig 1b). We use this to investigate alternative salt chemistries including NaFSI, KFSI and Mg[FSI]2.
The PolyIL/NaFSI systems at two (P+: Na+) ratios of 1:2 and 1:4 (termed Na12 and Na14) were firstly examined. The 1:2 cation ratio was chosen for the NaFSI system as it is close to the optimal 1:1.5 ratio of the lithium system, and Na+ requires to be coordinated with slightly higher number of FSI– than Li+. The 1:4 cation ratio as a higher salt concentration was also examined. The snapshots in Fig 1d-1e show the selected structure of the equilibrated Na12 and Na14 systems, including FSI– only in two single-coordination states (purple for P+-FSI– and green for Na+-FSI–), the cationic polymer matrix without H atoms (aqua blue) and the Na+ ions (yellow). The lowest percentage of P+-FSI– or Na+-FSI– coordination occurs in the Na12 system, with the highest state being the co-coordination state consisting of 97% FSI–, which is even higher than that (92%) in the optimal PolyIL/LiFSI (1: 1.5) system. This high percentage of co-coordination suggests that the 1:2 ratio can be selected as the initial concentration for the Na+ system. For the Na14 case, with the added salt doubled, the percentage of FSI– in the Na+-FSI– coordination state increases significantly to 19.7%, indicating the growth of the Na+ - FSI– rich domains at higher Na+ concentration, and this change is consistent with the LiFSI case26.
The 1:2 cation ratio was examined next for both K+ and divalent Mg2+ systems to investigate anion states and the possible optimal salt concentrations. As shown in Fig 1f and the snapshots in Supplementary Fig 1, there is an obvious increase in the percentage of FSI– in FSI–-Me+ coordination state with Na+ < K+ < Mg 2+ (2.2%, 8% and 23.1% respectively), while the co-coordination state decreases to 91% and 73.9% for K12 and Mg12, respectively. This suggests us to test a lower salt concentration, with the expectation of increasing the co-coordination state of FSI– as confirmed in Fig 1f. We found that, at 1:1 ratio, the co-coordination state increases to 94% and 85.4% in the K11 and Mg11 systems, respectively. In the K11, the percentage of FSI– in both single-coordination states is relatively low, while in the Mg11, the percentage of both is 7.3%, which may be due to the difficulty in completely dissociating the Mg salt.
Interestingly, when comparing cation ratios that can obtain nearing the maximum co-coordination anion state, i.e. 1:1.5 for Li+, 1:2 for Na+ and 1:1 for K+, they are ordered differently to alkali metals in the periodic table (Li < Na < K). This variation in cation ratio could be related to the difference in cation-anion coordination abilities. For example, for a given 1:2 cation ratio, the coordination number (CN) (Supplementary table 1) of P+- FSI– in the Na12 (8.7) is greater than that in the K12 (8.1), but an opposite situation is observed for the CN of Me+- FSI– coordination which is greater in the K12 (6.4) than in the Na12 (5.6). Therefore, the number of metal ions in the co-coordination state bridged by FSI– should also be less in the K12 and more in the Na12 system. The CN also helps to understand the extensive Mg-FSI enriched domains in the Mg12 system. Since the amount of FSI– is doubled whereas the CN of both Mg2+-FSI– (5.7) and P+-FSI– (9.2) are not, this results in excessive FSI– anions that mainly coordinate with Mg2+. A more detailed structural analysis of the cation-anion coordination by radial distribution function (RDF, supplementary Fig 2) can be found in the supplementary information.
The impact of chemical environmental on ion diffusion
In this section, the diffusion of ions is investigated to evaluate the ionic conductivity in sodium, potassium and magnesium electrolytes. The ion diffusion is compared among three electrolytes of K12, Na12 and Mg11 having the same concentration of anions via the mean square displacement (MSD) of metal ions and anions (Supplementary Fig 3). In general, ion diffusion in the alkali metal ion systems is 1-2 orders of magnitude higher than that in the magnesium system. At a temperature of 353 K, the diffusion of metal ions and anions in the K12 system is faster than that in the Na12 system. Similar results were also reported in an ionic liquid electrolyte N-methyl-N-propylpyrrolidinium bis(fluorosulfonyl)imide (C3mpyrFSI). When mixed with the same amount of KFSI or NaFSI salt, the IL system with the KFSI salt has higher ionic conductivity.37 This may be due to the lower Lewis acidity of potassium, making its interaction with the anion weaker than sodium, which is confirmed by the density functional theory (DFT) calculation of the binding energy (Eb) between metal ions and anions.
The relative diffusivities of metal ions and anions in monovalent alkali metal salts are different from those in divalent magnesium salt. The diffusivity of alkali metal metal ions is higher than that of anions. Taking into account that polycations are immobile, this indicates that a higher alkali metal ion transference number can be obtained, which is in sharp contrast with traditional PEO-based electrolyte systems. This was confirmed in the case of Na via experimental validation in next section. On the other hand, the opposite behaviour is observed for the Mg2+ system, which shows a lower Mg2+ diffusivity relative to FSI– in addition to being the least diffusing cations amongst all simulated systems. Next, we conducted a detailed structural analysis of the targeted ions having distinct fast and slow diffusion, and studied the environmental factors that affect the fast ion diffusion. The Na12 and Mg11 systems were selected for this analysis.
Firstly, the mode of interaction between the arbitrary selected fast and slow diffusive ions (defined in supplementary Fig 4) was studied by calculating the RDF between the two selected groups (i.e. fast vs. fast or fast vs. slow group) as presented in Figure 2. Considering firstly the correlation between FSI– and Na+ diffusion, Fig 2c and 2d show that fast FSI– tends to coordinate with fast Na+, and similarly, slow FSI– tends to coordinate with slow Na+, which is inferred from their more prominent RDF peaks and higher CN in comparison to the coordination between the fast and slow ions. This also indicates that the movement of Na+ and FSI– has a strong correlation. Secondly, Fig 2e and 2f show the RDF of polycation - FSI– and polycation - Na+, respectively, revealing that the fast FSI– (or Na+) has less coordination with polycations, which is in contrast to the slow FSI– (or Na+). This suggests that the diffusion of ions is slower near polycations or when coordinating with more polycations. Finally, Fig 2g and 2h investigate the number of Na+ and FSI– ions surrounding the selected fast and slow Na+. Interestingly, there are always more Na+ and FSI– ions surrounding the fast Na+, i.e. the environment rich in Na+ and FSI– facilitates the faster metal cation dynamics, as demonstrated by a snapshot in Fig 2b. The fast Na+ (green) has the molten Na-FSI salt-like environment which is commonly seen in superconcentrated ionic liquid electrolytes.32, 36, 38 In contrast, the environment of slow Na+ (blue) contains fewer FSI– and Na+ within the same distance range. This result suggests that the evolution of molten salt regions at high salt concentrations can promote metal ion diffusion in the polyIL, thus providing a theoretical guideline for future polymer electrolyte design, i.e. by enhancing the molten salt like domains in a polyIL matrix. We experimentally validate this concept later in the paper.
In the case of the Mg11 system we observe some contrasting behaviour. The fast FSI– preferentially coordinates with the polycationic backbone (Supplementary Fig 5a - 5b). Furthermore, the environment of the fast Mg2+ is also dominated by more polycations, with a lower CN for both FSI– and Mg2+ (Supplementary Fig 5d - 5f). In other words, the faster Mg2+ and FSI– are not found in domains rich in Mg2+ and FSI- ions, which is the opposite situation to that shown in the alkali metal salt systems. We can understand this by considering the relative binding interaction between the metal ions and FSI– (Fig 2b); for Mg2+ this is so strong that it dramatically reduces the diffusion of both ions. Therefore, in these high salt content PolyIL electrolytes, the ionic conductivity enhancement requires the metal ion-anion interaction to be weak enough to allow fast exchange in the coordination environment, which we discuss in more detail below. For the development of polyIL electrolytes that conduct multivalent ions (e.g. Mg2+, Zn2+ etc.), designing metal ion-anion interaction of similar magnitude to the alkali metal ion -FSI binding energy is the key.
Metal ion diffusion mechanisms
In this section, we show that the dominant metal ion diffusion mechanism in the polyIL-in-salt system, is a structural diffusion mechanism as opposed to the vehicular diffusion mechanism suggested in the low salt concentration polyIL31, or the diffusion through polymer backbone motion in the traditional PEs (e.g. PEO). If we consider the nearest coordination structure of a given metal cation to be its ‘cage’ we can then quantitatively determine the cage reorganisation dynamics and their relationship to ion diffusion processes. Three types of metal ion diffusion mechanisms can be distinguished; vehicular, structural and hopping. We define the period from the initial formation of an old cage A to the formation of a new cage B as the ‘cage-restructuration period’ (Fig 3a). This period will be longer for the vehicular mechanism since the metal ions will diffuse with the cage. In contrast, in the case of the structural mechanism, this time period will be much shorter since the diffusion of ions occurs through the cage restructuring. The cage restructuration in the hopping mechanism occurs instantaneously and causes ion migration via an ion hop from one cage environment to another.
In Fig 3b The ion cage restructuring in Na12 and K12 systems is analysed through the cage-restructuration correlation function C(t) (with a specific definition described in the supplementary file). The decay in this function is caused by the occurance of cage restructuring, i.e. new anions participate in metal ion’s coordination shell. C(t) decays to 0.9 almost immediately for both K12 and Na12 systems. It takes nearly 3.5 ns to decay to close to 0.01 in the K12, which is significantly faster than the 12.7 ns required in Na12 system for most metal ions to complete one cage restructuration process.
Fig 3c and 3e analyse the frequency of different cage-restructuration periods in the 30 ns trajectory, and accumulates the deviation distances of metal ions in all the same cage-restructuration periods. The average deviation distance of each restructuration period is given in Fig 3d. The accumulative distance profile and the frequency profile present a similar pattern, indicating that the distance the metal ion travels is positively correlated to the frequency of the cage-restructuration period, but has no relationship with the cage-restructuration time, as suggested by Fig 3d. This further supports the structural diffusion mechanism for metal ions, and demonstrates that the key to improve the metal ion conduction is to accelerate the reorganization of their coordination shells.
Finally, the frequency of an instant hopping event was identified in order to evaluate the hopping diffusion mechanism. This is when the cage restructuration is completed between two consecutive frames in the MD trajectory. Only 7 events were detected in more than 2 million inspections, indicating that his mechanism would only rarely contribute to diffusion, and thereby ionic conductivity, in these electrolytes.
Experimental validation of in-silico predictions
Based on simulation results, we conducted experimental validations to verify the composition dependence and high ion conductivity in these PolyIL PISE electrolytes. The Na11, Na12, Na14, K11 and K12 electrolytes were prepared and their thermal properties, ionic conductivities and electrochemical behaviour were studied. Fig 4a presents the Tg of those systems with that of the Li11 and Li12 systems from our previous study.26 Tg is normally used as an indicator of ion dynamics for those traditional electrolytes whose ion dynamics are highly correlated with the relaxation process of the polymer backbones. For all polyIL systems studied here in Fig 4a, Tg decreases with the increase of metal ion concentration, for example, Tg follows Na11 >Na12 > Na14. The decrease in Tg is due to the increase in metal ion-anion-polycation co-coordination so that the polycation-anion interaction and the ionic crosslinking of polymer chains is reduced, thereby increasing the local dynamics that controls Tg. Excitingly, the lowest Tg of Na14 could support the modelling prediction regarding the positive role of the molten salt-like regions on ion dynamics in such systems.
The conductivity was measured for Na 11, Na 12, K11, and K12 systems in Fig 4b. Na12 has higher conductivities than Na11 over the whole temperature range, and sNa12 reaches a value around 1.0 mS cm-1 at 80 oC, indicative of fast ion dynamics in this electrolyte. The conductivities of K11 and K12 systems measured during the 1st heating cycle, however, do not reflect their difference in Tg. The conductivity of K11 is higher than that of K12 except at a temperature of 100 oC, although the latter has a lower Tg. Both are also less conductive than their sodium counterparts, which is also contrary to MD predictions. The DSC measurements on the two K-based systems (Supplementary, Fig 6a) suggest that this lower conductivity is likely due to the presence of a second crystalline phase, as we previously noted for Li systems. Encouragingly, a dramatic increase in the conductivity of K12 is seen during the 2nd heating cycle, when the second crystalline phase is completely eliminated, as indicated by the disappearance of the melting peak in DSC (Supplementary Fig 6b). The conductivity of K12 is then significantly higher than K11 above 60 oC and also surpasses Na12 above 80 oC, and this result is consistent with our MD predictions at 353 K. It also confirms the conclusion that maintaining a single, amorphous polymer phase is essential to ensure high ionic conductivity in such high salt concentration PEs. It should also be noted that despite all PolyILs having the same FSI anion, their measured Tg is not completely correlated with the ionic conductivity across PolyILs with different alkali metal salts. For example, although the Li12 has the lowest Tg in Fig 4a, its ionic conductivity is around 0.6 ´ 10-4 S cm-1 at 80 oC, which is lower than that of Na12 (1.02 ´ 10-3 S cm-1) or K12 (1.2 ´ 10-3 S cm-1) (compared at 2nd heating) .
Fig 4c shows the chronoamperometry profile of a symmetric cell assembled with the Na12 electrolyte and the Nyquist plots before and after DC polarization (inset). Although there is a fast decay of current in the first hour, the current is still maintained at a high level, indicative of relatively high Na+ current (initial current 0.054 mA vs steady state current 0.044 mA). Further analysis via the Bruce-Vincent method confirms an Na+ transference number of 0.57 at 80 oC, which is significantly higher than that of PEO-NaFSI system (tNa+=0.16) as previously reported by Zhou et al. 39
The Na plating/stripping behaviours were also evaluated by symmetric Na cell cycling tests. As shown in Fig 4d, the Na12 electrolyte can sustain long term, stable Na plating and stripping under a high current density of 0.5 mA cm-2, with areal capacity of 0.5 mAh cm-2. It should be noted that, for solvent-free electrolytes, the current density applied in this case is much higher than that of state-of-the-art electrolyte systems. For example, a current density of 0.1 mA cm-2 was used for PEO/NaFSI electrolytes by Hu et al.40 More surprisingly, the polarization voltage for the Na12 electrolyte is relatively low, around 100 mV (inset in Figure 5b), even at such a high current density. We ascribe this superior and stable plating/stripping performance to a high ionic conductivity and high transference number in the PISE. The preliminary test on the Na14 system (Supplementary Fig 7) also show very promising results with the higher steady state Na+ current at 0.064 mA (initial current 0.084 mA) and higher Na transference number of 0.63 obtained at 80 oC, and overpotential when cycling at high current densities. However, in this case the molten salt regime is metastable and transforms through crystallization after a few days, and thus future work will focus on stabilizing such systems.