---
title: 'SeqBattNet: Adaptive Battery Voltage Forecasting'
url: https://www.emergentmind.com/topics/seqbattnet
type: topic
---

# SeqBattNet: Adaptive Battery Voltage Forecasting

SeqBattNet is a discrete-state physics-informed neural network for battery modeling that predicts the terminal-voltage trajectory of a lithium-ion battery during discharge from the applied current profile. It is presented as a battery-management-system submodel for forecasting quantities such as remaining discharge time and remaining discharge energy, both of which depend on accurately predicting when the voltage will cross the end-of-discharge cutoff. The model combines a cycle-adaptive encoder with an equivalent-circuit-model decoder, requires only three basic battery parameters—\(V_0\), \(V_{\mathrm{EOD}}\), and \(C\)—and is designed to remain computationally efficient while generalizing across cells even when trained on data from a single cell [2509.17621].

## 1. Problem setting and design rationale

SeqBattNet addresses the specific task of forecasting terminal voltage during the discharge process. In the formulation used by the paper, this task is not treated as an isolated regression problem: voltage prediction is the mechanism by which downstream battery-management quantities are obtained, especially remaining discharge time (RDT) and remaining discharge energy (RDE). If the predicted voltage trajectory is inaccurate, the predicted cutoff time is also inaccurate, which directly degrades time-to-empty and remaining-energy estimates [2509.17621].

The model is motivated by limitations in three established method families. Classical model-based approaches, including electrochemical models and equivalent circuit models (ECMs), are physically grounded but require many parameters that are difficult to identify and calibrate. High-fidelity electrochemical models such as P2D/DFN or SPM are especially costly for real-time use. Purely data-driven methods avoid explicit physics but require large labeled datasets spanning many operating conditions and degradation states, and the paper characterizes them as less interpretable. Prior physics-informed neural networks partially bridge these two regimes, but the paper identifies three remaining shortcomings: some depend on continuous-time ODE state generation, some still require additional physical parameters such as \(R_p\), \(R_s\), \(C_{sp}\), and \(C_s\), and some do not explicitly adapt to aging [2509.17621].

Against this background, SeqBattNet is defined by three coupled design choices. It is **physics-informed** because its decoder is structured as a discrete-time ECM with explicit state transitions and an output equation following battery circuit dynamics. It is **discrete-state** because internal states are propagated step by step in discrete time rather than regenerated by integrating an ODE over a preset horizon. It is **aging-adaptive** because the model infers cycle-specific latent physical parameters from an initial observed current-voltage segment at the start of each discharge cycle. The paper claims four main advantages: only three basic battery parameters are required; battery evolution is modeled through discrete physical states; the method is robust even when trained on data from only a single cell; and it remains computationally efficient [2509.17621].

A recurrent point in the paper is that discrete-state rollout is operationally aligned with battery-management deployment. The decoder can advance one sample at a time and stop as soon as the predicted voltage falls below the cutoff, rather than rolling continuously to an externally chosen terminal time. This is presented as especially relevant for proactive stopping in RDT and RDE estimation [2509.17621].

## 2. Encoder-decoder architecture and aging adaptation

SeqBattNet is an encoder-decoder model. For each discharge cycle, the encoder consumes the first \(n\) time steps of the current and voltage sequences and outputs cycle-specific adaptation parameters. The decoder then uses those adaptation parameters together with the subsequent current sequence to roll forward the voltage trajectory [2509.17621].

The encoder input is
\[
\mathbf{x} = [I_{\mathrm{init}}, V_{\mathrm{init}}] \in \mathbb{R}^{n \times 2},
\]
where \(I_{\mathrm{init}}\) and \(V_{\mathrm{init}}\) are the first \(n\) current and voltage samples of the cycle. This sequence is embedded by a linear layer,
\[
\mathbf{e} = \mathbf{W}_{\mathrm{emb}}\mathbf{x} + \mathbf{b}_{\mathrm{emb}} \in \mathbb{R}^{n \times d_{\mathrm{emb}}}.
\]

