Selecting the Regularization Parameter in the Distribution of Relaxation Times

Abstract

Electrochemical impedance spectroscopy (EIS) is used widely in electrochemistry. Obtaining EIS data is simple with modern electrochemical workstations. Yet, analyzing EIS spectra is still a considerable quandary. The distribution of relaxation times (DRT) has emerged as a solution to this challenge. However, DRT deconvolution underlies an ill-posed optimization problem, often solved by ridge regression, whose accuracy strongly depends on the regularization level λ . This article studies the selection of λ using several cross-validation (CV) methods and the L-curve approach. A hierarchical Bayesian DRT (hyper- λ ) deconvolution method is also analyzed, whereby λ 0 , a parameter analogous to λ , is obtained through CV. The analysis of a synthetic dataset suggests that the values of λ selected by generalized and modified generalized CV are the most accurate among those studied. Furthermore, the analysis of synthetic EIS spectra indicates that the hyper- λ approach outperforms optimal ridge regression. Due to its broad scope, this research will foster additional research on the vital topics of hyperparameter selection for DRT deconvolution. This article also provides, through pyDRTtools, an implementation, which will serve as a starting point for future research.