On the Electrochemical Impedance Spectroscopy of Marcus-Hush-Chidsey Heterogeneous Kinetics

Abstract

The last decade’s research and development of porous intercalating Li-ion electrodes have resulted in, among other, optimized mass-transport capabilities. This has recently allowed researchers to attribute the reorganization of the lithium ion’s solvation shell, when intercalated, to contribute to a significant overpotential (50-100 mV) at fast charging rates1,2. A further development of analytical techniques is therefore required to understand and quantify the solvent reorganization energy related to the power losses observed at higher current densities. Kinetic theories that describe heterogeneous charge-transfer where electrons are transferred from an electrode to a reactant in solution without adsorption and breaking of bonds has been intensively studied over the past century. Today three main theories are used to describe aprotic electron transfer: Butler-Volmer (BV), Marcus-Hush (MH), and Marcus-Hush-Chidsey (MHC). The Butler-Volmer theory3,4 is widely recognized and applied to electrochemical data to explain the evolution of rate constants as a function of overpotential. However, the BV theory provides little physical insight into understanding the limitations of the electron transfer whereas the Marcus-Hush theory has been able to predict the evolution of homogeneous rate constants as function of activation energy decades before they were experimentally confirmed - especially for the so-called inverted region5–10. The Marcus-Hush theory includes the reorganization energy (λ), the parameter causing the inverted region, and accounts for the energy change in the solvation shell of the redox active molecule. The Marcus-Hush-Chidsey theory is based on the original work by Marcus and Hush and was combined with the Fermi-Dirac distribution describing the probability of electrons around the Fermi level in a conducting electrode11,12. MHC theory therefore provides a microscopic description of the heterogeneous rate constants that, among other, captures the reorganization energy during an electron transfer and depending on the reorganization energy forces the heterogeneous rate constants to significant deviate from the classical BV theory as shown in figure 1. However, extracting kinetic contributions from any electrochemical system through voltammetry has previously proven to be difficult, as fast kinetics and mass-transport often complicate matters. Electrochemical impedance spectroscopy (EIS) is a technique that relies on a small periodic ac amplitude around a dc input. The ac input can be measured over a wide frequency range allowing for the careful separation of kinetic-, double-layer, and mass-transport contribution as each process often occurs at different time scales. Kinetic contributions are referred to as charge-transfer resistances (RCT), double-layer as capacitors and their derivatives, while mass-transport have a wide range of terms. Herein, we present formulas for the charge-transfer resistance that is governed by infinite or finite MHC kinetics and compare these to the classical BV representations. This allows for the numerical simulation and fitting of the EIS response for any electrochemical system that is governed by MHC kinetics. In addition, we have also solved the RCT governed by the simplified MHC terms developed by Zeng et al.13, which simplify the otherwise complex Fermi-Dirac integrals. Consequently, this opens the possibility of predicting and determining the reorganization energy by EIS.