The encoder is the paper’s proposed **HRM-GRU** module, where HRM denotes Hierarchical Reasoning Model. The architecture adapts a hierarchical-reasoning design that originally used Transformers, but replaces them with GRUs to reduce computational complexity and, according to the reported experiments, improve performance. HRM-GRU maintains low-level hidden states \(\mathbf{z}^L_s\) and high-level hidden states \(\mathbf{z}^H_s\), and introduces gradient-free preparatory micro-updates before each main update. For each sequence index \(s\), it performs \(Q = NT - 1\) micro-updates,
\[
(\mathbf{z}^L_{s^-}, \mathbf{z}^H_{s^-}) = \mathcal{U}_{N,T}\!\big(\mathbf{z}^L_{s-1}, \mathbf{z}^H_{s-1}; \mathbf{e}_s\big),
\]
followed by the main update,
\[
\mathbf{z}^L_s = \mathrm{GRU}_L\!\big([\mathbf{e}_s,\mathbf{z}^H_{s^-}], \mathbf{z}^L_{s^-}\big),
\qquad
\mathbf{z}^H_s = \mathrm{GRU}_H\!\big(\mathbf{z}^L_s, \mathbf{z}^H_{s^-}\big).
\]
This two-level low-to-high recurrence is intended to capture more structured temporal dependencies than a one-level GRU [2509.17621].

After dropout and layer normalization of the high-level state, a linear head produces
\[
\mathbf{o}_s = \mathbf{W}_{\mathrm{out}}\,\mathrm{Norm}(\mathbf{z}^H_s) + \mathbf{b}_{\mathrm{out}} \in \mathbb{R}^{d_{\mathrm{out}}},
\]
with
\[
d_{\mathrm{out}} = 2e_{\mathrm{rc}} + 3.
\]
Using the final output \(\mathbf{o}_n\), the encoder predicts
\[
\big[R_0,\boldsymbol{\tau},\mathrm{SOC}_0,\mathrm{SOH},\mathbf{w}\big] = \mathbf{o}_n,
\]
where \(R_0\) is the ohmic resistance, \(\boldsymbol{\tau} = \{\tau_1,\dots,\tau_{e_{\mathrm{rc}}}\}\) are RC branch time constants, \(\mathrm{SOC}_0\) is the initial state of charge, \(\mathrm{SOH}\) is the state of health, and \(\mathbf{w}\) are weights used to initialize the RC branch voltages. The number of RC branches is fixed to \(e_{\mathrm{rc}} = 2\) [2509.17621].

The raw outputs are mapped into physically plausible ranges:
\[
R_0 \mapsto \mathrm{Aff}\text{-}\sigma(R_0;10^{-3},0.5),
\]
\[
\boldsymbol{\tau} \mapsto \mathrm{Aff}\text{-}\sigma(\boldsymbol{\tau};10^{-2},10^{5}),
\]
\[
\mathrm{SOC}_0 \mapsto \sigma(\mathrm{SOC}_0), \qquad
\mathrm{SOH} \mapsto \sigma(\mathrm{SOH}), \qquad
\mathbf{w} \mapsto \mathrm{Softmax}(\mathbf{w}),
\]
with
\[
\mathrm{Aff}\text{-}\sigma(z;\ell,h)=\ell+(h-\ell)\sigma(z).
\]
The paper presents these mappings as a key physics-consistency mechanism because they enforce positivity, boundedness, and normalized RC initialization weights by construction [2509.17621].

The encoder also computes the initial RC branch voltages:
\[
\mathrm{OCV}_0 = V_{\mathrm{EOD}} + (V_0 - V_{\mathrm{EOD}})\,\sigma(g(\mathrm{SOC}_0)),
\]
\[
s_0 = \mathrm{OCV}_0 - R_0 I_{\text{last}} - V_{\text{last}},
\]
\[
\mathbf{v}_{RC,0} = \mathbf{w}\, s_0 \in \mathbb{R}^{e_{\mathrm{rc}}}.
\]
Here \(I_{\text{last}}\) and \(V_{\text{last}}\) are the final current and voltage samples in the initial observed segment, and \(g(\cdot)\) is the decoder’s feedforward network mapping SOC to OCV. The paper interprets \(s_0\) as the residual between OCV and measured loaded voltage after subtracting the ohmic drop, with \(\mathbf{w}\) partitioning that residual across the RC branches [2509.17621].

