Importance of liquid fragility for energy applications of ionic liquids.

Ionic liquids (ILs) are salts that are liquid close to room temperature. Their possible applications are numerous, e.g., as solvents for green chemistry, in various electrochemical devices, and even for such "exotic" purposes as spinning-liquid mirrors for lunar telescopes. Here we concentrate on their use for new advancements in energy-storage and -conversion devices: Batteries, supercapacitors or fuel cells using ILs as electrolytes could be important building blocks for the sustainable energy supply of tomorrow. Interestingly, ILs show glassy freezing and the universal, but until now only poorly understood dynamic properties of glassy matter, dominate many of their physical properties. We show that the conductivity of ILs, an essential figure of merit for any electrochemical application, depends in a systematic way not only on their glass temperature but also on the so-called fragility, characterizing the non-canonical super-Arrhenius temperature dependence of their ionic mobility.

Various new application-relevant material classes have been discovered during recent years. Among them, ILs maybe sparked the largest interest due to their numerous applications12345, and the number of publications treating these materials by far exceeds those on such prominent topics as high-T superconductivity, colossal-magnetoresistance or multiferroicity. For example, as ILs are exclusively composed of ions, they are good candidates for any applications where high ionic conductivity is needed. They also fulfil various other requirements as a broad electrochemical stability window, low vapour pressure or non-flammability34. The vast number of possible combinations of anions and cations, nowadays known to form ILs, opens up many possibilities for finding compounds optimized for application. Various systematic investigations of the dependence of the conductivity on composition have been performed, e.g., by varying the anion for ILs with the same cation or vice versa678910. However, when considering the fact that there is an estimated number of one million binary ILs1, obviously these investigations barely have scratched the surface.

Under cooling, ILs usually solidify via a glass transition, a continuous increase of viscosity instead of the abrupt crystallization via a first-order phase transition found in most electrolytes. Thus, the non-canonical dynamic properties of glass-forming matter1112, dominate many of their physical properties13141516. This is also true for their dc conductivity σdc (refs. 7, 8, 9, 10,13,14,16, 17, 18, 19), which exhibits the typical temperature dependence known to govern the dynamics of molecules, ions or any other constituents forming glassy matter. In most ionic conductors, the Arrhenius law, σdc ∝ exp[-E/(kBT)], characteristic for thermally activated ion hopping over an energy barrier E, provides a good description of this temperature dependence. However, in ILs it can be parameterized by the empirical Vogel-Fulcher-Tammann (VFT) formula known from glass physics, usually written in the modified form11:

Here D is the so-called strength parameter, which is used in the classification scheme for glass formers, introduced by Angell11, to distinguish between so-called strong and fragile glass formers. While the latter exhibit marked deviations from Arrhenius behaviour (small D; solid lines in Fig. 1), strong glass formers (dashed lines) more closely follow the Arrhenius law. Another, nowadays more common way to parameterize these deviations is the fragility index m (ref. 12), defined as the slope at T in the Angell plot11, log y vs. T/T, where T is the glass temperature and y is any quantity coupling to the glassy dynamics (viscosity, relaxation time or, in our case, conductivity). D and m are related by m = 16 + 590/D (ref. 12).

Figure 1

Influence of glass temperature and fragility on the room-temperature conductivity of ionic liquids.

The figure shows the temperature dependence of the conductivity (in Arrhenius representation) of four hypothetical ionic glass formers A–D with different glass temperatures and fragilities. Samples A and B, having identical, high fragilities (m = 150 or D = 4.4), demonstrate the effect of different glass temperatures. Glasses A and C (or B and D) having identical glass temperatures but different fragilities (glasses C and D have m = 30 and D = 42) demonstrate the significant influence of the latter quantity on the conductivity. The dotted line indicates room temperature.

