Non-negatively constrained least squares and parameter choice by the residual periodogram for the inversion of electrochemical impedance spectroscopy data

Abstract

The inverse problem associated with electrochemical impedance spectroscopy requiring the solution of a Fredholm integral equation of the first kind is considered. If the underlying physical model is not clearly determined, the inverse problem needs to be solved using a regularized linear least squares problem that is obtained from the discretization of the integral equation. For this system, it is shown that the model error can be made negligible by a change of variables and by extending the effective range of quadrature. This change of variables serves as a right preconditioner that significantly improves the condition of the system. Still, to obtain feasible solutions the additional constraint of non-negativity is required. Simulations with artificial, but realistic, data demonstrate that the use of non-negatively constrained least squares with a smoothing norm provides higher quality solutions than those obtained without the non-negative constraint. Using higher-order smoothing norms also reduces the error in the solutions. The L-curve and residual periodogram parameter choice criteria, which are used for parameter choice with regularized linear least squares, are successfully adapted to be used for the non-negatively constrained Tikhonov least squares problem. Although these results have been verified within the context of the analysis of electrochemical impedance spectroscopy, there is no reason to suppose that they would not be relevant within the broader framework of solving Fredholm integral equations for other applications.

Analysis of Electrochemical Impedance Spectroscopy on Zinc-Air Batteries Using the Distribution of Relaxation Times

“…The determination of the regularization parameter for the ridge regression is the current state of research and leads to different approaches [39,51,55,63]. Without regularization, the DRT would tend to oscillate, but excessive values for λ can suppress characteristic peaks (see Figure 4).…”

Monitoring and Diagnostics of SOFC Stacks and Systems

“…To overcome this problem, spectral analysis tools have been successfully applied, yielding distribution of relaxation time constants (DRT), in which the different processes appear as distinct peaks. The deconvolution problem is solved by means of Tikhonov regularization (8)(9). The DRT de-convoluted spectrum is used to understand which process is responsible for the degradation or the fault, providing metrics for the diagnostic and lifetime algorithm.…”

Fitting procedure based on Differential Evolution to evaluate impedance parameters of metal–coating systems

“…Even today, the analysis and interpretation of EIS data is an arduous task, though there are several commercial software packages available for the modelling of the experimental data, such as EQUIVCRT developed by Boukamp (1986), ZSimpWin (Yeum, 2001) as well as others included in the software of acquisition of experimental data such as NOVA (Metrohm Autolab, 2015). The adjustment procedures of the experimental data are always based on non-linear regression algorithms such as the Gauss-Newton method, its modified variant and the Levenberg-Marquardt method (Macdonald and Garber, 1977;Macdonald et al, 1982;Macdonald, 1990;Hansen et al, 2015). The resulting fits can be very accurate if a good initial estimation of the adjustment parameters is introduced.…”