Paper deep dive
Physics-informed distribution of relaxation times estimation and latent-space condition monitoring of solid oxide fuel and electrolysis cells from electrochemical impedance spectroscopy
Žan Gorenc, Žiga Gradišar, Felix Mütter, Vanja Subotić, Pavle Boškoski
Intelligence
Status: succeeded | Model: Gemma-4-26B-A4B | Prompt: intel-v1 | Confidence: 94%
Last extracted: 8/16/2026, 2:14:45 AM
Summary
This paper proposes a physics-informed convolutional autoencoder to estimate the Distribution of Relaxation Times (DRT) from Electrochemical Impedance Spectroscopy (EIS) data. By embedding the physical relationship between impedance and DRT into the training process, the model resolves overlapping relaxation processes and accurately reconstructs impedance spectra across multiple solid oxide fuel and electrolysis cell datasets with errors below 1.1%. The learned latent space is organized by relaxation timescale and effectively captures operating changes, hydrogen-shortage events, and long-term degradation, providing a robust basis for condition monitoring.
Entities (8)
Relation Signals (8)
Physics-informed convolutional autoencoder → estimates → Distribution of Relaxation Times
confidence 98% · We propose a physics-informed convolutional autoencoder that estimates DRT directly from EIS data
Physics-informed convolutional autoencoder → usesinput → Electrochemical Impedance Spectroscopy
confidence 98% · estimates DRT directly from EIS data
Physics-informed convolutional autoencoder → appliedto → Solid Oxide Fuel Cell
confidence 95% · accurately reconstructs measurements from three independent solid oxide fuel and electrolysis cell datasets
Physics-informed convolutional autoencoder → appliedto → Solid Oxide Electrolysis Cell
confidence 95% · accurately reconstructs measurements from three independent solid oxide fuel and electrolysis cell datasets
Latent Space → captures → long-term degradation
confidence 92% · Distances in this latent space capture operating changes, hydrogen-shortage events, and long-term degradation.
Latent Space → captures → hydrogen-shortage events
confidence 92% · Distances in this latent space capture operating changes, hydrogen-shortage events, and long-term degradation.
Jožef Stefan Institute → affiliatedwith → Žan Gorenc
confidence 90% · Affiliation: Jožef Stefan Institute... Žan Gorenc
→ →
Cypher Suggestions (0)
No Cypher suggestions yet.
Abstract
Abstract:Estimating the distribution of relaxation times (DRT) fromelectrochemical impedance spectroscopy (EIS) is an ill-posed inverse problem that is highly sensitive to regularisation choices. We propose a physics-informed convolutional autoencoder that estimates DRT directly from EIS data without spectrum-specific tuning. A discretised relation between impedance and the DRT is embedded in the training process, constraining the network to produce impedance-consistent distributions. The model resolves overlapping relaxation processes in synthetic two-ZARC spectra and accurately reconstructs measurements from three independent solid oxide fuel and electrolysis cell datasets, with range-normalised errors below 1.1%. Decoder-probe analysis shows that the learned latent representation is organised according to relaxation timescale. Distances in this latent space capture operating changes, hydrogen-shortage events, and long-term degradation. The same lightweight architecture is applied across all datasets without modification, providing consistent DRT estimation and an interpretable basis for condition monitoring.
Tags
Links
- Source: https://arxiv.org/abs/2608.13305v1
- Canonical: https://arxiv.org/abs/2608.13305v1
Trouble viewing inline? Open PDF directly →
Full Text
62,183 characters extracted from source content.
Expand or collapse full text
Physics-informed distribution of relaxation times estimation and latent-space condition monitoring of solid oxide fuel and electrolysis cells from electrochemical impedance spectroscopyJournal: International Journal of Hydrogen Energy Žan Gorenc Email: zan.gorenc@ijs.si Corresponding author: Corresponding author Affiliation: Jožef Stefan Institute, Jamova cesta 39, Ljubljana, 1000, Slovenia Affiliation: Jožef Stefan International Postgraduate School, Jamova cesta 39, Ljubljana, 1000, Slovenia Žiga Gradišar Email: ziga.gradisar@ijs.si Affiliation: Jožef Stefan Institute, Jamova cesta 39, Ljubljana, 1000, Slovenia Affiliation: Jožef Stefan International Postgraduate School, Jamova cesta 39, Ljubljana, 1000, Slovenia Felix Mütter Email: felix.muetter@tugraz.at Affiliation: Institute of Thermal Engineering, Graz University of Technology, Inffeldgasse 25/B, Graz, 8010, Austria Vanja Subotić Email: vanja.subotic@tugraz.at Affiliation: Institute of Thermal Engineering, Graz University of Technology, Inffeldgasse 25/B, Graz, 8010, Austria Pavle Boškoski Email: pavle.boskoski@ijs.si Affiliation: Jožef Stefan Institute, Jamova cesta 39, Ljubljana, 1000, Slovenia Affiliation: Faculty of Information Studies in Novo mesto, Ljubljanska cesta 30, Novo mesto, 8000, Slovenia Abstract Estimating the distribution of relaxation times (DRT) from electrochemical impedance spectroscopy (EIS) is an ill-posed inverse problem that is highly sensitive to regularisation choices. We propose a physics-informed convolutional autoencoder that estimates DRTs directly from EIS data without spectrum-specific tuning. A discretised relation between impedance and the DRT is embedded in the training process, constraining the network to produce impedance-consistent distributions. The model resolves overlapping relaxation processes in synthetic two-ZARC spectra and accurately reconstructs measurements from three independent solid oxide fuel and electrolysis cell datasets, with range-normalised errors below 1.1%. Decoder-probe analysis shows that the learned latent representation is organised according to relaxation timescale. Distances in this latent space capture operating changes, hydrogen-shortage events, and long-term degradation. The same lightweight architecture is applied across all datasets without modification, providing consistent DRT estimation and an interpretable basis for condition monitoring. Keywords: Electrochemical impedance spectroscopy , Distribution of relaxation times , Solid oxide cells , Physics-informed neural networks , Condition monitoring , Latent representation 1 Introduction In electrochemistry, characterization of a device is often performed via electrochemical impedance spectroscopy (EIS), which provides the impedance or frequency response of the system. While EIS is a powerful diagnostic method, interpreting raw complex-valued impedance across many frequencies remains challenging. To deconvolve the raw EIS spectrum, a widely known method is to fit an equivalent circuit model (ECM) to the data. ECMs require quite specific prior assumptions for the identification, which is a particularly acute problem for fractional-order systems where no systematic structural identification strategy exists. Moreover, as discussed by 29, impedance based analysis is subject to uncertainty arising from measurement noise, operating variability, and modelling assumptions, which can further obscure parameter identifiability and interpretation. To overcome these limitations, distribution of relaxation times (DRT) has emerged as a widely adopted method for analysing EIS data 1. DRT decomposes the impedance spectrum into contributions associated with different relaxation time constants, offering a physically interpretable representation that replaces discrete circuit selection with a continuous distribution over an implicit Voigt structure 6. The relation between DRT and impedance Z(jω)Z(jω) is given by 23 Z(jω)=∫−∞g(logτ)1+jωτdlogτ,Z(jω)= _-∞^∞ g( τ)1+jωτ\,d τ, (1) where Z(jω)Z(jω) denotes the impedance, ω the angular frequency, logτ τ the logarithmic relaxation time, and g(logτ)g( τ) the DRT function. The motivation for using DRT is twofold. First, the kernel 1/(1+jωτ)1/(1+jωτ) represents the impedance of a simple (parallel) RC element, which relaxes exponentially with time constant τ under an impulse perturbation. Thus, g(logτ)g( τ) can be interpreted as a superposition of many such elementary processes, revealing the dominant relaxation modes of the system 6. This is analogous to representing complex dielectric relaxation behaviour using a distribution of Debye relaxation times 6. Second, unlike complex impedance, DRT maps each frequency to a single scalar value, simplifying visualisation and interpretation. Note that (1) omits explicit resistive and inductive terms to focus solely on the relaxation behaviour. However, extracting the DRT from measured impedance data is an ill-posed inverse problem. Small variations or noise in the impedance measurements can produce significantly different DRTs that nonetheless yield similar Nyquist curves. To obtain physically meaningful results, regularisation is typically introduced 7. The choice of regularisation strongly influences the resulting DRT, making the analysis sensitive to the selected method. Regularisation may be chosen based on self-consistency criteria 25 or cross-validation approaches 16. To avoid hand-tuned regularisation, data-driven approaches offer an attractive alternative 13, 12. Physics-based neural networks provide a way of employing the powerful approximation mechanisms of neural networks while at the same time constraining the model output to satisfy known physical relationships 24, 4. This concept is directly transferable to the estimation of the DRT. In this setting, the unknown function g(logτ)g( τ) is represented by a neural network while being constrained such that its integral satisfies the impedance relation in (1). A notable example is the approach proposed by 12, which estimates DRTs directly from individual EIS spectra using a Deep Prior neural network. The method solves the integral in (1) through a Fredholm approximation and proves to be effective even in the presence of noise and overlapping relaxation processes. However, because a separate network is trained for each spectrum, the implicit regularisation, induced by the neural network, varies from measurement to measurement, making it perilous to attribute differences in estimated DRTs to genuine electrochemical changes rather than fitting artefacts. Recent studies have demonstrated the effectiveness of learning shared latent representations from EIS data. Autoencoder architectures have been successfully applied to battery diagnostics, where latent variables extracted from impedance spectra capture information related to degradation and improve state of health (SOH) estimation 19, 14, 15. Several works have also shown that transforming impedance spectra into DRT representations via radial basis function (RBF) deconvolution prior to learning improves interpretability and predictive performance by separating overlapping electrochemical processes into distinct relaxation features 10, 11. However, these approaches either treat the latent space as a black-box feature extractor or rely on a separate deconvolution step. None embed the DRT integral relation directly into the learning objective. We propose a physics-constrained convolutional neural network (CNN) framework that estimates DRTs from multiple EIS spectra simultaneously. In contrast to approaches that fit a separate model to each spectrum, our method learns a shared latent representation across the entire dataset, enforcing consistent regularisation and enabling direct comparison of the resulting DRTs. The learned latent space provides a compact, interpretable basis for analysing temporal evolution and condition-dependent behaviour, while decoder probe analysis reveals that distinct latent channels correspond to specific relaxation time regions. We validate the framework on three independent solid oxide cell datasets, comprising steam electrolysis (603 spectra), co-electrolysis (150 spectra), and reversible operation (1248 spectra), demonstrating that predicted DRTs reproduce degradation trends established by classical analysis. 2 Convolutional formulation of the inverse problem Solving for g(logτ)g( τ) can be recast as a convolution, given in general by (h∗f)(t)=∫−∞h(τ)f(t−τ)τ.(h*f)(t)= _-∞^∞h(τ)\,f(t-τ)\,dτ. (2) The relation between DRT and impedance can be thus written as a convolution Z^(logτ^) Z( τ) =∫−∞g(logτ)dlogτ1+jexp[−(logτ^−logτ)] = _-∞^∞ g( τ)\,d τ1+j \! [-( τ- τ) ] =(g∗f)(logτ^). =(g*f)( τ). (3) where Z^(logτ^)=Z(jω) Z( τ)=Z(jω), τ^=1/ω τ=1/ω and f(x)=1/(1+jexp(−x))f(x)=1/ (1+j (-x ) ). Convolution theorem states that Fourier transform of a convolution of two functions is the product of their Fourier transform 21, i.e. ℱZ^=ℱgℱfF\ Z\=F\g\F\f\. Rearranging the equations we get ℱg=ℱZ^1ℱfF\g\=F\ Z\ 1F\f\, which we can again write out as a convolution g(logτ)=(Z^∗h)(logτ), g( τ)= ( Z*h )( τ), (4) where h=ℱ−11ℱfh=F^-1 \ 1F\f\ \ can be seen as a filter that filters out DRT from impedance. Since filters can be decomposed to what is know as cascade form [27, p. 61] (Z^∗h)(logτ) ( Z*h)( τ) ≃(((Z^∗h1)∗h2)∗h3∗⋯)(logτ), ((( Z*h_1)*h_2)*h_3·s )( τ), (5) where filters h1,h2,h3,…h_1,h_2,h_3,… can perform several different tasks such as denoising 28, compression-decompression 27 and in the end return a regularized solution [9, Chapter 6]. The cascade representation in (5) suggests that the inverse DRT problem can be approximated by a sequence of filtering operations. Since convolutional layers perform precisely such operations, CNNs are a natural choice for this task. This motivates the encoder-decoder architecture described in next section, where convolutional layers learn the cascade of filters directly from data. 3 Physics-informed convolutional neural network framework The proposed framework employs a CNN to estimate DRTs from EIS data. CNNs exploit local correlations through trainable kernels that slide along the frequency axis, reducing the number of parameters while capturing frequency-dependent features within a defined receptive field 22. Each kernel operates within a defined receptive field (RF) and captures relevant spectral patterns while substantially reducing the number of trainable parameters compared to fully connected architectures. The feature extraction process is controlled by the kernel size, stride, and padding, which determine the output dimensionality and receptive field 22. In this work, the real and imaginary components of the impedance are treated as two input channels (Figure 1). Each kernel spans both channels and slides along the frequency axis, allowing the network to jointly learn local patterns across both components of the complex impedance spectrum. Figure 1: CNN first layer schematic representation for kernel size 2, padding 1 and stride 1. Within the proposed framework, the CNN does not predict impedance directly. Instead, it estimates the DRT, g(logτ)g( τ), which is subsequently inserted into a discretised form of (1). To compute the corresponding impedance, the continuous integral is numerically approximated over a grid of logarithmic relaxation times using the trapezoidal rule, ∫abf(x)x≈∑k=1Nf(xk−1)+f(xk)2Δxk, _a^bf(x)\,dx≈ _k=1^N f(x_k-1)+f(x_k)2 x_k, (6) where f(x)f(x) corresponds to the integrand of (1) and xkx_k denotes sampled values of logτ τ. The complete training procedure is illustrated in Figure 2. The measured real and imaginary components of the impedance spectrum are provided as input to the CNN, which predicts the corresponding DRT. In addition to the DRT, a separate neural network predicts the series resistance RsR_s, while the inductive contribution L is represented by a trainable parameter. The predicted distribution is then inserted into the discretised electrochemical model and transformed back into the impedance domain using the trapezoidal approximation. The reconstructed impedance is compared with the measured spectrum to compute the training loss, ℒ=MSE(Z^′,Z′)+MSE(Z^′,Z′)+λLL2,L=MSE\! ( Z ,Z )+MSE\! ( Z ,Z )+ _LL^2, (7) where primes denote the real components, double primes denote the imaginary components, and λL=106 _L=10^6 penalises large inductive contributions. measuredEIS (Z′,Z′)(Z ,Z )CNN outputs g(logτ)g( τ)numerical solution ofZ(ω)=Rs+∫g(logτ)1−jωτ1+(ωτ)2dlogτ+jωLZ(ω)=R_s+\! g( τ) 1-jωτ1+(ωτ)^2d τ+jω Lreconstructedimpedance Z Zloss:Z Z vs. ZZRsR_s networkpredicts RsR_s from Z′Z trainableparameter L latent space analysis (post training) 2 channelsg(logτ)g( τ)RsR_sjωLjω L Figure 2: Schematic representation of the proposed physics-informed CNN framework. Measured impedance spectra are provided as two input channels and mapped by the CNN to the corresponding DRT, g(logτ)g( τ). In parallel, an auxiliary branch estimates the series resistance RsR_s from the real part of the impedance, while the inductive contribution L is represented by a trainable parameter. The predicted DRT, RsR_s, and L are inserted into the physics-informed reconstruction layer to obtain the reconstructed impedance Z Z. The model is trained by comparing Z Z with the measured impedance Z, enabling physics-informed learning through impedance reconstruction error. After training, the learned latent representation is used for interpretability and condition monitoring analyses. Because the impedance reconstruction layer consists entirely of differentiable operations, gradients can be propagated through the numerical integration procedure during backpropagation. Consequently, the parameters of the CNN, the resistance model RsR_s, and the inductive parameter L are jointly optimised using the Adam optimiser. As a result, the network is trained to produce DRTs that are consistent with the underlying electrochemical relationship rather than merely fitting the impedance spectrum. After training, the encoder additionally provides a compact latent representation of the impedance data. This latent space can be analysed independently of the reconstruction task and is used in the subsequent sections to investigate latent space organisation, condition monitoring, and long term system evolution. By embedding the governing electrochemical relation directly into the optimisation procedure, the network is constrained to produce physically consistent DRTs while retaining the flexibility of a data driven model. 3.1 Mitigating checkerboard artifacts Transposed convolution layers can produce checkerboard artifacts when kernel size, stride, and padding lead to uneven overlap during upsampling [20]. In DRT estimation, such artifacts manifest as spurious oscillations that may be misinterpreted as physical relaxation processes (Figure 3a). To eliminate this effect, each transposed convolution was replaced by an upsampling layer followed by a standard convolution, producing smooth reconstructions without artificial structures (Figure 3b). The checkerboard artifacts are effectively removed while preserving the overall shape and locations of the reconstructed relaxation processes. Figure 3: Comparison of decoder architectures. (a) Reconstructions obtained using transposed convolutions exhibit characteristic checkerboard artifacts. (b) Replacing transposed convolutions with an upsampling layer followed by a standard convolution removes these artifacts and produces smoother DRT estimates. Beyond improving reconstruction quality, this modification promotes smoother and more physically plausible DRT estimates, which is important for the reliable interpretation of the resulting distributions. 3.2 Architecture selection In solving the inverse DRT problem, the primary challenge is not merely minimizing the reconstruction error but preventing the amplification of measurement noise. Models with excessive capacity tend to fit stochastic fluctuations in the impedance spectra, producing non-physical oscillations in the reconstructed DRTs. Conversely, models with insufficient capacity fail to resolve overlapping relaxation processes. To determine an appropriate balance between reconstruction fidelity and model complexity, a grid search across 216 network architectures was performed. The search space included variations in network depth, channel width, bottleneck size, pooling strategy, and receptive field. The relationship between model complexity and validation loss is shown in Figure 4. A distinct elbow is observed at approximately 10410^4 trainable parameters. Below this threshold, the model underfits the impedance spectra and exhibits elevated reconstruction error. Beyond approximately 1.5×1041.5× 10^4 parameters, improvements become marginal, indicating diminishing returns and an increased risk of overfitting. Figure 4: Validation loss as a function of model complexity. A more detailed analysis of the architectural search is provided in Section S1 of the Supplementary Material. Based on the complete study, a receptive field of 57 and a bottleneck size of 160 were selected as the best compromise between reconstruction accuracy, stability, and resistance to overfitting. The resulting architecture contains approximately 12,000 trainable parameters and was used throughout the remainder of this work. Importantly, the selected architecture was not optimized for a specific experiment. Instead, the same network configuration was subsequently applied to all investigated datasets without any dataset-specific redesign or hyperparameter tuning. This provides a stringent test of the framework’s ability to generalize across different operating conditions, degradation mechanisms, and experimental campaigns. 3.3 Experimental datasets To evaluate the generalisation capability of the proposed framework, the selected architecture was applied to three independent experimental datasets covering different operating regimes, degradation mechanisms, and experimental objectives. Solid-oxide fuel cell (SOFC) stack monitoring campaign. The first dataset consists of 603 EIS measurements acquired during a 3600-hour operation of a six-cell SOFC stack, including several fuel starvation events, emergency shutdowns due to H2 shortage, and an unplanned power loss 18, 3, 8. This dataset was used to evaluate impedance reconstruction accuracy, latent space condition monitoring, and event detection. Operating condition study. The second dataset corresponds to the experimental study presented in 2, where an industrial-scale solid-oxide electrolyser cell (SOEC) (100 cm2) was operated under three representative conditions: high hydrogen content at moderate current density (Condition 1, 700 mA cm-2, 20% H2, 60 h), low hydrogen content at the same current density (Condition 2, 700 mA cm-2, 10% H2, 85 h), and elevated current density with unstable steam supply (Condition 3, 900 mA cm-2, 120 h). Changes were originally characterised using frequency-resolved Kullback–Leibler (KL) divergence analysis, and here the dataset is used to examine whether the learned latent representation captures the same operating-state transitions. Multi-regime degradation campaign. The third dataset consists of 1248 EIS spectra collected during a 2650-hour experimental campaign on a commercial 4×44×4 cm2 electrolyte-supported cell 17. The campaign comprised six sequential degradation phases: steam electrolysis at 300 mA cm-2 (P1, 1000 h), co-electrolysis at 300 mA cm-2 (P2, 400 h), and 500 mA cm-2 (P3, 200 h), a return to steam electrolysis (P4, 300 h), reversible EC/FC operation at 300/−-150 mA cm-2 (P5, 264 h), and a final steam electrolysis phase (P6, 300 h). This dataset was used to investigate long term DRT evolution and to assess whether latent space trajectories capture transitions between different operating regimes and degradation processes. 4 Results and discussion The results are presented in four stages. First, the framework is evaluated on experimental data to verify impedance reconstruction accuracy, physics consistency, and comparison with existing DRT estimation methods. Second, the proposed framework is validated on synthetic spectra with analytically known distributions to assess its ability to resolve overlapping relaxation processes. Third, the learned latent representation is analysed to investigate its physical organisation and suitability for condition monitoring. Finally, the framework is applied to independent experimental datasets to assess the generalisability of the complete analysis pipeline. 4.1 Core validation and physics consistency The primary objective of the proposed physics-informed framework is to estimate DRTs that remain consistent with the governing electrochemical relation in (1). Since the predicted DRTs are transformed back into the impedance domain through the embedded physical model, accurate reconstruction of the measured impedance spectra provides direct evidence that the network satisfies the underlying electrochemical constraints. Representative reconstruction results obtained for the SOFC stack monitoring dataset introduced in Section 3.3 are shown in Figure 5. The reconstructed Nyquist spectra closely follow the measured impedance response, while the residual errors remain small across the investigated frequency range. Figure 5: Representative reconstruction results for the SOFC stack monitoring dataset. (a) Measured and reconstructed Nyquist spectra for selected measurements. (b) Frequency-dependent reconstruction error, |Z−Z^||Z- Z|. The low residuals confirm accurate impedance reconstruction from the predicted DRTs. The corresponding Nyquist and residual plots for the operating condition study and the multi-regime degradation campaign are provided in Section S2 of the Supplementary Material. Reconstruction accuracy across all three datasets was quantified using the root mean squared error (RMSE) between measured and reconstructed impedance spectra, complemented by a range-normalised RMSE to enable comparison across datasets with different impedance magnitudes. Dataset RMSE (Re) RMSE (Im) RMSE (total) R-RMSE [mΩ ] [mΩ ] [mΩ ] [%] SOFC stack monitoring 0.536 0.191 0.402 0.291 Operating condition study 1.461 1.741 1.607 1.062 Multi-regime degradation 0.504 0.453 0.479 0.451 Table 1: Impedance reconstruction accuracy for the investigated datasets. RMSE values are reported together with the range-normalised RMSE (R-RMSE), expressed relative to the impedance range of each dataset. Despite substantial differences in operating conditions, degradation mechanisms, and experiment duration, the reconstruction errors remain consistently low across all datasets. In all cases, the range-normalised RMSE remains below approximately 1.1%1.1\%, demonstrating excellent agreement between the measured and reconstructed impedance spectra. These results provide empirical validation of the physics-informed nature of the proposed architecture. To further assess the quality of the estimated DRTs, the proposed approach was compared with conventional Tikhonov regularisation and a previously reported neural network-based method. Figure 6: Evolution of DRTs over time using (a) Tikhonov regularisation 26, (b) a feed-forward neural network-based approach 12, 8, and (c) the proposed CNN framework. Figure 6 compares DRT estimates obtained using Tikhonov regularisation 7, a feed-forward neural network approach 12, and the proposed CNN framework. All three methods identify the dominant relaxation process, visible as the pronounced band between approximately 1 and 10 Hz, but differ substantially in how they represent its temporal evolution. Tikhonov regularisation (Figure 6a) produces a smooth, nearly stationary distribution, suppressing much of the temporal variation observed throughout the experiment. The feed-forward neural network (Figure 6b) recovers considerably richer structure, including the gradual emergence of a secondary lower-frequency relaxation process, but at the cost of several abrupt transitions between neighbouring experiment indices, most notably early in the experiment and again around index 200. The proposed CNN framework (Figure 6c) captures the same temporal features as the feed-forward approach while evolving more continuously between successive measurements. Abrupt transitions still appear at a few experiment indices, but they are markedly less pronounced. This continuity follows directly from training a single shared model across the entire dataset. Rather than fitting each spectrum with an independently chosen regularisation strength, as Tikhonov regularisation and per-spectrum neural networks both do, the proposed framework applies the same implicit regularisation to every measurement, allowing neighbouring spectra to be compared directly rather than through the lens of separately tuned fits. The value of this consistency lies not in superior point-wise accuracy, since all three methods agree closely on the dominant relaxation process, but in ensuring that differences observed across the dataset reflect genuine electrochemical changes rather than fitting artefacts. This is confirmed quantitatively by the sub-1.1%1.1\% reconstruction error achieved throughout, showing that the network has learned the governing electrochemical relationship rather than a statistical shortcut. 4.2 Validation on synthetic overlapping relaxation processes The ability to distinguish closely spaced relaxation processes was further assessed using synthetic two-ZARC spectra with analytically known DRTs. A ZARC element, consisting of a resistor in parallel with a constant-phase element, is described by ZZARC(ω)=R1+(jω)ϕRQ.Z_ZARC(ω)= R1+(jω)^φRQ. (8) The total impedance is given by Z(ω)=Rs+ZZARC,1(ω)+ZZARC,2(ω)Z(ω)=R_s+Z_ZARC,1(ω)+Z_ZARC,2(ω), where R is the process resistance, Q is the constant-phase-element parameter, ϕφ is the dispersion parameter, and RsR_s is the series resistance. For a known impedance transfer function, the corresponding analytical DRT was obtained by analytic continuation using the Fuoss–Kirkwood formulation 5, G(u) G(u) =−1π[ℑZ(e−u−jπ/2)+ℑZ(e−u+jπ/2)], =- 1π [ \! \Z\! (e^-u-jπ/2 ) \+ \! \Z\! (e^-u+jπ/2 ) \ ], (9) where u=logτu= τ. This analytical solution provides an exact reference against which the recovered peak structure and relaxation-time positions can be compared. To evaluate the proposed framework, the architecture selected in Section 3.2 was trained on 1,000 synthetic two-ZARC spectra generated by varying the resistances, dispersion parameters ϕφ, and the logarithmic separation between the characteristic relaxation times. Peak separations spanning 0.2–2.5 decades were included to cover both strongly overlapping and well-separated relaxation processes. The complete parameter ranges and data generation procedure are provided in Appendix A. The validation spectra presented below were generated separately and were not included in the training dataset. Following the validation protocol of 12, two representative test cases were considered, with relaxation times separated by two decades and one decade, respectively. To evaluate the network under non-ideal conditions, independent Gaussian noise with a standard deviation of σ=0.5Ωσ=0.5~ was added to both the real and imaginary impedance components. This corresponds to the upper limit of the noise range used during training and therefore represents a more demanding validation than the noise-free spectra considered in the original study. The results are shown in Figure 7. For the two-decade separation, the two relaxation processes are already reflected by partially separated arcs in the Nyquist spectrum (Figure 7a). The corresponding network-estimated DRT accurately reproduces the analytical two-peak structure, with both relaxation processes clearly identified (Figure 7b). Although the higher-τ peak is shifted by 0.241 decades, the two processes remain clearly resolved. The one-decade case provides a considerably more challenging example. Here, the two ZARC contributions overlap strongly and appear as a single broad depressed arc in the Nyquist representation (Figure 7c), such that the presence of two underlying relaxation processes is not readily apparent from visual inspection of the impedance spectrum. Nevertheless, the network-estimated DRT reproduces the analytical two-peak structure with good agreement (Figure 7d), with peak-position errors of only 0.065 and 0.025 decades. These results demonstrate that the proposed framework accurately recovers the analytical DRT for overlapping two-ZARC spectra. In the one-decade case, the network reproduces the correct two-peak structure even though the corresponding Nyquist spectrum appears as a single broad depressed arc. This indicates that the learned impedance-to-DRT mapping exploits subtle information contained across the full impedance spectrum rather than relying solely on features that are visually apparent in the Nyquist representation. Figure 7: Validation on synthetic two-ZARC spectra with additive Gaussian noise (σ=0.5Ωσ=0.5~ ), using the test configurations considered by 12. (a)–(b) Two-decade separation with τ1=0.1 _1=0.1 s and τ2=10 _2=10 s. (c)–(d) One-decade separation with τ1=0.1 _1=0.1 s and τ2=1 _2=1 s. The left panels compare the analytical impedance spectra with the network reconstructions, while the right panels compare the analytical and network-estimated DRTs. 4.3 Interpretable latent space dynamics The question is whether the learned latent representation goes further by containing interpretable information about the electrochemical system beyond what is needed for reconstruction. In particular, we examine whether the latent space exhibits a meaningful internal structure and whether its evolution can be used to monitor changes in system operation. At time t, the encoder maps the input impedance spectrum to a latent tensor t∈ℝC×PZ_t ^C× P, where C=40C=40 denotes the number of output channels produced by the final convolutional layer and P=4P=4 the pooled positions retained after adaptive pooling. Each channel represents one learned component of the encoded impedance spectrum, while the pooled positions preserve coarse information about where that component originated along the frequency axis before pooling. Flattening tZ_t row-wise yields the latent vector t∈ℝCPz_t ^CP used throughout the remainder of this section. To investigate the physical meaning of these latent components, a decoder probe analysis was performed. For every channel–position pair (c,p)(c,p), an artificial latent tensor ~∈ℝC×P Z ^C× P was constructed with all entries set to zero except Z~c,p Z_c,p, which was assigned a fixed activation. Passing this tensor through the decoder reveals the relaxation-time response associated with activating only that individual latent component. Repeating this procedure for all 40×440× 4 channel–position pairs produces the activation maps shown in Figure 8, where each panel corresponds to one pooled position and each row to one latent channel. The channels are ordered according to the location of their peak response at pooled position 0, and this ordering is retained across all panels to facilitate comparison. Figure 8: Decoder probe analysis of the learned latent basis. For each pooled position, every latent channel was activated individually while all remaining latent entries were fixed to zero. Each row corresponds to one latent channel, ordered according to the location of its peak response at pooled position 0. This ordering is preserved across all panels. Colour indicates the amplitude of the reconstructed relaxation-time response produced by activating the corresponding channel–position pair. The decoder associates the four pooled positions with progressively different regions of the relaxation-time axis. For each channel, the peak of the reconstructed response was identified, and the median peak location across all channels was computed for each pooled position (the median was used because relaxation times span several decades). Pooled positions 0 and 1 predominantly activate long-timescale processes, with median peak responses at approximately τ=24τ=24 s and τ=83τ=83 s, respectively, whereas positions 2 and 3 shift towards shorter timescales, with median peak responses near τ=3.5τ=3.5 s and τ=0.3τ=0.3 s. Within each pooled position, different channels respond to different parts of the corresponding relaxation-time range, providing a finer subdivision of the represented electrochemical processes. The latent representation therefore exhibits a hierarchical organisation: the pooled position determines a coarse relaxation-time region, while the channel identity provides additional resolution within that region. Importantly, this organisation was not imposed during training. The network was optimised solely to reconstruct the measured impedance, without any supervision encouraging latent components to specialise by relaxation timescale. The emergence of this structure indicates that the learned latent space reflects the physical organisation of the underlying electrochemical processes rather than constituting an arbitrary compressed representation. This physically organised latent representation also provides a natural space for quantifying changes in system behaviour over time. Rather than comparing impedance spectra directly, system evolution can be monitored by measuring distances between the corresponding latent representations. The Euclidean distance between the latent vector at time t and a chosen reference vector refz_ref is defined as Dref(t)=‖t−ref‖2.D_ref(t)= _t-z_ref _2. (10) To demonstrate this, the SOFC stack monitoring dataset was analysed using D10(t)D_10(t), where the 10th measurement serves as the nominal reference state. Figure 9: Latent space condition monitoring for the SOFC stack monitoring dataset. D10(t)D_10(t) is used as an indicator of system evolution. Pronounced peaks correspond to major operational events during the campaign. The resulting distance profile exhibits pronounced peaks that coincide with known operational events during the experimental campaign, summarised in Table 2. Event number Experiment time [h] Description E1 0–40 Start up E2 270–342 Increased fuel utilisation (decreased fuel flow rate) E3 582–648 Increased fuel utilisation (increased current) E4 654–672 Emergency shutdown due to H2 shortage E5 1242 Power loss due to thunderstorms E6 1434 Emergency shutdown due to H2 shortage E7 1884–1890 Moving into the new building E8 2202 Emergency shutdown due to H2 shortage Table 2: Operational events corresponding to the peaks observed in Figure 9 [18, same as Table 1]. Increases in D10(t)D_10(t), seen as peaks in Figure 9, coincide with fuel starvation events, shutdowns, and other major perturbations, demonstrating that deviations from the nominal operating state are reflected directly and immediately in the latent representation. These observations are consistent with the conclusions of 8, where significant events were detected through comparisons of reconstructed DRTs using unbalanced optimal transport (UOT). While UOT provides a principled metric for comparing DRTs, the latent space distance (10) offers a computationally efficient alternative obtained directly from the encoder output, avoiding ambiguous DRT curve comparisons. The decoder probe analysis and the condition monitoring results confirm that the learned latent representation is both physically interpretable and operationally informative. The network not only reconstructs impedance spectra accurately, but also spontaneously organises electrochemical information into a structured latent space, a physically interpretable representation that was never explicitly supervised and yet reliably identifies degradation and operational events. 4.4 Pipeline generalisation The previous sections established that the proposed framework produces physically consistent DRT estimates and learns an interpretable latent representation from a single experimental campaign. The critical remaining question is whether the same pipeline can extract equally meaningful information from fundamentally different experiments, without any modification to the architecture, hyperparameters, or analysis workflow. To answer this, the framework was applied to two independent datasets that differ from the SOFC stack monitoring campaign in cell type, operating regime, degradation mechanism, and experimental objective. The results therefore constitute a direct test of architectural generalisation, as the same 12,000-parameter network, trained independently on each dataset, required no redesign, retuning, or workflow modification. Operating condition study 2 Figure 10: Latent space evolution for the three operating conditions investigated in 2. Heatmaps show the smoothed absolute difference between each latent representation and the initial spectrum, |t−0| |z_t-z_0 |, where brighter regions indicate larger cumulative deviation from the initial operating state. Dashed lines denote characteristic time points reported in the original analysis. This dataset (Section 3.3) comprised three operating conditions investigated over 60–125 h each on an industrial-scale SOEC. In the original study, changes in system behaviour were identified using KL divergence computed frequency by frequency relative to an initial reference spectrum, identifying frequency bands where the impedance changed statistically significantly 2. Here, the same intuition is applied in latent space. Recall from Section 4.3 that each spectrum’s latent tensor is flattened into a vector t∈ℝCPz_t ^CP. The scalar distance Dref(t)D_ref(t) from (10) summarises the overall change between tz_t and a reference refz_ref as a single number, but does not indicate which part of the latent representation changed. To recover this information, we instead compute the element-wise absolute difference |t−0| |z_t-z_0 |: rather than combining all CPCP entries into one value, this keeps a separate difference for each entry, so that changes can be localised to specific channels and pooled positions, in the same way that the original study localises changes to specific frequency bands. This correspondence is temporal rather than spectral: the original study attributes each transition to a specific frequency band, whereas here we only identify the same transition times through changes in latent activation, without establishing which latent channels or positions correspond to which frequency band. Unlike the original approach, this requires neither frequency-band selection nor numerical integration over the impedance spectrum, since the diagnostic information is extracted directly from the encoder output as a byproduct, with no additional analysis required. The resulting activation maps are shown in Figure 10. For Condition 1, diffusion-related losses emerging after approximately 20 h (0.3–2 Hz) and air-electrode changes after approximately 40 h (10–40 Hz), originally identified through KL divergence 2, coincide in time with changes in latent activation visible in Figure 10a. The same holds for Condition 2, where the gradual convergence toward a new equilibrium, marked by spectral changes shifting from 200–800 Hz near 40 h to 800 Hz–4 kHz near 70 h 2, coincides with two distinct transitions in Figure 10b. Condition 3 differs in character: the persistent pulsating instability caused by inconsistent steam supply at 900 mA cm-2 2 does not produce discrete transition events. Instead, Figure 10c shows sustained, recurring latent activity throughout the experiment, mirroring the continuous operational instability rather than isolated degradation onsets. Multi-regime degradation campaign 17 Figure 11: Difference between the network-estimated DRTs and the first estimated DRT within each operating regime (Δgt(logτ) g_t( τ), (11)). Horizontal axis shows elapsed experiment time within each phase. This dataset (Section 3.3) comprised 1248 EIS spectra collected on a commercial 4×44×4 cm2 electrolyte-supported cell, with degradation characterised phase by phase through bi-hourly DRT diagnostics 17. Table 3 summarises the six operating phases and their correspondence to the panels in Figure 11 and Figure 12. Phase Panels Experiment time [h] Operating regime P1 a) 22–1022 Steam electrolysis, 300 mA cm-2 P2 b) 1055–1455 Co-electrolysis, 300 mA cm-2 P3 c) 1473–1673 Co-electrolysis, 500 mA cm-2 P4 d) 1701–2001 Steam electrolysis, 300 mA cm-2 P5 e)–f) 2048–2312 Reversible EC/FC operation P6 g) 2335–2635 Final steam electrolysis, 300 mA cm-2 Table 3: Operating phases of the degradation campaign investigated in 17 and their correspondence to the panels in Figure 11 and Figure 12. To enable direct comparison with the phase by phase degradation trends reported in 17, temporal evolution within each phase is quantified using the same metric employed in the original study, Δgt(logτ)=gt(logτ)−g1(logτ), g_t( τ)=g_t( τ)-g_1( τ), (11) where g1(logτ)g_1( τ) is the first network-estimated DRT of that phase, serving as the baseline specific to that phase. While the latent representation introduced in Section 4.3 offers a better-conditioned basis for tracking system evolution, adopting the original study’s own metric here allows the estimated DRTs to be benchmarked directly against their published results. This comparison is further strengthened by the constant regularisation of the proposed framework, since all 1248 spectra share the same network weights and Δgt g_t therefore reflects genuine electrochemical change and nothing else. A method applied per spectrum would impose a different regularisation on each measurement, making it impossible to determine whether differences in the estimated DRTs reflect genuine electrochemical changes or fitting artefacts. The evolution of Δgt(logτ) g_t( τ) tells a coherent degradation story. Through the first three phases, changes accumulate at the low-frequency end of the distribution: co-electrolysis (P2) adds a new contribution there, and the higher current density of P3 shifts it without amplifying it [17, Figure 5a–c], both captured in Figure 11a–c. Returning to steam electrolysis (P4) partially undoes this, with the low-frequency contribution retreating and confirming that not all co-electrolysis-induced losses are permanent [17, Figure 5d] (Figure 11d). The most pronounced behaviour occurs during reversible EC/FC operation (P5): the alternating modes leave a strongly periodic imprint on the DRT, with pronounced high-frequency modulation tied to FC segments [17, Figure 5e] (Figure 11e–f). P6 closes the campaign quietly, with only moderate evolution reflecting the slow pace of ohmic-dominated ageing [17, Figure 5f]. Importantly, the same degradation trajectory is independently recovered by the latent space, without any DRT computation. Just as in the SOFC stack monitoring campaign, D10(t)D_10(t) (10) reduces the entire 2650-hour campaign to a single scalar signal shown in Figure 12, where the slope steepens at the onset of co-electrolysis, relaxes upon recovery, oscillates through the EC/FC phase, and settles into a slow drift in P6, tracing every regime transition identified in the DRT analysis. Figure 12: Evolution of D10(t)D_10(t) during the 2650-hour campaign. Coloured regions denote the operating regimes in 17. Summary Table 4 places the key findings of the proposed framework alongside the original results to make the correspondence explicit and verifiable. Across all three experimental campaigns, spanning different cell types, operating regimes, degradation mechanisms, and campaign durations from 60 h to 3600 h, the proposed framework reproduces the key findings of the original studies without any dataset-specific tuning. Notably, the framework required neither manual regularisation parameter selection, frequency-band segmentation, equivalent circuit model assumptions, nor dataset-specific architectural modifications. Original study This work SOFC stack monitoring campaign 18, 3, 8 Fuel-starvation and shutdown events via UOT 8 Peaks in D10(t)D_10(t) (Figure 9) No prior latent space analysis Distinct channels per relaxation region (Figure 8) Operating condition study 2 Diffusion onset ∼ 20 h, 0.3–2 Hz (Cond. 1) Latent activation change at ∼ 20 h (Figure 10a) Air-electrode changes ∼ 40 h, 10–40 Hz (Cond. 1) Latent activation change at ∼ 40 h (Figure 10a) New equilibrium at ∼ 40 h and ∼ 70 h (Cond. 2) Two transitions in latent activation (Figure 10b) Pulsating instability throughout (Cond. 3) Sustained latent activity throughout (Figure 10c) Multi-regime degradation campaign 17 Low-freq. DRT grows, case (P2) [17, Figure 5b] Δgt g_t increase at low frequencies (Figure 11b) Low-freq. DRT retreats, case (P4) [17, Figure 5d] Δgt g_t decrease at low frequencies (Figure 11d) Periodic high-freq. modulation, (P5) [17, Figure 5e] Periodic DRT pattern (Figure 11e–f) Moderate evolution, ohmic ageing (P6) [17, Figure 5f] Slow latent drift (Figure 12) Table 4: Correspondence between findings reported in the original studies and signatures identified by the proposed framework. All results in the This work column were obtained with the same 12,000-parameter architecture and zero dataset-specific tuning. 5 Conclusion This work addressed the problem of estimating the DRT from EIS data, an ill-posed inverse task that traditionally relies on careful regularisation or strong modelling assumptions. We proposed a physics-informed CNN that directly predicts the DRT and reconstructs the corresponding impedance through a discretised formulation of the governing integral relation. By embedding the physical forward model into the training process, the network is constrained to produce impedance-consistent DRTs without imposing a fixed number or shape of relaxation peaks, allowing it, in principle, to represent overlapping relaxation processes. A systematic architectural study showed that physically meaningful reconstructions depend strongly on the imposed inductive bias. Separating upsampling from convolution in the decoder eliminated checkerboard artifacts, while appropriate choices of bottleneck size and receptive field prevented both over-smoothing and noise amplification. The resulting 12,000-parameter architecture achieves a favourable balance between reconstruction accuracy and generalisation, with range-normalised reconstruction errors below 1.1% across all three datasets. What distinguishes this framework from prior approaches is not reconstruction accuracy alone, but what emerges from the shared latent representation. Because all spectra in a dataset share the same network weights, differences in estimated DRTs are constrained to reflect genuine electrochemical changes. This consistency, absent in per-spectrum methods, is what makes the latent space meaningful, as decoder probe analysis revealed that individual channels correspond to distinct relaxation-time regions. Such a structure emerged solely through optimisation of the physics-informed training objective, directly from the raw EIS data. The resulting latent space enables condition monitoring as a direct byproduct. Euclidean distances in latent space reliably identified every major operational event in a 3600-hour SOFC campaign, matching conclusions previously obtained through explicit DRT comparison using UOT 8. Most significantly, the same architecture was applied without modification to three independent experimental datasets spanning fundamentally different cell types, operating regimes, degradation mechanisms, and campaign durations from 60 h to 2650 h. In every case, the framework reproduced the key degradation signatures, operational transitions, and DRT peak evolutions previously established by dedicated per-dataset analyses, including manual regularisation tuning, frequency-band segmentation, expert interpretation, and the frequency-resolved KL divergence analysis of the operating condition study. None of those steps were required here. Physics-informed learning can therefore replace the manual preprocessing pipeline that has long been the bottleneck in large-scale EIS diagnostics. A single, unified, and fully reproducible model is sufficient to extract physically meaningful DRTs, monitor system health, and generalise across the diverse operating conditions encountered in electrochemical energy conversion research. CRediT authorship contribution statement Žan Gorenc: Conceptualization, Methodology, Software, Validation, Formal analysis, Investigation, Visualization, Writing – original draft. Žiga Gradišar: Methodology, Writing – review and editing. Felix Mütter: Investigation, Resources. Vanja Subotić: Supervision, Investigation, Resources, Funding acquisition. Pavle Boškoski: Conceptualization, Methodology, Software, Resources, Supervision, Funding acquisition, Writing – review and editing. Data and code availability The source code, pretrained models, reproducibility notebooks, and instructions for downloading the datasets required to reproduce the results presented in this work are publicly available at https://repo.ijs.si/e2pub/cnn_drt.git. Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Acknowledgements The authors gratefully acknowledge support for the project ”Probabilistic and explainable data-driven modelling of solid-oxide systems”, jointly financed by the Slovenian Research and Innovation Agency (ARIS), project number J2-4452, and the Austrian Science Fund (FWF), project number I 6251-N. The authors alsoacknowledge financial support from the Slovenian Research and Innovation Agency through research programme P2-0001. Appendix A Synthetic two-ZARC dataset The synthetic dataset used for the overlapping-process validation in Section 4.2 comprised 1,000 two-ZARC impedance spectra. The ZARC parameters were sampled over the ranges summarised in Table 1. The resistances and dispersion parameters were sampled independently from uniform distributions, while the first characteristic relaxation time was sampled uniformly in logarithmic space. For the ZARC formulation introduced in Section 4.2, the characteristic relaxation time is related to the constant-phase-element parameter by τk=(RkQk)1/ϕk. _k= (R_kQ_k )^1/ _k. (12) The second relaxation time was determined from a randomly sampled logarithmic separation according to τ2=τ110Δlog10τ. _2= _110 _10τ. (13) Parameter Range Sampling RsR_s [Ω ] 5–15 Uniform RkR_k [Ω ], k=1,2k=1,2 30–70 Uniform ϕk _k, k=1,2k=1,2 0.70–1.00 Uniform log10(τ1/s) _10( _1/s) −2-2–1 Uniform Δlog10τ _10τ 0.2–2.5 Uniform σ [Ω ] 0–0.5 Uniform Table 1: Parameter ranges and sampling distributions used to generate the synthetic two-ZARC training dataset. Each spectrum was evaluated at 200 logarithmically spaced frequencies between 10−410^-4 and 10410^4 Hz. Independent zero-mean Gaussian noise was added to the real and imaginary impedance components, Znoisy=Z+ϵre+jϵim,ϵre,ϵim∼(0,σ2),Z_noisy=Z+ _re+j _im, _re, _im (0,σ^2), (14) with σ sampled independently for each spectrum from the range reported in Table 1. The relaxation-time grid used for physics-based impedance reconstruction comprised 200 logarithmically spaced points between 10−510^-5 and 10410^4 s. The dataset was divided into training and validation subsets using an 80:20 split with a fixed random seed. The two test configurations evaluated in Section 4.2 were generated separately and excluded from the training dataset. Declaration of generative AI and AI-assisted technologies in the manuscript preparation process During the preparation of this work, the authors used Claude to assist with selected coding tasks, code documentation, and streamlining the development workflow. Grammarly was subsequently used to review the manuscript for grammar, spelling, and writing style. All outputs generated with the assistance of these tools were critically reviewed, verified, and edited by the authors. The authors take full responsibility for the accuracy, integrity, and final content of the manuscript and the accompanying code. References Barsukov et al. (2012) Y. Barsukov J. R. Macdonald et al. Electrochemical impedance spectroscopy. Characterization of materials 2, p. 898–913. Cited by: §1. Boskoski et al. (2024) P. Boskoski, B. Königshofer, G. Nusev, A. Ostrež, and V. Subotić Extracting quantitative insights from electrochemical impedance spectra using statistical methods. In 16th Annual Conference of the Prognostics and Health Management Society, PHM 2024, Cited by: §3.3, Figure 10, §4.4, §4.4, §4.4, Table 4. Boškoski et al. (2024) P. Boškoski, L. Žnidarič, Ž. Gradišar, and V. Subotić Probabilistic deconvolution for electrochemical impedance through variational bayesian inference. Journal of Power Sources 622, p. 235359. External Links: ISSN 0378-7753, Document, Link Cited by: §3.3, Table 4. Cuomo et al. (2022) S. Cuomo, V. S. di Cola, F. Giampaolo, G. Rozza, M. Raissi, and F. Piccialli Scientific machine learning through physics-informed neural networks: where we are and what’s next. arXiv. External Links: Document, Link Cited by: §1. Fuoss and Kirkwood (1941a) R. M. Fuoss and J. G. Kirkwood Electrical properties of solids. viii. dipole moments in polyvinyl chloride-diphenyl systems. Journal of the American Chemical Society 63 (2), p. 385–394. Cited by: §4.2. Fuoss and Kirkwood (1941b) R. M. Fuoss and J. G. Kirkwood Electrical properties of solids. viii. dipole moments in polyvinyl chloride-diphenyl systems*. Journal of the American Chemical Society 63 (2), p. 385–394. External Links: ISSN 1520-5126, Link, Document Cited by: §1, §1. Gavrilyuk et al. (2017) A. L. Gavrilyuk, D. A. Osinkin, and D. I. Bronin The use of tikhonov regularization method for calculating the distribution function of relaxation times in impedance spectroscopy. Russian Journal of Electrochemistry 53 (6), p. 575–588. External Links: ISSN 1608-3342, Document, Link Cited by: §1, §4.1. Gradišar et al. (2026) Ž. Gradišar, Ž. Gorenc, V. Subotić, and P. Boškoski Towards robust comparison of distributions of relaxation times. Electrochimica Acta, p. 148339. External Links: ISSN 0013-4686, Document, Link Cited by: §3.3, Figure 6, §4.3, Table 4, Table 4, §5. Hansen et al. (2006) P. C. Hansen, J. G. Nagy, and D. P. O’Leary Deblurring images: matrices, spectra, and filtering. SIAM. External Links: Document Cited by: §2. Kim et al. (2025) K. Kim, K. Park, K. Yoon, H. Moon, H. Pyeon, J. Kim, and S. Han A light-weight electrochemical impedance spectroscopy-based soh estimation method for lithium-ion batteries using the distribution of relaxation times with grad-cam analysis. Journal of Power Sources 631, p. 236169. External Links: ISSN 0378-7753, Document Cited by: §1. Li et al. (2025) Y. Li, Z. He, M. Ye, Q. Wang, G. Lian, Y. Sun, and M. Wei A semi-supervised learning strategy for lithium-ion battery capacity estimation with limited impedance data. Energy 319, p. 135129. External Links: ISSN 0360-5442, Link, Document Cited by: §1. Liu and Ciucci (2020a) J. Liu and F. Ciucci The deep-prior distribution of relaxation times. Journal of The Electrochemical Society 167 (2), p. 026506. External Links: ISSN 1945-7111, Document Cited by: §1, §1, Figure 6, Figure 7, §4.1, §4.2. Liu and Ciucci (2020b) J. Liu and F. Ciucci The gaussian process distribution of relaxation times: a machine learning tool for the analysis and prediction of electrochemical impedance spectroscopy data. Electrochimica Acta 331, p. 135316. Cited by: §1. Liu et al. (2024) Y. Liu, Q. Li, and K. Wang Revealing the degradation patterns of lithium-ion batteries from impedance spectroscopy using variational auto-encoders. Energy Storage Materials 69, p. 103394. External Links: ISSN 2405-8297, Link, Document Cited by: §1. Lou et al. (2025) C. Lou, S. Wang, Z. Li, and K. Wang Extract features from lithium-ion battery electrochemical impedance spectra and estimate state of health based on improved convolutional autoencoder-temporal convolutional network. Ionics 31 (5), p. 4261–4279. External Links: ISSN 1862-0760, Link, Document Cited by: §1. Maradesa et al. (2023) A. Maradesa, B. Py, T. H. Wan, M. B. Effat, and F. Ciucci Selecting the regularization parameter in the distribution of relaxation times. Journal of The Electrochemical Society 170 (3), p. 030502. Cited by: §1. Mütter et al. (2026) F. Mütter, P. Boškoski, H. Schröttner, C. Hochenauer, and V. Subotić Bi-hourly eis-drt diagnostics over 2650 h benchmark degradation of a commercial soec under steam electrolysis, co-electrolysis and reversible operation. Journal of Power Sources 672, p. 239640. External Links: ISSN 0378-7753, Document, Link Cited by: §3.3, Figure 12, §4.4, §4.4, §4.4, §4.4, Table 3, Table 4, Table 4, Table 4, Table 4, Table 4. Nusev et al. (2021) G. Nusev, B. Morel, J. Mougin, Ð. Juričić, and P. Boškoski Condition monitoring of solid oxide fuel cells by fast electrochemical impedance spectroscopy: a case example of detecting deficiencies in fuel supply. Journal of Power Sources 489, p. 229491. External Links: ISSN 0378-7753, Document, Link Cited by: §3.3, Table 2, Table 4. Obregon et al. (2023) J. Obregon, Y. Han, C. W. Ho, D. Mouraliraman, C. W. Lee, and J. Jung Convolutional autoencoder-based soh estimation of lithium-ion batteries using electrochemical impedance spectroscopy. Journal of Energy Storage 60, p. 106680. External Links: ISSN 2352-152X, Document Cited by: §1. Odena et al. (2016) A. Odena, V. Dumoulin, and C. Olah Deconvolution and checkerboard artifacts. Distill. External Links: Link, Document Cited by: §3.1. Oppenheim and Schafer (2013) A. Oppenheim and R. Schafer Discrete-Time Signal Processing . Pearson Deutschland. External Links: ISBN 9781292025728, Document, Link Cited by: §2. O’Shea and Nash (2015) K. O’Shea and R. Nash An introduction to convolutional neural networks. External Links: 1511.08458, Link Cited by: §3. Plank et al. (2024) C. Plank, T. Rüther, L. Jahn, M. Schamel, J. P. Schmidt, F. Ciucci, and M. A. Danzer A review on the distribution of relaxation times analysis: a powerful tool for process identification of electrochemical systems. Journal of Power Sources 594, p. 233845. External Links: ISSN 0378-7753, Link, Document Cited by: §1. Raissi et al. (2019) M. Raissi, P. Perdikaris, and G.E. Karniadakis Physics-informed neural networks: a deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations. Journal of Computational Physics 378, p. 686–707. External Links: ISSN 0021-9991, Document Cited by: §1. Saccoccio et al. (2014) M. Saccoccio, T. H. Wan, C. Chen, and F. Ciucci Optimal regularization in distribution of relaxation times applied to electrochemical impedance spectroscopy: ridge and lasso regression methods-a theoretical and experimental study. Electrochimica Acta 147, p. 470–482. Cited by: §1. Wan et al. (2015) T. H. Wan, M. Saccoccio, C. Chen, and F. Ciucci Influence of the discretization methods on the distribution of relaxation times deconvolution: implementing radial basis functions with drttools. Electrochimica Acta 184, p. 483–499. External Links: ISSN 0013-4686, Document Cited by: Figure 6. Wanhammar and Saramaki (2020) L. Wanhammar and T. Saramaki Digital filters using matlab. External Links: ISBN 978-3-030-24062-2, Document Cited by: §2, §2. Wiener (1949) N. Wiener Extrapolation, interpolation, and smoothing of stationary time series: with engineering applications. The MIT Press. External Links: ISBN 9780262257190, Document Cited by: §2. Žnidarič et al. (2021) L. Žnidarič, G. Nusev, B. Morel, J. Mougin, Đ. Juričić, and P. Boškoski Evaluating uncertainties in electrochemical impedance spectra of solid oxide fuel cells. Applied Energy 298, p. 117101. External Links: ISSN 0306-2619, Document, Link Cited by: §1.