The aging-adaptation mechanism operates once per discharge cycle rather than at every decoder step. For each cycle, the first \(n\) measured current-voltage samples yield one set of adaptation parameters,
\[
R_0,\;\boldsymbol{\tau},\;\mathrm{SOC}_0,\;\mathrm{SOH},\;\mathbf{v}_{RC,0},
\]
which are then held fixed during decoder rollout, except for decoder-internal quantities such as \(\mathrm{SOC}_m\), \(\mathbf{r}_{RC,m}\), and \(\mathbf{v}_{RC,m}\), which evolve over time. The paper provides qualitative support for this mechanism using PCA of encoder embeddings that form smooth trajectories across cycle index [2509.17621].

## 3. Discrete-time ECM decoder and physics-informed formulation

The decoder is a discrete-time, physics-informed ECM that uses the cycle-adaptive parameters
\[
R_0,\;\boldsymbol{\tau},\;\mathrm{SOC}_0,\;\mathrm{SOH},\;\mathbf{v}_{RC,0}
\]
together with the future current sequence \(\mathbf{I} = \{I_{n+1},\dots,I_{t_{\mathrm{EOD}}}\}\) to generate future terminal voltage [2509.17621].

The RC decay coefficients are
\[
\alpha = \exp(-\Delta t/\boldsymbol{\tau}),
\]
with the exponential applied elementwise over the RC branches. At each step \(m=n+1,\dots,t_{\mathrm{EOD}}\), the open-circuit voltage is predicted from the current SOC:
\[
\mathrm{OCV}_m = V_{\mathrm{EOD}} + (V_0 - V_{\mathrm{EOD}})\,\sigma(g(\mathrm{SOC}_m)).
\]
This constrains OCV to the interval between cutoff voltage \(V_{\mathrm{EOD}}\) and nominal full-charge voltage \(V_0\) [2509.17621].

The dynamic RC branch resistances are generated by a second neural network:
\[
\mathbf{r}_{RC,m} = r_{\min} + (r_{\max}-r_{\min})\,\sigma\!\big(f(\mathrm{SOC}_m,\mathrm{SOH})\big),
\]
with \(r_{\min}=10^{-4}\) and \(r_{\max}=1\). In this formulation, polarization resistance depends on both the present charge state and health state [2509.17621].

The RC branch voltages obey the discrete first-order relaxation equation
\[
\mathbf{v}_{RC,m} = \alpha \mathbf{v}_{RC,m-1} + (1-\alpha)\,\mathbf{r}_{RC,m} I_m.
\]
This is the central state transition of the model. Because the decoder explicitly carries forward \(\mathbf{v}_{RC,m-1}\), the state is not regenerated from scratch; it evolves causally in discrete time [2509.17621].

SOC is propagated by Coulomb counting:
\[
\mathrm{SOC}_{m+1} = \Pi_{[0,1]}\!\left(\mathrm{SOC}_m - \frac{I_m \Delta t}{3600\,C_{\mathrm{eff}}}\right),
\]
where \(\Pi_{[0,1]}\) clips to \([0,1]\), and \(C_{\mathrm{eff}}\) is the degradation-adjusted effective capacity,
\[
C_{\mathrm{eff}} = \beta C_{\mathrm{EOL}} + \big(C_{\mathrm{rated}} - \beta C_{\mathrm{EOL}}\big)\,\mathrm{SOH}.
\]
The paper states that \(C_{\mathrm{EOL}}\) is typically \(0.8\,C_{\mathrm{rated}}\), and uses \(\beta=0.8\) “to account for capacity fade.” For the NASA dataset, \(C_{\mathrm{EOL}}=0.6\,C_{\mathrm{rated}}\) is used instead [2509.17621].

