Analytical–Computational Integration of Equivalent Circuit Modeling, Hybrid Optimization, and Statistical Validation for Electrochemical Impedance Spectroscopy

Background: Electrochemical impedance spectroscopy (EIS) is indispensable for disentangling charge-transfer, capacitive, and diffusive phenomena, yet reproducible parameter estimation and objective model selection remain unsettled. Methods: We derive closed-form impedances and analytical Jacobians for seven equivalent-circuit models (Randles, CPE, Warburg variants), enforce physical bounds, and fit synthetic spectra with 2.5% and 5.0% Gaussian noise using hybrid optimization (Differential Evolution → Levenberg– Marquardt). Uncertainty is quantified via non-parametric bootstrap; parsimony is assessed with RMSE, AIC, and BIC; physical consistency is checked by Kramers–Kronig diagnostics. Results: Rs and Rct are consistently identifiable across noise levels. CPE parameters (Q, n) and diffusion amplitude (σ) exhibit expected collinearity unless the frequency window excites both processes. Randles suffices for ideal interfaces; Randles+CPE lowers AIC when non-ideality and/or higher noise dominate; adding Warburg reproduces the 45◦ tail and improves likelihood when diffusion is present. The (Rct + ZW ) ∥ CPE architecture offers the best trade-off when heterogeneity and diffusion coexist. Conclusions: The framework unifies analytical derivations, hybrid optimization, and rigorous statistics to deliver traceable, reproducible EIS analysis and clear applicability domains, reducing subjective model choice. All code, data, and settings are released to enable exact reproduction.