It is obvious and well known that the glass temperature T has an immediate effect on the conductivity of ILs20. It marks the solidification of a glass former when its viscosity reaches about 1012 Pa s or when the time scale of the molecular motions becomes larger than 100–1000 s. Low glass temperatures imply low viscosities at room temperature and in ILs viscosity and conductivity are coupled, even though the degree of coupling can vary1316192122. Thus the conductivity becomes enhanced for small T. This is illustrated by curves A and B in Fig. 1, showing temperature-dependent conductivities of four hypothetical ILs. (For many ionic conductors, σdc(T) ≈ 10−15 Ω−1 cm−1 (ref. 23) and, thus, in Fig. 1 T is determined by the point where the curves reach the abscissa). However, based on the phenomenology of glass physics, the fragility should also play an important role (cf. curves A and C in Fig. 1). Indeed, the dependence of the room-temperature viscosity on fragility was already pointed out1316. In the present work we demonstrate that the fragility also has strong impact on the conductivity of ILs.

Results and Discussion

We have measured the dielectric response of 13 ILs. Covering frequencies from about 10−1 to 109 Hz and a wide temperature range, extending deep into the liquid regime and approaching T at low temperatures, enables a thorough analysis of the dc conductivity and the fragility. As a typical example, Fig. 2 shows spectra of the real and imaginary part of the dielectric permittivity (ε′ and ε″, respectively), the conductivity (σ′) and the imaginary part of the dielectric modulus24 (M″) of Omim PF6 (for a definition of the sample abbreviations, see Table 1) for selected temperatures. It should be noted, that partly the information contained in these plots is redundant (e.g., σ′ ∝ ε″ ν). However, the various dielectrically active processes of ILs are differently emphasized in these plots, making their separate discussion helpful.

Figure 2

Dielectric spectra of Omim PF6.

Spectra are included for a variety of temperatures. The shown quantities are: Dielectric constant (a), dielectric loss (b), conductivity (c) and the imaginary part of the dielectric modulus (d). The lines in (a,b) are fits assuming a distributed RC circuit to model the blocking electrodes25, dc conductivity and three relaxational processes described by the Cole-Davidson or Cole-Cole functions. ε′(ν) and ε″(ν) were simultaneously fitted. The lines in (c,d) were calculated from the fits to ε′ and ε″.

Table 1

Parameters of the ionic liquids included in Fig. 4.

Ionic liquid	Tg (K)	ρdc (Ω cm)	m	ref.
1-Methyl-3-octylimidazolium hexafluorophosphate (Omim PF6)	194	3.6 × 103	78
1-Hexyl-3-methylimidazolium hexafluorophosphate (Hmim PF6)	192	2.1 × 103	84
1-Butyl-3-methylimidazolium hexafluorophosphate (Bmim PF6)	189	7.7 × 102	92
1-Butyl-3-methylimidazolium bis-(trifluoromethylsulfonyl)imide (Bmim TFSI)	181	2.1 × 102	88
1-Butyl-3-methylimidazolium tetrafluoroborat (Bmim BF4)	182	2.3 × 102	93
1-Butyl-3-methylimidazolium chloride (Bmim Cl)	228	4.2 × 104	97
1-Butyl-3-methylimidazolium tetrachloroferrate (Bmim FeCl4)	182	1.5 × 102	144
1-Butyl-3-methylimidazolium bromotrichloroferrate (Bmim FeCl3Br)	180	1.1 × 102	146
1-Benzyl-3-methyl-imidazolium chlorid (Benzmim Cl)	253	1.4 × 106	78
1-Ethyl-3-methyl-imidazolium tricyanomethanide (Emim TCM)	183	4.6 × 101	158
1-Butylpyridinium tetrafluoroborat (Bpy BF4)	195	4.6 × 102	117
(1-Butylpyridinium)0.6(1-Butyl-3-methylimidazolium)0.4 tetrafluoroborat ([Bpy+Bmim]BF4)	191	3.6 × 102	111
1,3-Dimethylimidazolium(Li 1.0m) bis-(trifluoromethylsulfonyl)imide ([Li+Dimim]TFSI)	202	6.9 × 102	145
1-Propyl-3-methylimidazolium tetrafluoroborate	175	1.9 × 102	78	8
1-Butyl-3-methylimidazolium tetrafluoroborate	178	2.8 × 102	93	8
1-Pentyl-3-methylimidazolium tetrafluoroborate	183	6.1 × 102	78	8
1-Hexyl-3-methylimidazolium tetrafluoroborate	188	8.5 × 102	66	8
1-Hepyl-3-methylimidazolium tetrafluoroborate	186	1.5 × 103	68	8
1-Octyl-3-methylimidazolium tetrafluoroborate	190	1.7 × 103	62	8
1-Nonyl-3-methylimidazolium tetrafluoroborate	191	2.4 × 103	55	8
1-Hexyl-3-methylimidazolium bis-(trifluoromethylsulfonyl)imide	187	5.9 × 102	57	8
1-Hexyl-3-methylimidazolium chloride	220	3.7 × 104	64	9
1-Hexyl-3-methylimidazolium bromide	216	1.1 × 104	71	9
1-Hexyl-3-methylimidazolium iodide	208	6.3 × 103	79	9
1-Hexyl-3-methylimidazolium hexafluorophosphate	194	2.4 × 103	94	9
1-Hexyl-3-methylimidazolium tetrafluoroborate	187	6.9 × 102	89	9
1-Butyl-3-methylimidazolium bromide	221	1.3 × 104	69	10
1-Butyl-3-methylimidazolium iodide	215	3.6 × 103	55	10
1-Butyl-3-methylimidazolium thiocyanate	195	2.4 × 102	56	10
1-Butyl-3-methylimidazolium tetrafluoroborate	189	3.1 × 102	56	10
1,3-Dimethylimidazolium dimethylphosphate	201	4.5 × 102	110	17
1,5-Bis(3-benzyl-2-methylimidazolium)pentane di-bis(trifluoromethanesulfonyl)imide	250	7.1 × 104	173	18
1,10-Bis(2,3-methylimidazolium)decane di-bis(trifluoromethanesulfonyl)imide	225	4.3 × 104	141	18
1,10-Bis(3-methylimidazolium)decane di-bis(trifluoromethanesulfonyl)imide	213	1.7 × 104	141	18
1,5-Bis(3-methyl-2-phenylimidazolium)pentane di-bis(trifluoromethanesulfonyl)imide	251	4.0 × 104	168	18