The ECM output equation is
\[
\widehat{V}_m = \mathrm{OCV}_m - R_0 I_m - \sum_{k=1}^{e_{\mathrm{rc}}} v_{RC,m,k}.
\]
Thus the predicted terminal voltage is the equilibrium voltage minus an instantaneous ohmic drop minus the polarization drops from the RC branches. The complete predicted sequence is written as
\[
\widehat{\mathbf{V}} = [V_0,\dots,V_n,\widehat{V}_{n+1},\dots,\widehat{V}_{t_{\mathrm{EOD}}}],
\]
where the measured early segment is concatenated with the predicted future segment [2509.17621].

The paper uses this decoder structure to distinguish SeqBattNet from unconstrained sequence models. The model does not directly map \((I,V)\)-history to future voltage with a black-box neural network. Instead, it predicts physically interpretable latent quantities—\(R_0\), \(\tau_k\), SOC, SOH, RC voltages, OCV, and RC resistances—and evolves them through explicit state equations. Neural subnetworks are inserted only where nonlinear constitutive relations are learned: OCV as a function of SOC via \(g(\cdot)\), and RC resistance as a function of SOC and SOH via \(f(\cdot)\) [2509.17621].

A related conceptual point is the paper’s distinction between **physics-informed** and merely **physics-inspired** modeling. In SeqBattNet, the physical structure is not a post hoc regularizer but part of the forward model itself. The exact effect of removing the ECM decoder is not quantified by a direct decoder-ablation table, but the reported comparisons against plain LSTM, GRU, and Transformer sequence models are used to argue that the combination of learned adaptation and ECM structure yields materially lower error on TRI and RT-Batt [2509.17621].

## 4. Training objective, implementation, and battery metadata requirements

SeqBattNet is trained using a weighted Huber loss over the voltage trajectory:
\[
\mathcal{L}_{\mathrm{wL1}} =
\frac{1}{\sum_{i=1}^{t_{\mathrm{EOD}}} a_i w_i}
\sum_{i=1}^{t_{\mathrm{EOD}}} a_i w_i \,\ell_\beta(\hat V_i,V_i),
\]
where \(a_i \in \{0,1\}\) is a mask for zero-padded sequences and \(w_i\) is a time-dependent weight. The Huber loss is
\[
\ell_{\beta}(\hat{V}_i,V_i)=
\begin{cases}
\frac{1}{2\beta}(V_i-\hat{V}_i)^2, & |V_i-\hat{V}_i|<\beta,\\[4pt]
|V_i-\hat{V}_i|-\tfrac{1}{2}\beta, & \text{otherwise},
\end{cases}
\]
with \(\beta=0.1\). The adaptive temporal weights are designed to emphasize both the beginning of the trajectory and especially the final prediction step near end of discharge [2509.17621].

The loss-function study reports that this weighted Huber loss outperforms both plain MSE and the baseline loss from Wang et al. On TRI-b1c4, RMSE drops from \(0.1299\) with MSE and \(0.0624\) with the baseline loss to \(0.0280\) with the proposed loss. On RT-Batt-1-2, RMSE improves from \(0.0902\) and \(0.0579\) to \(0.0346\). On NASA-RW5, RMSE improves from \(0.0329\) and \(0.0301\) to \(0.0253\). MAE and MAPE are reported to follow the same trend [2509.17621].

The model is implemented in PyTorch and optimized with AdamW. Training uses 1000 epochs, batch size 128, and an NVIDIA 4080 GPU. The learning rate starts at \(2\times 10^{-3}\) and is reduced by ReduceLROnPlateau when validation loss stagnates, with a floor at \(10^{-4}\). Each experiment is repeated five times, and final predictions are averaged across runs. The HRM-GRU hyperparameters selected by grid search are \(N=1\), \(T=2\), \(d_{\mathrm{emb}}=32\), \(d_L=128\), and \(d_H=64\). The OCV network \(g(\cdot)\) is a 3-layer FNN with dimensions \(1\to 32 \to 32 \to 1\), using SiLU, then SiLU, then Sigmoid. The RC-resistance network \(f(\cdot)\) is \(2\to 32 \to e_{\mathrm{rc}}\), using SiLU and then Sigmoid [2509.17621].

