Nyquist Geometry and a Two-Stage Optimizer Slash Battery Model Errors by Over 90 Percent

Lithium-ion batteries quietly power nearly every electric vehicle on the road and buffer gigawatts of renewable electricity, yet the software that keeps them safe still depends on mathematical models that are, in many respects, crude approximations of the electrochemical reality inside the cell. Now, researchers Shun-Chung Wang and Yi-Hua Liu report in Results in Engineering a parameter identification framework that could make those models dramatically more faithful. By reading the geometry of impedance spectra the way a cardiologist reads an ECG, and then letting a marine-predator-inspired algorithm refine the picture in two carefully staged passes, their method cuts fitting errors by more than 90 percent compared with conventional approaches while slashing the computational budget by roughly 80 percent in iteration count. The work targets a long-standing bottleneck in battery management systems: extracting accurate model parameters from noisy, high-dimensional data without getting lost in a landscape riddled with false optima.

The heart of the problem lies in how batteries are modeled. Traditional equivalent circuit models represent a cell as a network of resistors and capacitors, which is computationally cheap and therefore attractive for embedded battery management systems, but integer-order components carry hidden assumptions. Real batteries do not behave like ideal capacitors; their charge transfer and diffusion processes are distributed across time scales and across the physical volume of the electrodes, producing frequency-dependent impedance that lumped elements simply cannot replicate. The result is systematic voltage error, poor performance under varying temperatures, and the need for constant recalibration as the cell ages. Fractional-order models solve this by allowing non-integer exponents in the calculus that describes the system, capturing memory effects and non-exponential relaxation naturally. A key element is the constant phase element, whose exponent alpha slides continuously between the behavior of a resistor and a capacitor, and whose value of one half corresponds to the classic Warburg diffusion element seen in electrochemical impedance spectroscopy.

But fractional-order models come with a price: identifying their parameters is brutally difficult. The team adopted an R(RQ)(RQ)W structure, an advanced extension of classical second-order and PNGV models, in which an ohmic resistance is followed by two parallel resistor-constant-phase-element branches and a Warburg element. Each component maps onto a distinct region of the Nyquist plot, the standard graphical representation of impedance data. The ohmic resistance appears as the high-frequency intercept on the real axis, the two RQ branches form overlapping depressed semicircles in the mid-to-high frequency range corresponding to charge transfer and slower interfacial reactions, and the Warburg element produces the characteristic 45-degree low-frequency tail of solid-phase diffusion. Eight parameters must be extracted simultaneously, and they are strongly coupled, meaning an error in one drags the others off course. Metaheuristic optimizers, which explore search spaces without gradients, are the usual tool, but they frequently converge to local optima that look good numerically yet produce physically meaningless parameters.

To choose the best optimizer fairly, the researchers staged a systematic tournament of ten metaheuristic algorithms: three established methods, particle swarm optimization, grey wolf optimization, and artificial bee colony, against seven recent bio-inspired contenders including the sea lion, golden jackal, fox, honey badger, seahorse, sand cat, and marine predator algorithms. Each ran with an identical budget of 10,000 swarm members and 100 iterations across 20 independent trials. Scoring by root mean square error, mean absolute percentage error, and runtime, the Marine Predators Algorithm, which mimics the Lévy flight and Brownian motion strategies of foraging marine animals, emerged as the most robust performer. A Friedman statistical ranking placed it first with an average rank of 2.00, narrowly ahead of the Honey Badger Algorithm, and its mean runtime of 114 seconds comfortably beat both the honey badger’s 150 seconds and the artificial bee colony’s punishing 2,518 seconds. Yet even the winner showed troubling discrepancies for the second semicircle’s resistance and capacitance and the Warburg coefficient, confirming that raw optimization power alone was insufficient.

The first fix is elegantly physical. Instead of letting the optimizer start from arbitrary bounds, the team derived initial parameter estimates directly from the geometry of the Nyquist plot. A parallel resistor-capacitor circuit traces a perfect semicircle of radius R over two; replacing the capacitor with a constant phase element rotates that semicircle by an angle theta that depends analytically on the fractional exponent alpha, following the relation theta equals 90 degrees times one minus alpha. By extending the measured arcs to full circles, reading off their real-axis intercepts, and measuring the rotation angles of their apexes, the researchers could calculate the ohmic resistance, both branch resistances, both constant phase element exponents and capacitance offsets, and even the Warburg value from the lowest-frequency data point. Six of the eight parameters were thus estimated before optimization even began, shrinking the search space from bounds spanning twenty-fold in each direction to a tight window around physically consistent values.