Glass temperature (T), dc resistivity at room temperature (ρdc) and fragility (m) are listed. For the materials taken from literature, the reference number is provided. The other liquids were measured in the present work.

In ε′(ν) (Fig. 2a) the huge increase and the approach of a plateau at low frequencies is due to electrode polarization or blocking electrodes arising from the trivial fact that the ions cannot penetrate into the metallic capacitor plates25. The steep low-frequency increase in the loss, ε″(ν), up to about 105 (Fig. 2b), is caused by the dc conductivity σdc and corresponds to the plateau observed in σ′(ν) (Fig. 2c). The levelling off of ε″(ν) at the lowest frequencies and the corresponding decrease in σ′(ν) also arises from blocking electrodes. Curiously, in ionic conductors ac measurements are necessary to determine the dc conductivity.

A closer inspection of the rather weakly frequency-dependent regions, observed in ε′(ν) and ε″(ν) at the higher frequencies, reveals signatures of relaxational processes, i.e. peaks in the loss, accompanied by steps in the real part, both shifting towards higher frequencies with increasing temperatures. This resembles the reorientational relaxations found, e.g., in molecular dipolar glass formers26. Such behaviour has also been reported for ILs791427 where it can be ascribed to the reorientational motions of dipolar ions and the corresponding secondary relaxations14. In σ′(ν) (Fig. 2c) these modes cause the observed increase at high frequencies, following the dc plateaus. A more detailed assessment of these dipolar modes is out of the scope of the present work.

M″ spectra as shown in Fig. 2d are often used to define a so-called conductivity-relaxation time τ, thought to characterise the ionic dynamics24. The main peak in M″(ν) is ascribed to the translational ionic motions and τ is determined from the peak frequency v via τ = 1/(2πv). The modulus representation suppresses the effects of blocking electrodes showing up at low frequencies in the other quantities. It should be noted that this main peak in M″(ν) is the conductivity-relaxation peak. It is not related to the reorientational relaxation processes (treated in the previous paragraph and revealed by the peaks in ε″(ν)) but to the translational ionic motions. Any sizable charge transport in a material generates a prominent peak in M″(ν) [24]. In ionic liquids with dipolar ions, there is both, translational charge dynamics (leading to a peak in M″(ν) and a plateau in σ′ (ν)) and reorientational dynamics (leading to one or several peaks in ε″(ν)). However, it is known that relaxation peaks observed in ε″(ν) (which arise from reorientational motions) also give rise to peaks in M″(ν), but with somewhat shifted peak frequency28. This explains the additional peaks and shoulders observed in Fig. 2d at higher frequencies, which have the same origin as those found in Fig. 2b. As mentioned above, these dipolar reorientational modes are out of the scope of the present work and in the following we concentrate on τ, determined from the prominent M″ peak and on the dc conductivity σdc. One should be aware that the significance of the modulus representation is rather controversial but, nevertheless often used in ILs714171921.