A central deployment claim is that SeqBattNet requires only three basic battery parameters:
\[
V_0,\qquad V_{\mathrm{EOD}},\qquad C.
\]
More concretely, the capacity term is the rated capacity \(C_{\mathrm{rated}}\), together with an end-of-life capacity setting for each dataset. \(V_0\) sets the upper limit of the OCV range, \(V_{\mathrm{EOD}}\) sets the lower limit and the discharge stopping criterion, and capacity provides the scale for SOC propagation by Coulomb counting. The paper explicitly emphasizes that users do not need to pre-identify additional ECM parameters such as \(R_p\), \(R_s\), \(C_{sp}\), or \(C_s\); instead, \(R_0\), \(\tau_k\), and RC-resistance profiles are inferred from data in a constrained way [2509.17621].

## 5. Datasets, evaluation protocol, and empirical results

The evaluation is deliberately organized around cross-cell generalization. On each dataset, training is performed on one cell, validation on a different cell, and testing on different cells again. The paper presents this as stricter than splitting cycles within a single cell, and explicitly criticizes Wang et al.’s BattNN protocol for splitting within a cell because training and test fragments then share the same cell and cycle-trajectory context [2509.17621].

| Dataset | Cell split | Battery settings |
|---|---|---|
| TRI | train b1c0; val b2c13; test b1c4, b1c2 | \(C_{\mathrm{rated}}=1.1\) Ah, \(C_{\mathrm{EOL}}=0.88\) Ah, \(V_0=3.6\) V, \(V_{\mathrm{EOD}}=2.0\) V |
| RT-Batt | train 1-1; val 1-4; test 1-2, 1-3 | \(C_{\mathrm{rated}}=1.1\) Ah, \(C_{\mathrm{EOL}}=0.88\) Ah, \(V_0=3.6\) V, \(V_{\mathrm{EOD}}=2.0\) V |
| NASA | train RW3; val RW6; test RW4, RW5 | \(C_{\mathrm{rated}}=2.22\) Ah, \(C_{\mathrm{EOL}}=1.33\) Ah, \(V_0=4.2\) V, \(V_{\mathrm{EOD}}=3.2\) V |

The three datasets differ in discharge-profile variability. TRI contains 124 LFP/graphite A123 APR18650M1A cells with highly standardized constant 4C discharge down to \(2\) V, while charge protocols vary. RT-Batt also uses A123 APR18650M1A LFP/graphite cells, but each cell has its own multistage discharge profile that remains fixed across that cell’s cycles. NASA contains four 18650 lithium-ion cells RW3–RW6 from the NASA Ames repository, charged at constant 2 A to \(4.2\) V and discharged to \(3.2\) V with randomized current profiles between 0.5 A and 4 A; both cells and cycles vary in discharge profile. The initial adaptation window length is chosen empirically: \(n=80\) works best on TRI and RT-Batt, while \(n=30\) works best on NASA [2509.17621].

The main baseline set includes BattNN, plain LSTM, plain GRU, plain Transformer, and multiple SeqBattNet variants that keep the same decoder and change only the encoder: LSTM, GRU, Transformer, FNN, HRM-LSTM, HRM-GRU, and HRM-Transformer. This design isolates the contribution of the aging encoder from that of the physics-informed decoder [2509.17621].

Relative to BattNN, classical sequence models already improve performance substantially. On TRI-b1c4, BattNN has RMSE \(0.5781\), whereas LSTM and GRU are around \(0.071\). On RT-Batt-1-2, BattNN has RMSE \(0.5536\), while LSTM achieves \(0.1303\) and Transformer \(0.1450\). On NASA-RW4, BattNN has RMSE \(0.0502\), while Transformer reaches \(0.0251\) [2509.17621].

Within the SeqBattNet family, the proposed HRM-GRU encoder is strongest on TRI and RT-Batt, and second-best on NASA, where the simple FNN encoder performs slightly better. The paper explicitly identifies this as dataset-specific encoder suitability rather than a universal dominance claim [2509.17621].