The second fix attacks dimensionality. Rather than asking the optimizer to juggle all eight parameters at once, the two-stage refinement alternates between metaheuristic optimization of six parameters and analytical recovery of the remaining two. Once the first stage converges on R0, R1, Q1, alpha1, alpha2, and W1, the residual impedance, obtained by subtracting those known contributions from the measured spectrum, is fitted by global least squares to extract R2 and Q2 over the second branch’s frequency band. A Trust Region Reflective algorithm then polishes the estimates, and the cycle repeats for ten iterations with bounds that tighten by a factor of two at each stage. The swarm size drops from 10,000 in the brute-force baseline to just 2,000 in the first stage and 100 thereafter, cutting the iteration count by roughly 80 percent while improving precision.

The validation was unusually thorough. On 50 synthetic impedance profiles generated by Latin Hypercube Sampling, where the true parameter values were known exactly, the full framework achieved an overall mean absolute percentage error of just 7.27 percent, compared with a catastrophic 24,462 percent for standard marine predator optimization and an even worse result for two-stage refinement without informed initialization. The mean root mean square error fell to 1.74 times ten to the minus five, a 92.4 percent reduction relative to the baseline, with a median near ten to the minus six. Scatter plots of estimated versus true parameters clustered tightly along the ideal 45-degree line for all eight quantities, while competing methods showed dispersion spanning orders of magnitude. Crucially, the team demonstrated a subtle pitfall: methods 1 and 2 sometimes achieved low numerical error while their simulated Nyquist plots visibly failed to match the target data, proof that a low objective function value does not guarantee a physically accurate model.

Robustness testing added further confidence. Under 1 percent Gaussian noise the overall parameter error rose only from 7.27 to 12.92 percent, and under 2 percent noise to 18.41 percent, with no instability or divergence. Removing the five lowest-frequency data points, simulating instruments that cannot capture the Warburg region, barely changed the error at 7.38 percent, because the initialization relies on mid-to-high-frequency arc geometry. A transfer test on the simpler R(RQ)W structure showed the initialization concept generalizes across model topologies when the correct structure is known, though the authors are careful to note that automatic model-structure selection from raw data remains an open challenge. Against three state-of-the-art published methods, chaotic adaptive fractional particle swarm optimization, hybrid genetic-particle swarm optimization, and complex-order beetle swarm optimization, the new framework delivered curve-fitting errors one to two orders of magnitude lower, with the rivals showing parameter errors often exceeding several hundred percent for the second branch and Warburg elements.

Finally, the method proved itself on real hardware. The team measured electrochemical impedance spectra from commercial LG INR21700M50LT nickel-manganese-cobalt cells in three states: pristine, initially degraded at 90 percent state of health, and significantly aged at 80 percent state of health and 50 percent state of charge. Across 20 independent runs per condition, the full framework reduced mean root mean square error by 38.0, 32.5, and 76.6 percent respectively compared with the standard optimizer, and achieved parameter coefficients of variation as low as 0.016 percent, meaning repeated identifications landed in the same parameter basin with remarkable consistency. That repeatability is what makes the parameters trustworthy inputs for state of charge and state of health estimation, the foundation of reliable battery diagnosis, degradation tracking, and second-life applications. The authors acknowledge their framework assumes the two impedance arcs remain geometrically distinct, and extending it to overlapping arcs and adaptive topology selection is future work. But the central lesson stands: injecting electrochemical knowledge into the starting point of an optimization can accomplish what brute-force computation cannot, at a fraction of the cost.

Subject of Research: Physics-informed parameter identification for fractional-order lithium-ion battery models using electrochemical impedance spectroscopy

Article Title: Physics-informed parameter identification of fractional-order battery models using Nyquist-based initialization and two-stage metaheuristic refinement

Article References: Wang, S.-C., & Liu, Y.-H. (2026). Physics-informed parameter identification of fractional-order battery models using Nyquist-based initialization and two-stage metaheuristic refinement. Results in Engineering, 32, Article 113425. https://doi.org/10.1016/j.rineng.2026.113425

Image Credits: AI Generated

DOI: 10.1016/j.rineng.2026.113425

Keywords: lithium-ion batteries, fractional-order models, electrochemical impedance spectroscopy, parameter identification, metaheuristic optimization, Marine Predators Algorithm, Nyquist plot, battery management systems, constant phase element, state of charge, state of health, equivalent circuit model