We have fitted the spectra of Fig. 2a,b assuming a distributed RC circuit for the blocking electrodes25, dc conductivity and three reorientational modes (the lines in Fig. 2c,d were calculated from these fits). For the latter, the empirical functions usually applied to molecular glass formers as the Cole-Davidson (for the main reorientational mode) or Cole-Cole function (for the secondary relaxations)2629 were used. An excellent description of the experimental data over ten decades of frequency is possible in this way (lines in Fig. 2). All 13 ILs investigated in the present work (see Table 1 for a list) exhibit qualitatively similar behaviour as shown in Fig. 2. In literature, similar spectra are also documented for other ILs7914172022. In the context of the present work, the most interesting quantities determined from these data are σdc and τ. The resistivity values ρdc = 1/σdc obtained by us are listed in Table 1, also including literature data for various additional ILs.

Figure 3 shows an Arrhenius representation of the temperature dependences of σdc and τ for the 13 ILs investigated by us. The found non-linear behaviour reveals different degrees of deviations from the Arrhenius law. From τ the glass temperature can be determined using τ (T) = 100 s (sometimes 250 or 1000 s are assumed, leading to a shift of T of few K only). The obtained values are listed in Table 1. It should be noted that this defines the glass temperature of the ionic subsystem. It should agree with the “structural” glass temperature defined, e.g., via the viscosity if charge transport and viscosity are closely coupled. In any case, it is just this glass temperature which determines the conductivity at room temperature (Fig. 1). It reasonably matches the glass temperature read off from the σdc plot of Fig. 3a using the above-mentioned value at T of 10−15 Ω−1 cm−1 (ref. 23).

Figure 3

Temperature dependence of ionic dynamics.

Data are shown for 13 ILs using an Arrhenius representation. Frame (a) presents the dc conductivity determined from the dielectric spectra. In (b) the conductivity relaxation time deduced from M″(ν) is provided. The solid lines in (a,b) are fits with the VFT formulae, equations (1, 2), respectively. The dashed lines indicate room temperature. For the meaning of the sample abbreviations in the legend, see Table 1.

The data shown in Fig. 3a were fitted with equation (1) (solid lines). The τ data (Fig. 3b) were fitted with the corresponding VFT formula for the relaxation time:

At high temperatures, more data points are available for σdc than for τ because the modulus peaks start shifting out of the available frequency window (Fig. 2d) and the peak frequency cannot be unequivocally determined. Therefore we used D as obtained from the fits of σdc to calculate m. In some cases the phenomenological VFT fits show significant deviations from the experimental data at low temperatures (e.g., for Bmim TFSI). In such cases we used the original definition of m by the slope in the Angell plot for its determination. The deduced values are provided in Table 1, together with data determined from literature.

Finally, Fig. 4a shows the dependence of the dc resistivity at 300 K on m and T for all ILs investigated in the present work and for those from literature. It reveals a clear trend of ρdc to become strongly reduced at high values of m and low values of T. The filled black circles are the projections to the two vertical planes of this 3D plot. Obviously, considering the dependence on T or m alone, without taking care of the other parameter, is not sufficient to account for the observed resistivity variation. This becomes especially obvious when considering cases where substantial variations of ρdc show up, despite the corresponding T or m values are nearly identical (see, e.g., variation for m ≈ 145 at the ρdc-m plane or for T ≈ 180 K at the ρdc-T plane).

Figure 4

Correlation of the room-temperature dc resistivity of ionic liquids with the glass temperature and fragility.