| Test case | HRM-GRU RMSE | HRM-GRU MAE / MAPE |
|---|---:|---:|
| TRI-b1c4 | 0.0280 | 0.0113 / 0.0042 |
| TRI-b1c2 | 0.0416 | 0.0168 / 0.0063 |
| RT-Batt-1-2 | 0.0346 | 0.0120 / 0.0044 |
| RT-Batt-1-3 | 0.0491 | 0.0171 / 0.0063 |
| NASA-RW4 | 0.0257 | 0.0191 / 0.0055 |
| NASA-RW5 | 0.0253 | 0.0191 / 0.0055 |

Two ablation trends are especially clear. First, hierarchical recurrent encoding matters on TRI and RT-Batt: on TRI-b1c4, plain SeqBattNet with a GRU encoder gives RMSE \(0.0742\), while HRM-GRU gives \(0.0280\); on RT-Batt-1-2, plain SeqBattNet with GRU gives \(0.0827\), while HRM-GRU gives \(0.0346\). Second, the specific hierarchical design matters: HRM-GRU is both more accurate and faster than HRM-Transformer on TRI and RT-Batt. On TRI-b1c4, HRM-Transformer gives RMSE \(0.0701\) and takes \(109.5\) s, whereas HRM-GRU gives \(0.0280\) and takes \(68.8\) s [2509.17621].

The paper also reports favorable runtime. SeqBattNet-HRM-GRU takes roughly \(46\)–\(85\) s per experiment on TRI and RT-Batt and under 1 s on NASA in the reported setup, while remaining competitive with or faster than several sequence baselines and much faster than HRM-Transformer. The model does not provide a total neural parameter count, so its “low parameterization” is explicitly framed as low **physical-parameter requirements** rather than a quantified network-size comparison [2509.17621].

## 6. Interpretation, scope, and limitations

The paper’s interpretation of SeqBattNet’s performance is that it balances inductive bias and adaptability. The encoder extracts degradation-sensitive information from a short initial discharge fragment, while the decoder constrains prediction through explicit battery physics: SOC decreases by Coulomb counting, RC voltages relax exponentially, terminal voltage is assembled from OCV minus ohmic and polarization drops, and latent variables are confined to physically valid ranges. This suggests that the model avoids spending representational capacity on rediscovering known structure while retaining flexibility where constitutive mappings are nonlinear [2509.17621].

Several clarifications follow from the reported evidence. First, aging adaptation is inferred once per sequence or cycle, not continuously at every decoder step. Second, the paper argues for the value of the discrete-state design conceptually—particularly proactive stopping at cutoff and causal state carryover—but does not report a numerical ablation against a continuous-state version. Third, qualitative PCA evidence is provided for aging-aware embeddings, but the physics-informed decoder is not isolated in a direct “decoder removed” ablation table [2509.17621].

The paper also states its limitations clearly. The model uses only voltage and current as inputs and explicitly ignores temperature. All datasets come from controlled laboratory conditions with minimal temperature variation, so the method has not been validated under realistic environmental variation where temperature influences internal resistance, capacity fade, and voltage response. The work focuses on discharge rather than charge dynamics, and its chemistry and operating-condition coverage is limited to the benchmark datasets used: mainly A123 APR18650M1A LFP/graphite cells in TRI and RT-Batt, and the specific NASA randomized-use cells. Future work is identified as incorporating ambient and cell temperature as additional inputs and validating on real-world datasets with more diverse environmental and operational conditions [2509.17621].

In the form presented, SeqBattNet occupies a specific position in battery modeling. It is more constrained and physically structured than a black-box sequence model, yet less dependent on hand-calibrated physical parameter sets than traditional ECM or electrochemical modeling. Its principal contribution is a cycle-adaptive, discrete-time voltage predictor that embeds ECM state evolution directly into the prediction loop, infers aging-sensitive parameters from a short initial discharge segment, and demonstrates low RMSE with single-cell training across TRI, RT-Batt, and NASA benchmarks [2509.17621].

Source: https://www.emergentmind.com/topics/seqbattnet