Ionic conductivity calculator from impedance (EIS)

Conductivity from resistance. The electrolyte is sandwiched between two ion-blocking electrodes (stainless steel, Au or Pt) and its resistance R is read from the impedance spectrum. The conductivity is (Front Chem 2021, doi:10.3389/fchem.2020.619832)

σ = L / (R·A)

with L the electrolyte thickness and A the area of the blocking electrode; the same relation, σ′ = (1/R)·(d/A) with d the pellet thickness, gives the bulk conductivity of garnet pellets (J Phys Chem C 2020, doi:10.1021/acs.jpcc.0c03558). With L in cm, R in Ω and A in cm², σ is in S/cm; ×1000 gives mS/cm. A disc electrode of diameter d has A = π·d²/4.

Which resistance. In the Nyquist plot (−Z″ against Z′) of a solid electrolyte, a high-frequency semicircle comes from the bulk (and grain boundaries), and the steep rise at low frequencies from ions piling up at the blocking electrodes (J Phys Chem C 2020; Materials 2021, doi:10.3390/ma14071671). The high-frequency intercept with the real axis is taken as the resistance of the solid electrolyte (ACS Appl Mater Interfaces 2022, doi:10.1021/acsami.2c10853); when one semicircle covers grains and grain boundaries together, its resistance gives the total conductivity (Materials 2021). Each arc is modelled as a resistor in parallel with a constant phase element (CPE), and that resistor gives the conductivity; a CPE exponent n = 1 is an ideal capacitor and n < 1 a depressed arc (ACS Appl Mater Interfaces 2025, doi:10.1021/acsami.4c22390). Arcs with capacitances of a few pF are bulk; a grain-boundary arc in Li garnets had 8.8 nF (J Phys Chem C 2020), and electrode interfaces about 0.1 µF (Chem Mater 2016, doi:10.1021/acs.chemmater.6b00579).

How the tool reads R. High-frequency intercept: from the highest frequency down, the first pair of points where −Z″ changes from ≤ 0 to > 0, with Z′ interpolated linearly at −Z″ = 0 (if −Z″ never crosses zero, the two highest-frequency points are extrapolated, and the tool says so). Semicircle fit: a least-squares circle through the points in the frequency range you choose (the algebraic fit minimising Σ(x² + y² + D·x + E·y + F)²); its two crossings with the real axis are R₁ (high frequency) and R₂ (low frequency), and R₂ − R₁ is the arc’s own resistance. A resistor in parallel with a CPE traces an arc of a circle, so a series resistance R_b plus one such arc gives R₁ = R_b and R₂ = R_b + R_arc exactly; this is checked against synthetic spectra.

Activation energy. Two Arrhenius forms are in use: σ = σ₀·exp(−Ea/kT) (RSC Adv 2022, doi:10.1039/d1ra09138g) and σT = A·exp(−Ea/kT) (Nanomaterials 2021, doi:10.3390/nano11020414), k the Boltzmann constant. The tool fits a straight line to ln(σT) (default) or log₁₀ σ against 1000/T, so Ea = −slope × 1000 × k for ln(σT) and Ea = −slope × 1000 × k × ln10 for log σ, with k = 8.617333 × 10⁻⁵ eV/K (1.380649 × 10⁻²³ J/K ÷ 1.602176634 × 10⁻¹⁹ C). Ea per mole is Ea(eV) × F/1000 = Ea(eV) × 96.485 kJ/mol, with F = 96 485 C/mol (doi:10.1021/acsaem.5c01865). Because ln(σT) = ln σ + ln T, the two forms give Ea values that differ by about k·T (≈ 0.026 eV near room temperature), so report which one you used.

Limits. The circle only describes one arc: leave out the low-frequency electrode spike and inductive points above the axis crossing, and use the arc’s range of frequencies. Which arc is bulk and which is grain boundary has to come from its capacitance or frequency, not from this fit. An Arrhenius line holds only over a range without phase changes: a polymer electrolyte studied from −90 to 100 °C showed three regimes with Ea of 0.35 and 0.45 eV rather than one line (RSC Adv 2022).

The built-in spectrum is illustrative, not measured: it is computed from a series resistance R_b = 120 Ω, an arc of 230 Ω in parallel with a CPE (n = 0.85), a blocking-electrode CPE and 1.5 µH of cable inductance, from 3.16 MHz to 10 Hz. The pellet is 1.05 mm thick between 10 mm electrodes, so A = π × (1.0 cm)² / 4 = 0.7854 cm².

1. Semicircle fit over 10 kHz–1 MHz (13 points, RMS 0.42 Ω off the circle): R₁ = 122.48 Ω, R₂ = 357.04 Ω, arc diameter 234.56 Ω, centre 27.9 Ω below the axis (a depressed arc).
2. Total conductivity from R₂: σ = 0.105 cm / (357.04 Ω × 0.7854 cm²) = 3.744 × 10⁻⁴ S/cm = 0.3744 mS/cm. From R₁ (the series, bulk resistance): 1.092 × 10⁻³ S/cm.
3. High-frequency intercept: −Z″ crosses zero between 1.468 and 1 MHz, at Z′ = 123.81 Ω. The inductance puts it 3 % above the true 120 Ω; the circle fit (122.48 Ω) is closer.
4. Activation energy from five temperatures (25–65 °C, R = 357.0, 246.3, 174.1, 125.8, 92.8 Ω): the ln(σT) line gives Ea = 0.320 eV = 30.88 kJ/mol (R² = 1, A = 2.86 × 10⁴ S K/cm); the same data plotted as log σ give 0.293 eV, the k·T difference between the two forms.

For comparison, Nb-doped Li₇La₃Zr₂O₁₂ pellets measured between Au electrodes from 25 to 50 °C had activation energies of 0.28–0.35 eV (Materials 2021, doi:10.3390/ma14071671), and Al/Ga-substituted LLZO 0.31 eV, within the 0.26–0.37 eV reported before (Chem Mater 2016, doi:10.1021/acs.chemmater.6b00579).