In (a) experimental data for the 13 ionic liquids measured by us and for another 22 compounds taken from literature are included (see Table 1). The spheres are colour coded as indicated in the colour bar. Frame (b) shows a colour-coded plane calculated from the VFT law making a number of assumptions as noted in the text. The colours shown in the m-T plane of frame (a) correspond to the colours of the plane in (b). The colours of the columns, connecting the data points in (a) to the m-T plane, approximately match the colours of this plane. Thus, the calculation result shown in (b) roughly accounts for the experimentally observed absolute values of ρdc in (a).

Figure 4b shows the dependence of ρdc(300 K) on m and T as predicted on the basis of the VFT law. To calculate this colour-coded plane, we used eq. (2) with D replaced by 590/(m-16) (ref. 12). When assuming τ(T) = 100 s and τ0 = 10−14 s, the latter being a typical inverse attempt frequency, one arrives at:

From this the dc resistivity was calculated via ρdc = τ/(ε∞ε0) (ref. 24). Here ε0 is the permittivity of vacuum and ε∞ is the high-frequency limit of the dielectric constant. We used ε∞ = 3, which is a reasonable average value for ILs. The obtained plane shown in Fig. 4b nicely reproduces the experimentally observed trend in Fig. 4a. The colours shown on the m-T plane of Fig. 4a correspond to those of the calculated plane plotted in Fig. 4b. Obviously, this calculation can astonishingly well account even for the absolute values of the experimentally observed ρdc, despite the various simplifying assumptions made and some uncertainties in the determination of m for data sets that do not extend down to T.

Notably, the investigated ILs also include one example containing additional lithium ions [1,3-Dimethylimidazolium bis(trifluoromethane)sulfonimide with 0.1 mol Li]. Such systems are highly relevant for applications in batteries and obviously the correlation of the resistivity with T and m also holds for this system. However, more investigations are needed to check the universality of this approach for such ternary systems. Finally, it should be noted that the ILs considered in this work are non-protic. Protic systems are known to show a stronger variation of ionicity16. Thus it would be interesting to check if the found correlation holds anyhow.

Summary

In summary, the experimental data on 35 ILs, 13 of them investigated by dielectric spectroscopy in the present work, reveal a distinct dependence of the room-temperature resistivity on both glass temperature and fragility. We find that those ILs that combine high fragility with low glass temperature have the highest conductivity. Our work clearly demonstrates that fragility is just as essential for their optimisation for application as the glass temperature. Fragility is an old concept in glass physics. There are numerous approaches trying to explain what makes a glass former fragile, e.g., via energy-landscape variations or an increase of cooperativity1230. For ILs, the dependences of fragility on ionicity16 or anion size10 have been considered. Following these lines seems a promising task on the way to the development of better ILs, suitable for electrochemical applications.

Methods

Sample preparation

The samples were purchased from IoLiTec (Ionic Liquids Technologies GmbH, Heilbronn, Germany) with a minimum purity of 97%. To minimize water content, all samples were dried in N2-gas or vacuum at elevated temperatures for several hours right before measurement. The sample of [1,3-Dimethylimidazolium + Li 1.0 mol/kg] bis(trifluoromethane)sulfonamide was prepared by Dr. Xiao-Guang Sun from Oak Ridge National Laboratory and had a purity >98%.

Dielectric measurements

Different experimental techniques were combined to measure the dielectric properties over a broad frequency range of up to ten decades (0.1 Hz < v < 3 GHz). In the low-frequency range, v < 1 MHz, a frequency-response analyser (Novocontrol Alpha-analyser) and an autobalance bridge (Agilent 4980A) were used. Measurements at v > 1 MHz were performed by a I–V technique where the sample capacitor is mounted at the end of a coaxial line, bridging inner and outer conductor26. For these measurements, impedance analysers (Agilent E4991A or Hewlett-Packard HP4291A) were used. For both methods, the sample materials were filled into parallel-plate capacitors. For cooling and heating of the samples, a N2-gas cryostat was used.

Additional Information

How to cite this article: Sippel, P. et al. Importance of liquid fragility for energy applications of ionic liquids. Sci. Rep. 5, 13922; doi: 10.1038/srep13922 (2015).