Diffusion impedance

<< Back to: the constant phase element

Semi-infinite diffusion

The simplest and most common circuit element for modelling diffusion behaviour is the Warburg impedance, which models semi-infinite linear diffusion – that is, diffusion in one dimension which is only bounded by a large planar electrode on one side.

Warburg element

The equation for this element is relatively simple:

$$\displaystyle \pmb{Z}_W = \sigma \omega^{-1/2} - j \sigma \omega^{-1/2}$$

where \(\sigma\) is the Warburg coefficient, with units of \(\Omega \text{ s}^{-1/2}\) – I’ll come back to this shortly. In the Nyquist plot, the Warburg impedance gives a straight line with a phase of 45°, which is very recognisable in EIS. If you spot a 45° line in the Nyquist plot, it is usually associated with diffusion.

And if you’re looking at this thinking that this line look very much like a constant phase element but with a phase angle of 45°, then you’d be right – mathematically they are very much the same. In fact, my impedance software of choice (ZView) doesn’t actually provide a semi-infinite Warburg element, so I will ordinarily use a CPE element with n fixed at 0.5. This gives slightly different numbers back when fitting – if you use a CPE, you will get back the Q0 value, rather than the Warburg coefficient \(\sigma\). The two values are related though, by the equation:

$$\displaystyle \sigma = \frac{1}{\sqrt{2} \cdot Q_0} $$

But what does \(\sigma\) actually mean? This is somewhat more complicated, but if you consider the case where you have a soluble redox couple in solution, with a nice reversible electrochemical reaction like:

$$\displaystyle \text{O + n e}^- \rightleftharpoons \text{ R} $$

then \(\sigma\) is related to the diffusion coefficients and concentrations of those species:

$$\sigma = \frac{RT}{n^2 F^2 A \sqrt{2}}\left(\frac{1}{D^{1/2}_O c^\infty_O} + \frac{1}{D^{1/2}_R c^\infty_R}\right) $$

where, in addition to the usual constants, \(D_O\) and \(D_R\) are the diffusion coefficients and \(c^\infty_O\) and \(c^\infty_R\) are the bulk concentrations for the species O and R respectively.

The transmission line

Conventional "Fickian" diffusion is not the only process which gives rise to this type of impedance in electrochemical systems. In batteries, the porosity of electrodes also gives rise to a similarly characteristic 45° line in the Nyquist plot.

This was described in detail by de Levie, who proposed the transmission line model for an electrode with cylindrical pores filled with electrolyte:

transmission line

Consider the two parallel "rails" as being the electronic resistance in the electrode material itself and the ionic resistance in the electrolyte respectively. The capacitors represent the double layer capacitance.

This equivalent circuit, when infinitely long, gives an impedance response which is identical to the Warburg impedance above. Practically, however, porous electrodes have a finite length, and so show a 45° line only in a certain frequency range. The impedance response due to finite diffusion is discussed below.

A final point on the transmission line: if the transmission line shows the same response as Fickian diffusion, can the case of a porous electrode be considered diffusion as well? In a sense, yes. The movement of ions through the pores is coupled to the movement of electrons through the pore walls. This is an example of ambipolar diffusion.

Finite diffusion

Often in the "classic" electrochemical setups, diffusion often appears semi-infinite because the timescale of the experiment is not long enough for the system to reach a steady state. However, in many real systems and in some standard experiments, diffusion is either naturally, or by design, limited. This gives rise to finite diffusion behaviour, which shows a different response than the standard Warburg impedance.

There are two important equivalent circuit elements for finite diffusion. They are the finite length Warburg (FLW) and the finite space Warburg (FSW), sometimes called the "short" and "open" Warburg elements respectively. Their responses in a Nyquist plot look like this:

Let's look at the FLW first. Mathematically, it can be written as:

$$\displaystyle \pmb{Z}_\text{FLW} = Z_0 \left(j \omega \tau \right)^{-1/2} \tanh \left(j \omega \tau \right)^{1/2}$$

The two variables \(Z_0\) and \(\tau\) reflect the properties of the system or the process giving rise to the FLW behaviour. The impedance of the FLW tends to the value of \(Z_0\) at low frequency. At high frequencies, the response is almost exactly that of the Warburg impedance. The shape of the FLW in the Nyquist plot therefore looks like a 45° line at higher frequencies, and transitions into a semi-circle shape at low frequencies. It can be thought of as a Warburg being "shorted" by a resistor - although a parallel W-R circuit will not give the same response.

This response is typically associated with diffusion (or more generally mass transport) through a layer with a finite length. A classic example of this is the response of the rotating disk electrode - the point of which being to reduce the distance from the electrode to the bulk by controlling convection, rather than diffusion in this case.

The FSW has a rather similar definition:

$$\displaystyle \pmb{Z}_\text{FSW} = Z_0 \left(j \omega \tau \right)^{-1/2} \coth \left(j \omega \tau \right)^{1/2}$$

In this case, the response tends towards capacitive-like behaviour at low frequencies, where \(Z'(\omega = \infty) = Z_0/3\). This response is typically associated with diffusion where one of the boundaries is blocking, such as in a porous electrode, as discussed previously on this page. This response is also associated with the diffusion of ions within a storage electrode, such as in lithium-ion batteries.

Be warned! You can find a number of different definitions of these finite elements and it's easy to get confused. One source of confusion can be from defining these elements in terms of admittance. For example, the FLW is frequently defined like this:

$$\displaystyle \pmb{Y}_\text{FLW} = Y_0 (j \omega)^{-1/2} \coth \left[B (j \omega)^{-1/2}\right] $$

You might have noticed that this definition contains \(\coth\) instead of \(\tanh\). In this case it's important to remember that:

$$\coth(x) = \frac{1}{\tanh(x)}$$

Next: Simulating the Randles circuit >>

Thank you for reading!

These pages will remain open to all and 100% ad-free, now and always. If you like this page and want to see more, please consider buying me a cup of coffee, to keep me motivated and to help with the running costs of this site!

Buy Me a Coffee at ko-fi.com
comments powered by Disqus