Anirban Roya,
Bula Duttab and
Subhratanu Bhattacharya*a
aDepartment of Physics, University of Kalyani, Kalyani, Nadia-741235, West Bengal, India. E-mail: subhratanu1@gmail.com
bDepartment of Chemistry, P. D. Women's College, Jalpaiguri-735101, West Bengal, India
First published on 5th July 2016
Herein, a physical model, based on impedance spectroscopy and space charge polarization with the consideration of low frequency capacitance dispersion is presented to evaluate the parameters describing the electrical transport properties of some previously studied polymer electrolytes. Implementing the model, the complete frequency response of complex conductivity within the measured frequency range can be imitated, which enables us to appropriately evaluate the macroscopic DC conductivity, hopping frequency and double layer frequency within the measured temperature range. The temperature dependent mobile ion concentration, free-ion diffusivity and average ion hopping length have also been estimated from the analysis of the frequency dependent real and imaginary parts of the complex permittivity spectra, using the same model. Reasonable agreement of the ion diffusivity obtained from the present analysis to that obtained previously from the pulsed-field gradient (PFG) NMR measurements strongly justifies the applicability of the model to a wide variety of ion conducting systems.
The fundamental problem inhibiting the understanding of ionic conductors is the incapability to explicitly estimate the free-ion diffusivity and mobile ion concentration from the measured dc conductivity. Moreover, it is also difficult to unambiguously determine the macroscopic dc conductivity of the ionic conductors in presence of the above mentioned EP effect in the experimental ac impedance data. Usually, the dc conductivity has been determined from the narrow frequency-independent region of the real part of the frequency dependent conductivity or from the linear part of the imaginary complex permittivity. But, both the processes can give erroneous results as the contribution of overlapping electrode polarization is generally not considered in these methods. Evaluation of free-ion diffusivity and concentration of mobile ions from the analysis of EP has been modeled by Coelho,4 Schutt,5 and Trukhan.9 Although, all these models do produce reasonable estimates for diffusivity, mobile ion concentration and diffusion coefficients in a few cases, it fails for most of the materials.10,11 For systems with high ion concentration, the above models overvalues the free-ion diffusivity while undervalues the mobile-ion concentration with respect to that obtained from NMR.10,11 These may be attributed to the oversimplified consideration of smooth electrode–electrolyte interfaces5,7 and the method of determination of dc conductivity values used in these models.
In the current study, a physical model is presented that enables us to simulate the complete frequency response of complex conductivity within the measured frequency range to unambiguously determine the macroscopic dc conductivity, mobile-ion diffusivity and the concentration of mobile ions of some lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in poly(ethylene oxide) (PEO) polymer electrolyte system over a wide range of EO/Li ratio which had been studied earlier by PFG-NMR measurement.12 The rough electrode–electrolyte interface has taken into account by considering the frequency dependent complex phase element instead of frequency independent double layer capacitance.10,13 The ion diffusivity evaluated according to the proposed model shows reasonable agreement to the reported PFG-NMR value.12 The average hopping length of mobile ions, evaluated from the proposed model has been correlated to the ion conductivity and diffusivity of the materials to get better insight to the general understanding of the conduction mechanism within the ion conducting materials.
The admittance of the sample can be evaluated from the equivalent electrical circuit shown in Fig. 1a which follows the classical PNP model. Here the bulk conductivity σb = (d/A)(1/Rb) (conductivity at much higher frequency than the region where the electrode polarization phenomena dominates i.e. at the frequency independent region of the conductivity spectra), is connected in series with the frequency independent double layer interface capacitance Cdl = Aε0εs/2LD and this series combination is in parallel with the bulk capacitance Cdl = Aε0εs/d where Rb is the bulk resistance and εs the dielectric constant of the electrolyte. Here d and LD are sample thickness and the Debye length (the length scale of the electrostatic double layer) respectively, satisfying the condition dLD−1 ≫ 1 and A is the surface area of the blocking electrodes. The total number density of free ions (Nc), which is the sum of the number densities of free cations (n+) and anions (n−), can be related the Debye length by the following equation:9
![]() | (1) |
![]() | (2) |
It is worth evident from the experimental data that the electrode polarization is frequency dependent i.e. show capacitance dispersion with frequency. To explain the observed behaviour, the double layer capacitance must be considered as frequency dependent, which also replaces the simplified assumption of plane electrode–electrolyte interfaces5,13 in circuit Fig. 1a by the interface with roughness and pores. The model can be demonstrated by the equivalent circuit Fig. 1b, in which the double layer capacitance term has been replaced by the constant phase element (CPE) term with impedance ZCPE(ω) = Cdl−1τ∞(iωτ∞)−α with 0 < α ≤ 1 and the Maxwell–Wagner relaxation time τ∞ = RbCb = ε0εs/σb = ε0εs/Ncqμ to obtain the effective complex conductivity as
![]() | (3) |
Considering eqn (3) with ZCPE, the expression for the frequency dependent real and imaginary part of effective conductivity with the contribution of electrode polarization become
![]() | (4) |
![]() | (5) |
This is the expression for the frequency dependent conductivity up to the high frequency limit as the onset of electrode polarization where σ′eff ∼ σb ∼ σdc.
When the applied frequency becomes higher than a characteristic hopping frequency ωH, the free charge carriers may behave as dipoles with non-Debye-like absorption, the conductivity starts to disperse with frequency, characteristic of the “universal dielectric response”.11,14 The frequency dispersion of the conductivity occurs at about the frequency ω ≈ ωH. At higher frequency ω > ωH both the real and imaginary part of the ac conductivity increase with frequency.15 Thus, above the onset frequency of electrode polarization, it is illustrative to express the frequency dependence of the real part of the conductivity, σ′(ω) in terms of the dc conductivity, σb and hopping frequency ωH, using the Jonscher power law15 and Almond–West formalism16 as
| σ′(ω) = σb[1 + (ω/ωH)n] | (6) |
Similarly, the linear frequency dependence of the imaginary part of conductivity can be expressed as15
| σ′′(ω) = Aωs | (7) |
![]() | (8) |
Considering the equivalent circuit of Fig. 1a the complex dielectric function of electrode polarization is typically demonstrated by a Debye type relaxation as4,5,7,13
![]() | (9) |
In eqn (9), Δε = εs,EP − εs is the dielectric strength and εs,EP = εsd/2LD = εsM is the low frequency dielectric constant in presence of EP. The characteristic double layer relaxation time of electrode polarization is given by
![]() | (10) |
Considering the equivalent circuit in Fig. 1b and expressing ZCPE in terms of τ∞ the expression for the complex permittivity becomes7
![]() | (11) |
The real part of eqn (11) is
![]() | (12) |
![]() | (13) |
It is important to notice that the above expression should give the frequency dependence of the complex permittivity up to the onset frequency of electrode polarisation, above which the observed frequency variation should be explained by high frequency non-Debye type dielectric relaxation phenomena6 of the electrolytes which are not considered here.
In the case of α = 1, the CPE reduces to a double layer capacitance (eqn (9)) and the expressions for real and imaginary part of the complex permittivity, considering the circuit of Fig. 1a, reduces to respectively.
![]() | (14) |
For the symmetric electrolytes, from the Nernst–Planck expression and from the definition of LD, the bulk conductivity can be expressed as9
| σb = [εs + ΔεEP]ε0D/LD2 | (15) |
Following the above equations the ion diffusion coefficient D can be expressed as
![]() | (16) |
From the microscopic description of diffusion and conduction proposed by Einstein and Smoluchowski17 for the ion conducting systems we get
![]() | (17) |
| 〈λ2〉 = 6D/ωH | (18) |
Polymer electrolytes, comprised of different poly(ethylene oxide) (PEO) (Sigma-Aldrich, MW ∼ 9 × 105) and lithium trifluoromethanesulfonylimide (LiTFSI) (Alfa-Aesar) were prepared by solution casting method in a vacuum dry box, according to the procedure described in the report.12 Four polymer electrolytes, corresponding to EO/Li ratios 6, 8, 10 and 20 respectively were prepared by dissolving the polymer and the lithium salt (LiTFSI) in acetonitrile (Sigma-Aldrich) as and subsequently removing the solvent in a vacuum oven at 60 °C for two days.
![]() | ||
| Fig. 2 (a and b) Temperature and (c and (d) compositional variation of frequency dependent real and imaginary components of the complex conductivity respectively. Different regions corresponding to the electrical response due to electrode polarization and conductivity relaxation are indicated for the guide of eye. Solid lines are best fits to the combination of eqn (4) with eqn (6) and eqn (5) with eqn (7) respectively. | ||
It is observed from the figures that the frequency dependence of ion conduction within the materials can be characterized by three mutually correlated phenomena viz. low frequency electrode polarization, intermediate frequency independent macroscopic charge conduction and high frequency conductivity dispersion. The most striking features of electrode polarization are manifested in the frequency variation of complex conductivity are emphasized by the arrows in Fig. 2b.
The phenomena come into play at a certain onset frequency (ωon), where a minimum in σ′′(ω) can be observed first. Below the onset frequency (ωon), σ′(ω) slowly decreases with lowering the frequency. At lower frequencies with the full development of the electrode polarization near the frequency ωmax, a peak in σ′′(ω) starts to develop (Fig. 2c and d) with a simultaneous sharp decrease in σ′(ω) (Fig. 2a and b). Above the onset frequency of electrode polarizations σ′(ω) remains frequency independent up-to a certain frequency, which, in accordance with other reports19–21 can be ascribed to the macroscopic dc conductivity of the materials. With further increase of frequency the slope of the conductivities change in frequency and σ′(ω) approaches power law behaviour. The crossover from the frequency independent conductivity to the conductivity dispersion indicates the commencement of the conductivity relaxation phenomena.19,20
It is observed from Fig. 2 that as the temperature decreases, the overall conductivity decreases and the onset of electrode polarization (ωon) shifts toward lower frequencies with a reduction in the frequency range of electrode polarization. Simultaneously, the onset of the conductivity dispersion from frequency independent region also shifts towards lower frequency. This indicates that both the phenomena are strongly correlated with each other. It is interesting to observe from the real and imaginary components of the frequency dependent complex conductivity have been analysed by the combination of eqn (4) with eqn (6) and eqn (5) with eqn (7) respectively. The absolute agreement of the simulated data by the combination of the above equations with the experimental data justifies the applicability of the model to analyse the frequency dependence of the total conductivity within the measured frequency and temperature range. In this fitting, the value of ‘s’ remains close to 1 for all the measured temperatures and the frequency exponent ‘n’ remains almost temperature independent.
The temperature dependence CDR for an electrolyte is shown in Fig. 3a. It is worthy to observe that for each of the temperatures in contrast to the frequency independent double layer capacitance Cdl considered in the PNP model (Fig. 1a), CDR gradually decreases with frequency and attains a frequency independent constant value after the nearby frequency of ωon. Also the magnitude of CDR at the lower frequency gradually decreases with the decrease of temperature.
![]() | ||
| Fig. 3 (a) Temperature and (b) composition dependency of the real part of constant phase element (CDR). | ||
The temperature dependence of the Maxwell–Wagner relaxation time τ∞ obtained from the above fitting are depicted in Fig. 4a. The decreasing value of τ∞ with increasing temperature clearly demonstrates that the electrode effect dominates at higher temperature. Furthermore.
The temperature dependence of the macroscopic dc conductivity σdc and hopping frequency ωH obtained from the above fitting are shown in Fig. 5a and b respectively. It is worthwhile to mention that the variation of the conductivity with respect to the EO/Li ration show similar trend to that observed previously.12,22,23 It is observed that all the samples exhibit a non-Arrhenius temperature dependence of the dc conductivity within the measured temperature range, similar to the previously studied identical systems.24
![]() | ||
| Fig. 5 Temperature variations of (a) dc conductivity and (b) hopping frequency for different electrolytes. Solid lines are best fits to VTF equations (eqn (18)). | ||
This non-Arrhenius temperature dependence of the conductivity has been fitted by the modified Vogel–Tammann–Fulcher (VTF) equation expressed as20,25,26
σdc = σ0 exp(−BT0/(T − T0))
| (19) |
| Samples | Dc conductivity | Hopping frequency | Ion diffusivity | ||||||
|---|---|---|---|---|---|---|---|---|---|
log σ0 |
T0 (K) | B | log ωH0 |
T0 (K) | B | log D0 |
T0 (K) | B | |
| (PEO)6 (LiTFSI) | 0.27 | 177 | 6.34 | 10.23 | 176 | 6.73 | −4.22 | 171 | 6.63 |
| (PEO)8 (LiTFSI) | 0.37 | 186 | 5.70 | 10.57 | 189 | 5.87 | −4.39 | 181 | 6.07 |
| (PEO)10 (LiTFSI) | 0.43 | 162 | 5.28 | 10.51 | 164 | 5.20 | −4.54 | 186 | 4.51 |
| (PEO)20 (LiTFSI) | 0.38 | 177 | 6.00 | 10.45 | 175 | 5.99 | −4.13 | 169 | 5.87 |
It can as well be mentioned that the hopping frequency ωH for all the electrolytes follows VFT temperature dependence, and the Vogel temperature and strength parameter for ωH closely follow those of the dc conductivity process (Table 1) for the electrolytes with different EO/Li ratio, confirming that for all the electrolytes there is a strong coupling between the macroscopic ion transport, ion pair dissociation and the hopping frequency.
The observed temperature independence of the power law exponents for these electrolyte systems prompts us to analyse the scaling behaviour of the conductivity spectra of the systems, considering the hopping frequency and dc conductivity as the scaling parameters as reported earlier.27 Fig. 6 depicts the scaling of the real and imaginary part of the conductivity spectra for an electrolyte at different temperatures respectively. The perfect overlap of both the conductivity spectra for different temperatures justifies the selection of the parameters. Similar nature of scaling for all the electrolytes with different EO/Li ratio clearly indicates that the physical mechanism involving the ion dynamics within the electrolytes is temperature independent.
![]() | ||
| Fig. 6 Temperature scaling of the real and imaginary part of the conductivity spectra for an electrolyte. | ||
![]() | ||
| Fig. 7 Frequency dependence of (a) and (c) real and (b) and (d) imaginary part of the permittivity at different temperatures for an electrolyte and at a particular temperature for different electrolytes. Solid lines are best fits to eqn (12) and (13) respectively. | ||
However, the fitting of both the components start to deviate at higher frequency and at lower temperature, where the dielectric relaxations of the Li ions coordinated polymer segments become prominent and the above equations cannot explain the same. Typically the empirical equation of Havriliak and Negami (HN) may be well fitted to experimental data for the polymer relaxations at the high frequency region6 which beyond the scope of this paper.
Plots of temperature dependent relaxation times for electrode polarization τEP obtained from the above fittings are presented in Fig. 8. As observed, τEP shows an identical temperature dependency to that of τ∞ for each of the electrolytes, but an opposite composition dependency. Moreover, the value of τEP is found to be much higher than that of τ∞ at each temperature for all the electrolytes.
![]() | ||
| Fig. 8 Temperature variation of relaxation times for electrode polarization τEP for different electrolytes. | ||
The ion diffusivity (D) at different temperatures for different electrolytes have been estimated from eqn (16) and are shown in Fig. 9 in an Arrhenius fashion. The average of the absolute translational diffusion coefficients Davg = DLi7 + DF19/2 estimated from the PFG-NMR measurement by Gorecki et al.12 has been incorporated in the figure for comparison. The diffusion coefficient from EP analysis is found to be slightly higher than the values obtained from NMR experiments, but shows similar trend with respect to EO/Li ratio. Also, the ‘D’ value obtained from the current analysis is much lower than that was estimated using Debye model for identical electrolyte system by Wang et al.10 Nevertheless, it is worthy to mention that the PFG-NMR measurement evaluates the average diffusivities of both free ions and ion pairs, whereas the EP phenomenon in theory only reflects the diffusivity of free ions. Also there is a difference in the experimental and other allied conditions. Considering the above aspects, the estimated value of ‘D’ from the current model is quite reasonable as compared to the PFG-NMR value.12 The observed non-linearity with temperature in contrast to the PFG-NMR data may be attributed to the difference in the measured temperature range.
![]() | ||
| Fig. 9 Temperature variation of ion diffusivity for different electrolytes determined from electrode polarization and PFG-NMR.12 Solid lines are best fits to VTF equation. | ||
Considering the relation between the conductivity (σb), mobile ion concentration Nc and the free ion diffusivity the temperature dependent mobile ion concentration of different electrolytes has been evaluated as shown in Fig. 10a. It is observed that with the increase of temperature the mobile ion concentration linearly increases for all the electrolytes.
The nonlinear variation of the data are well described by similar VFT equation as that of eqn (19). The VFT temperature dependence of diffusivity reflects the coupling of segmental motion and ion transport. Fitting parameters are included in Table 1. The estimated value of the free ion concentration is in agreement with the general expectations for polymer electrolytes, as reported earlier.10
Yet, for the electrolyte with EO/Li ratio 20, the value of Nc shows minimum suggesting incomplete dissociation of the salt within the polymer matrix. The temperature dependence of Nc is well described by an Arrhenius equation Nc = Ntot
exp(−Ea/kBT) (solid lines) with Ntot is the total ion concentration at full salt dissociation and Ea is the activation energy of the free ion concentration, signifies the energy needed to separate an ion pair (binding energy). The Arrhenius temperature dependence of ion concentration implies that the ion becomes unpaired (truly free) through a thermally activated process of fully overcoming electrostatic attraction. The concentration of unpaired ions will therefore increase with increasing temperature and the ion association seems therefore to be ridden by the electrostatic attraction between anion and cation. The activation energy for the electrolytes with EO/Li ratio 6, 8, 10 and 20 are found to be 0.1, 0.07, 0.05 and 0.13 eV respectively. Thus, for the electrolyte with EO/Li ratio 20 the ion dissociation requires much larger energy than that of the electrolyte with EO/Li ratio 10 (0.05 eV).
The above features of ion concentration are in accordance with earlier investigations for different electrolyte samples based upon the analysis of electrode polarization.1,26,28
In this regard, it is worthwhile to mention that a correction procedure in the estimation of diffusivity and mobile ion concentration has also been proposed by different researchers10,11,26 to accommodate the electrode roughness into the electrode polarization analysis. However, the applicability of the above correction procedure is under debate10 as that approach cannot extend to the reasonable estimation of the electrical parameters for polymer electrolytes where, the dissociation energy would change significantly with salt concentration.10 The current model pre-considers all the concerns related to the electrode roughness and hence the values of parameters estimated from the analysis are physically reasonable and far more accurate.
From eqn (18), the average hopping length of mobile ions
within the electrolytes has been evaluated as depicted in Fig. 10b with temperature. It is evidenced from Fig. 10b that 〈λ2〉 sharply decreases with increase of temperature for each of the electrolytes i.e. it follows reciprocal trend to that of Nc. For the electrolyte with EO/Li ratio 20,
is much higher than that of the electrolyte with EO/Li ratio 10. It has been reported earlier26 that the interactions between the ions and polymer backbones transforms the crystalline regions of PEO into amorphous regions by distorting or stretching its helical conformation and local chain organization. Furthermore, the phase diagram of the (PEO)m–(LiTfSi) system studied earlier showed that the electrolytes with m < 10 are more amorphous that that of the electrolytes with m > 10.23 Thus, the observed difference in the crystallinity for the above electrolytes strongly affects to the ion dissociation and formation of hopping cites within the electrolytes.12,23,26
As the amorphous region progressively increases at higher temperature; the polymer chains obtain quicker internal modes in which bond rotations enhance segmental motion. As a consequence, the inter-chain and intra-chain hopping movement of ions is favoured and the ions can hop with higher frequency to avail higher conductivity with shorter average hopping length. To examine this correlation of segmental motion and hopping frequency directly we plot both dc conductivity and free ion diffusivity against the frequency of segmental motion i.e. the characteristic hopping frequency ωH as shown in Fig. 11a and b respectively. For ion motion strongly coupled to polymer segmental motion, conductivity is expected to obey the Debye–Stokes–Einstein (DSE) equation26 σb ∝ ωHn. The conductivity for all the electrolytes behaves according to the ideal DSE equation with a slope nearly equal to 1 in the log–log plot (Fig. 11a, implying that the conductivity is strongly regulated by the segmental motion, i.e. the hopping frequency of the ionic clusters.
![]() | ||
| Fig. 11 (a) Dc conductivity and (b) free ion diffusivity against hopping frequency. Solid lines are best fits to DSE equations. | ||
Similar relation may be observed for the diffusivity also26 as shown in Fig. 11b. However, the ion diffusivity follows a fractional DSE type behavior with n < 1 as shown in Fig. 11b. The deviation from ideal DSE behavior for the ion diffusivity than for the conductivity can be attributed to the considerable temperature dependency of the mobile ion content, which is more significant for the electrolyte with EO/Li = 6 and 20 respectively.
| This journal is © The Royal Society of Chemistry 2016 |